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This  study  develops  a  model  for  the  source  of  radio 
recombination  line  emission  from  the  galactic  plane.   The 
model  is  based  mainly  on  recombination  line  spectra  ob- 
served in  the  direction  of  the  supernova  remnant  3C391  . 
The  H157a  recombination  line  was  observed  and  used  with  the 
H92a  spectrum  observed  by  Cesarsky  and  Cesarsky  (1973)  to 
obtain  a  model  consistent  with  both  spectra.   Other 
criteria— such  as  the  80  MHz  optical  depth,  the  contri- 
bution to  the  thermal  continuum  temperature,  and  ionization 
rates  —  were  also  considered. 

Four  models  were  studied:   a  cold  region  illumi- 
nated by  3C391  ;  an  isolated  cold  region;  a  group  of  hot, 
dense  H  II  regions;  and  a  low  density  H  II  nebula.   Of 
these,  the  best  model  is  a  small  (^85  pc  diameter)  H  II 

region  with  a  temperature  of  7500°  K  and  an  average  elec- 

-3    -3 
tron  density  of  10   cm 
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The  H157a  spectrum  in  the  direction  of  the  pulsar 
1858+03  was  also  observed.   It  was  found  that  it  is  prob- 
able that  one  of  the  recombination  line  nebulae  lies  along 
the  line  of  sight  to  the  pulsar.   In  this  case  the  pulsar 
would  be  at  a  distance  of  2.6  kpc  in  order  for  its  dis- 
persion measure  to  agree  with  the  measured  value  of  402 

-3 

pc  cm 

The  large  scale  distribution  of  the  H  II  nebulae 
was  also  studied.   The  results  of  available  radio  recom- 
bination line  spectra  surveys  from  the  galactic  plane  sug- 
gest  that  they  exist  within  an  annular  region,  which  ex- 
tends from  4  to  about  8.7  kpc  from  the  galactic  center. 
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CHAPTER  I 
INTRODUCTION 

The  Galactic  Ridge 

As  a  radio  telescope  is  scanned  across  the  plane 
of  the  Milky  Way  at  some  galactic  longitude,   it  is  seen 
to  be  a  source  of  radio  emission.   For  radio  frequencies 
of  100  MHz  or  greater,  the  intensity  increases  as  the 
galactic  equator  is  approached  and  peaks  at  the  equator 
itself.   Scans  through  the  galactic  plane  made  at  several 
different  frequencies  reveal  that  the  emission  has  two 
components,  thermal  and  nonthermal.   At  frequencies  below 
a  few  GHz,  the  emission  is  predominantly  nonthermal  with 
a  spectral  index  between  -2.7  and  -3.0  (Westerhout,  1958; 
Matthews,  Pedlar  and  Davies,  1973).   At  85  MHz  the  radi- 
ation is  almost  completely  nonthermal  and  the  intensity 


In  this  work,  the  revised  galactic  latitude,  b 
and  longitude,  l11,  as  adopted  by  the  IAU  (Blaauw  et  al  .  , 
I960)  will  be  used  throughout.   The  superscript  II,  how- 
ever, will  be  dropped  from  the  notation. 


distribution  across  the  galactic  plane  has  a  total  half-in- 
tensity width  of  about  5°  at  longitude  greater  than  12° 
from  either  side  of  the  galactic  center  (Hill,  Slee,  and 
Mills, 1958).   As  the  frequency  increases,  the  overall  in- 
tensity decreases  and  the  width  of  the  emission  becomes 
smaller,  approaching  a  thin  ridge  at  centimeter  wavelengths. 
At  1390  MHz  the  hal f - i ntensi ty  width  is  generally  about 
195.   The  thickness  of  the  nonthermal  emission  region  in- 
creases somewhat  near  the  galactic  center.   This  thin  con- 
tinuum emission  region  surrounding  the  galactic  equator 
has  been  named  the  "galactic  ridge." 

The  galactic  ridge  has  been  recently  mapped  with 
an  antenna  resolution  of  11  arc  minutes  at  frequencies  of 
1414,  2695,  and  5000  MHz  by  Altenhoff  et  aj_.  (1970). 
Matthews,  Pedlar,  and  Davies  (1973)  have  used  these  maps 
to  obtain  the  latitude  distribution  of  the  continuum  in- 
tensities at  three  longitudes  where  no  resolved  sources 
are  present.   They  plotted  the  distribution  with  the  in- 
tensity scale  at  each  of  the  frequencies  adjusted  so  that, 
if  the  emission  were  entirely  nonthermal,  the  three  dis- 
tribution curves  would  be  congruent.   The  plots  are  con- 
gruent at  latitudes  greater  than  about  0?5,  but  inside  this 
range,  the  intensity  increases  markedly  with  frequency  due 
to  thermal  emission.   These  curves  illustrate  the  two  com- 
ponents of  the  galactic  ridge  continuum  emission. 


It  is  known  that  hydrogen  is  concentrated  within 
a  thin  region  around  the  galactic  equator  (cf.  Westerhout, 
1958;  Baker  and  Burton,  1 975)  .   The   thermal  emission  comes 
from  the  free-free  process  between  electrons  and  protons 
from  ionized  hydrogen. 

At  frequencies  below  about  30  MHz,  a  scan  cutting 
through  the  ridge  shows  an  intensity  profile  with  two 
peaks,  one  on  each  side  of  the  equator,   The  decrease  of 
intensity  along  the  equator  results  from  absorption  of 
the  strong  nonthermal  emission  by  the  free-free  process 
within  the  thermal  emission  region.   The  map  by  Shain, 
Komesaroff,  and  Higgins  (1961)  at  19.7  MHz  shows  the  ab- 
sorption lane  prominently. 

The  galactic  ridge  took  on  an  interesting  aspect 
when  Gottesman  and  Gordon  (1970)  observed  low  intensity 
H157a  radio  recombination  line  spectra  in  two  directions 
in  which  the  Altenhoff  et^  a]_.  (1970)  maps  show  no  resolved 
sources.   Gordon,  Brown,  and  Gottesman  (1972)  found  that 
at  a  longitude  of  33°,  the  recombination  line  emission  is 
limited  to  a  latitude  of  less  than  0?5  above  the  equator. 
In  the  region  toward  the  galactic  center,  recombination 
line  emission  has  been  observed  at  a  latitude  of  -0.62 
(Lockman  and  Gordon, 1 973) .   Previous  to  Gottesman  and 
Gordon's  observations,  radio  recombination  lines  had  only 
been  observed  from  known  galactic  H  II  regions.   They  have 
not  been  detected  from  the  interstellar  medium  in  the  solar 
neighborhood . 


A  number  of  workers  have  studied  the  galactic 
ridge  in  order  to  determine  the  nature  of  the  radio  recom- 
bination line  sources  and  how  they  relate  to  the  thermal 
continuum  and  21  cm  spectral  line  emission.   The  analysis 
of  Cesarsky  and  Cesarsky  (1971,  1973b)  and  Lockman  and 
Gordon  (1973)  find  that  the  electron  temperature  of  the 

radio  recombination  line  emission  regions  is  close  to  20  °K. 

3      4 
Other  workers   derive  temperatures  ranging  from  10   to  10 

°K  (Gottesman  and  Gordon,  1 970 ;  Gordon  and  Gottesman, 1 971  ; 
Matthews,  Pedlar,  and  Davies,  1973;  Jackson  and  Kerr,  1971; 
Chaisson, 1 974a;  and  Downes  and  Wi 1  son,  1 974)  .   Very  recently 
two  papers  have  appeared,  whi ch  conclude  that  these  recom- 
bination lines  originate  from  low  density  H  II  regions 
(Jackson  and  Kerr,  1975;  Pankoni  n,  1  975  )  . 


Objective,  Method,  and  Scope  of  This  Study 

The  problem  to  which  this  study  is  directed  is  that 
of  determining  the  type  of  region  responsible  for  the  recom- 
bination line  emission  in  the  galactic  ridge.   The  first 
step  will  be  to  observe  recombination  line  emission  from 
regions  in  the  ridge  in  which  there  are  no  H  II  regions 
visible  on  the  Altenhoff  et  al_.  (1970)  maps  and  from  which 
other  data  have  previously  been  obtained.   The  recombination 
line  spectrum  will  provide  one  piece  of  datum,  the  power 
contained  in  the  spectral  line.   However,  there  are  several 


parameters  which  define  the  type  of  region:   electron  tem- 
perature, electron,  neutral  hydrogen,  ionized  hydrogen 
number  densities,  size  of  region,  and  distribution  of  the 
region,  or  regions,  within  the  galactic  ridge.   The  latter 
may  imply  a  relation  between  recombination  line  emission 
and  other  phenomena  such  as  low  frequency  radio  absorption 
and  galactic  structure. 

Since  the  observations  made  as  a  part  of  this  study 
can  only  supply  one  bit  of  information,  the  identity  of 
the  recombination  line  source  can  be  obtained  only  after 
other  observations  are  utilized.   Even  so,  there  are  still 
more  parameters  whose  values  are  to  be  found  than  there 
are  observed  quantities.   For  that  reason,  the  values  of 
some  parameters  must  be  chosen  at  the  onset  of  model  com- 
putation.  A  type  of  model  is  hypothesized  and  theory  or 
the  conclusions  from  previous  observation  suggest  what 
these  values  are.   The  second  step  will  be  to  obtain,  from 
one  region  in  the  ridge,  a  model  solution  which  is  consis- 
tent between  two  similar  data.   The  ideal  case  for  this 
is  the  region  toward  the  supernova  remnant  3C391  from  which 
Cesarsky  and  Cesarsky  (1973b)  have  observed  an   H92a  recom- 
bination line  spectrum.   This  region  has  undergone  a  con- 
siderable amount  of  scrutiny,  particularly  by  those  who 
have  attempted  to  measure  the  80  MHz  optical  depth  (Cas- 
well et  al.,  1971;  Dulk  and  Slee,  1972;  Goss,  1972). 


The  third  step  in  the  process  will  be  to  assume 
that  this  model  applies  in  the  mean  sense  to  the  entire 
ridge  and  use  the  remaining  data  obtained  for  this  study 
to  check  the  model.   The  method  is  to  calculate  a  parame- 
ter which  depends  on  the  model,  to  see  if  it  agrees  with 
what  is  observed.   Finally,  recombination  line  surveys 
will  be  combined  in  such  a  way  as  to  examine  the  large 
scale  structure  of  the  recombination  line  ridge. 

The  theoretical  basis  of  this  study  will  be  devel- 
oped in  Chapter  II.   All  equations  used  in  the  analysis 
will  be  derived  and  approximations  and  assumptions  will  be 
discussed.   The  third  chapter  will  discuss  theory  behind 
the  observing  procedure.   The  method  and  conditions  under 
which  the  observations  were  made  and  the  display  of  the 
processed  spectra  will  occupy  Chapter  IV.   Chapters  II, 
III,  and  IV  provide  a  basis  for  the  subsequent  work. 

Four  types  of  regions  on  which  to  base  the  models 
will  be  selected.   Three  of  them  result  from  previous  obser- 
vation of  3C391,  a  summary  of  which  begins  Chapter  V.   The 
Hot  Spot  model  is  based  on  unpublished  observation  of  Alten- 
hoff  et  aj.  (1973).   The  remainder  of  Chapter  V  develops 
the  four  models.   Chapter  VI  is  involved  with  the  selection 
of  the  most  probable  model  and  tests  this  by  computing  the 
distance  to  the  pulsar  1858+03  and  the  80  MHz  optical  depth 
expected  in  the  direction  of  3C396.   Chapter  VI  ends  with 


an  investigation  of  the  large  scale  structure  of  ridge 

recombination  line  emission.  The  last  chapter  concludes 

the  study  by  proposing  future  projects  to  study  the  galac- 
tic ridge. 


CHAPTER  II 
INTERSTELLAR  MEDIUM  ASTROPHYSICS 

Introduction 


It  is  the  purpose  of  this  chapter  to  present  a 
foundation  on  which  this  project  is  based,  particularly 
the  analysis  of  the  spectra.   The  starting  point  is  the 
equation  of  radiative  transfer.   The  coefficients  of  ab- 
sorption and  emission  are  parameters  of  the  solution  to 
the  transfer  equation. 

Since  this  study  is  based  on  radio  spectral  lines, 
the  line  absorption  coefficient  will  be  examined  in  de- 
tail.  The  line  profile  of  a  given  element  depends  on  the 
thermodynamic  state,  temperature,  abundance,  electron 
density,  and  motion  of  the  gas  creating  the  line. 

The  Transfer  Equation  and  a  General  Solution  of  It 

The  transfer  equation  can  be  set  up  by  considering 
the  radiation  passing  through  a  region  as  shown  in  Figure 
II-l.   The  radiation  entering  this  region  has  a  frequency- 
dependent  specific  intensity,  I  (v),  which  can  be  both 
thermal  and  nonthermal  in  origin.   The  line  of  sight  enters 


ds 


I(v) 


I0(y) 


OBSERVER 


s=L 


s=0 
r=0 


Figure  II-l.   Ray  Path  Geometry  Through  a  Line- 
emitting  Region.   The  region  is 
shown  shaded. 
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the  region  at  a  point  where  the  path  length,  s,  is  zero. 
The  total  path  length  through  the  region  is  L.   At  the 
near  boundary,  where  s  =  L,  the  specific  intensity  is 
I(v).   The  change  in  specific  intensity  along  a  differ- 
ential path  length  ds  in  a  homogeneous  medium  is 

dl(v)  =  -  kl(v)  I(v,s)  ds  -kt(v)  I(v,s)  ds  +  j,_(v)  ds  +  j'T(v)  ds 


+  Jm(v)  ds. 


II-l 


The  nonthermal  emission  coefficient  is  j»>(v).   In  the 
interstellar  medium  the  nonthermal  absorption  processes 
are  not  important.   The  thermal  process,  which  will 
generate  or  absorb  radiation,  is  the  free-free  process; 
it  will  be  discussed  later.   The  coefficients  of  ther- 
mal continuum  emission  and  absorption  are  jj(v,s)  andKT(v,s) 
respectively.   The  applicable  coefficients  for  spectral 
lines  are  denoted  by  the  subscript  L. 

Let  j(v)  =  j'T(v)  +  JL(v)  and  k  (v)  =  kj{v)   +  KL(v),then  at  any 
point, 

I(v)  =  IT(v)  +  IN(v)  +  IL(v)  =  Ic(v)  +  IL(v).  Equation  II-l  then 
becomes 


dl(v)    .   .   J(v)    JN 


1 1-2 
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The  optical  depth,  x,  is  defined  by 


x(v)  =  I  k(v)  ds  =  tc(v)  +  rL(v, 


The  thermal  source  function,  S(v),  is  defined  by 


S(v) 


j(v)   JL(v)  +  JT(v) 

<(v)    <L(v)  +  Kj(v) 


II-3 


Since  there  is  no  nonthermal  absorption  coeffi- 
cient and  the  thermal  and  nonthermal  processes  are  in- 
dependent, there  is  no  nonthermal  source  function. 

A  general  solution  to  the  transfer  equation  can 
be  obtained  by  multiplying  Equation  1 1  -  2  by  the  inte- 
grating factor   exp(-x)  and  integrating  along  the  path 
length  to  get 


<JN(v)> 


I(v)  =  IQ(v)  e"T  +  <S(v)>  (l-e"T)  +  ^y^  (1-e"1). 


The  brackets,  <  >,  denote  the  average  over  the  path 
length.   The  intensity  of  the  nonthermal  emission  gen- 
erated within  the  line-emitting  region  is  1^  =  <JN>L 
and  <k>L  =  t ;  so,  the  general  solution  is 


I.. 
Kv)  -   IQe"T+  <S(v)>  (l-e"T)  +-f   (l-e"T) 


II-4 


Nonthermal  emission  decreases  in  intensity  with  increas 
ing  frequency,  in  general . 
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In  order  to  derive  a  solution  for  a  specific 
region,  the  source  function — which  depends  on  the  thermo- 
dynamic state  of  the  region — must  be  found.   It  is  im- 
portant, therefore,  to  consider  the  conditions  for  thermo- 
dynamic equilibrium.   This  is  the  subject  of  the  next 
section. 


Thermodynamic  Equilibrium 

To  be  in  thermodynamic  equilibrium  (TE),  a  gas 
must  have  the  following  properties. 

1.  The  temperature  and  pressure  throughout  the 
region  are  constant  and  are  equal  to  that  of 
its  surroundings;  hence,  there  will  be  no 
flow  of  energy  or  mass  through  it. 

2.  The  gas  will  radiate  as  a  black  body  at  that 
temperature . 

3.  Detailed  balancing  must  exist;  that  is,  for 
each  process  occurring  in  the  gas,  its  inverse 
must  also  occur  by  an  equal  amount  (e.g. 
ionization  and  recombination). 

4.  The  velocities  of  the  gas  particles  have  a 
Maxwellian  distribution. 

5.  The  population  of  any  excited  state  of  an 
atom,  ion,  or  molecule  is  determined  by  the 
Bol tzmann  equati  on . 
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6.   The  relative  population  of  any  level  of  ion- 
ization of  a  species  to  that  of  the  neutral 
species  is  given  by  the  ionization  (Saha) 
equation. 

The  observed  spectral  line  intensity  depends  on 
how  well  these  conditions  are  met.   Radio  recombination 
line  observations  of  a  variety  of  regions  indicate  that 
the  interstellar  medium  is  generally  not  in  local  thermo- 
dynamic equilibrium  (LTE);  that  is,  the  conditions  of  TE 
are  not  met  even  on  a  small  (local)  scale  (Hjellming  and 
Davies,  1970;  Andrews,  Hjellming,  and  Churchwell,  1971; 
Gordon,  1971;  Brown  and  Balick,  1973).   If  LTE  does  not 
exist,  the  transfer  equation  parameters  must  be  modified 
to  account  for  this.   In  the  following  work,  a  super- 
script asterisk  denotes  an  LTE  quantity.   The  second  con' 
dition  of  LTE  dictates  that  S*(v)  =  B(v,T),  the  Planck 
function  for  a  temperature  T.   For  radio  frequencies, 
this  simplifies  to 


v,T 


2kT   2 

c 


1 1-5 


In  the  case  of  non-LTE,  a  correction  parameter  is  defined 
such  that 


S(v)  =  A(v)  S*(v)  =  A(v)  B(v,T).     II-6 
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In  any  steady-state  situation  (LTE  or  non-LTE), 
the  total  per-unit  volume  emission  and  absorption  must 
balance.   Only  in  LTE  do  emission  and  absorption  balance 
at  each  frequency  so  that,  by  Kirchhof f ' s  Law,  j(v)  = 
k(v)  B(v,T). 

The  Thermal  Continuum 


Thermal  emission  and  absorption  are  the  result  of 
an  encounter  between  a  free  electron  and  a  positively- 
charged  ion  after  which  the  electron  remains  free  (cf. 
Mitchell  and  Zemansky,  1934,  p.  154).   From  a  classical 
point  of  view,  an  electron  moves  in  the  field  of  an  ion 
(e.g.  H  ),  is  accelerated  and  absorbs  a  wave  or  is  decel- 
erated and  emits  a  wave.   As  shown  in  Figure  II-2,  the  tra- 
jectory of  an  electron  with  a  given  velocity  is  altered 
by  an  amount  determined  by  the  impact  parameter,  p.   For 
radio  frequencies,  p  is  large  enough  that  the  entire 
trajectory  can  be  approximated  by  a  straight  line. 

Oster  (1961)  has  derived  an  expression  for  the 
thermal  absorption  coefficient,  which  includes  a  quantum 
mechanical  correction  applicable  for  temperatures  in  ex- 
cess of  5.5  x  10   °K.   The  derivation  of  <(  v  )  was  ob- 
tained by  considering  emission  from  a  single  electron  with 
specified  velocity  and  impact  parameter.   Then  j,(  v  )  was 
found  for  an  ensemble  by  integrating  first  over  a  suitable 
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range  of  p  and  then  over  a  Maxwellian  velocity  distri 
bution.   Then,  <T(v)  =  jy(v)/B(v,T)  and,  for  hydrogen 


K  (v)  =  3.014  x  10"2  rK5  v~2  ln(4.955x!0"2  T"2"5  v"1  )  Ng  I   Z2  N. 


1 1-7 


The  numerical  constants  are  functions  of  fundamental  con- 
stants as  shown  in  the  Appendix,  Equation  A-l.   The  sum- 
mation is  over  all  ionic  species  present.   Z.  is  the 
effective  charge  number  (ZR+  -  1)  and  N.  is  the  number 
den i sty  for  the  i    species. 

Altenhoff  et  al ■  (1960)  have  derived  an  approx- 
imate expression  for  Kj(v),  which  is 


Ky  (v  ) 


.235  x  10"2  T"1'35  v"2'1  N   $  Z2  N . .  II 
e  e  i   i   i 


To  bring  this  approximation  into  agreement  with  Oster's 
more  exact  formulation,  it  is  multiplied  by  a  correction 
factor  a(T  ,v).   Using  it,  the  optical  depth  then  becomes 


tt(v)  =  I   (v)  ds  =  8.235  x  10"2  a(Te,v)  T"1-35  v"2J  E] 


II-9 


ET  is  the  continuum  emission  measure  defined  by 
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(N  E  Z.  N.)  ds. 

e  .   i   "i 
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The  emission  measure  is  expressed  in  units  of  pc  cm 

The  factor  a(T  ,v)  has  been  tabulated  by  Mezger 

5 
and  Henderson  (1967)  for  the  ranges  100  ±   v  <    10   MHz 

and  3  x  103  <   T  <  15  x  103°K.   For  these  ranges,  a(T  ,v) 

—       g   —  3  g 

varies  between  0.7680  and  0.9971. 

The  emission  measure  is  often  given  by 

X 


T 


N2  ds 
e 
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The  effect  of  heavy  elements  on  the  size  of  E  T 

depends  on  the  fractional  ionization  of  hydrogen;  x  = 

N  /  N ,. .   The  atomic  species  with  an  ionization  potential 
e   H  r 

less  than  that  of  hydrogen  (13.6  eV)  may  be  totally 
ionized.   If  x  ^  10   ,  the  contribution  of  these  elements 
to  the  electron  pool  will  be  an  appreciable  fraction  of 
the  total  unless  their  abundances  are  less  than  those 
adopted  by  Spitzer  (1968,  p.  122).   The  most  abundant 

contributor  to  the  electron  pool,  other  than  hydrogen, 

-4 
is  carbon  which  has  an  abundance  N  ~  -    4  x  10   N ,, . 

The  continuum  emission  spectrum  from  a  region  of 

the  interstellar  medium  can  provide  some  information 

about  the  conditions  there.   Consider  Equation  1 1  -  3  for 

the  continuum;  that  is,  I.  =  0.   Let  I   =  I..  =  0  and, 

L  ON 
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ssuming  LTE,  S(v)  =  B  (v,T).   Then,  using  Equation  II-5, 


if    \         2kT   2  n  Q"TT 
I(v)  =  — —  v   (1-e 

c 
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When  the  frequency  is  high  enough  that  the  gas  is  opti- 
cal ly  thin,  tt  <<  1 ,  and 

.,  v    2kT   2 
I(v)  =  — g-  v   tt- 

c 

Since,  by  Equation  1 1  -8 ,  tt  a  v"  '  ,  I(v)  *   v"  '   and  the 
spectrum  is  nearly  flat.   For  low  frequencies  such  that 
the  gas  is  optically  thick,  tt  X    1 • 5  and  !  ( v  )  «  v  .   The 
turnover  frequency,  \>j,    is  that  for  which  tt  =  1,  al- 
though sometimes  the  value  of  1.5  is  used.   Setting  tt 
equal  to  one  in  Equation  1 1 -8  and  solving  for  vT  gives 


vT  .  1.893  x  102  a(T  .v)-2'1  El2''  T"0'643 
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If  in  some  region  for  which  vT  can  be  measured,  either 
ET  or  T  is  already  known,  the  other  parameter  can  be 
found.   One  independent  source  of  information  is  the  in 
tegrated  radio  recombination  line  power.   Radio 
recombination  lines  will  be  considered  in  the  next 
section. 


Radio  Recombination  Lines 

Radio  recombination  lines  arise  from  transitions 
between  levels  with  large  principal  quantum  numbers; 
n  ?■  40.   They  result  from  the  capture  (recombination)  of 
an  electron  by  an  ion  into  an  upper  level.   The  nomen- 
clature for  such  a  transition  is  Xn(An)  where  X  designates 
the  recombining   species;  n  is  the  principal  quantum  num- 
ber for  the  transition  which  is  specified  by  a  for  An  =  1  , 
B  for  An  =  2,  etc.   The  predominant  element  from  which 
recombination  lines  are  observed  is  hydrogen;  however, 
radio  recombination  lines  of  helium,  carbon,  and  possibly 
heavier  elements  have  been  observed  (Gordon,  1974; 
Chaisson,  1973;  Pedlar  and  Hart,  1974). 

The  emission  frequency,  v   ,  for  the  transition 
between  the  m   and  n   levels  (m  >  n)  is  given  by  the 
Rydberg  formula  (cf.  Shore  and  Menzel  ,  1968,  p.  13) 


v    =  c  R   Z   (-0-— o) 
mn      xx  v  2  2' 
n   m 


2c  R   Z2  ^ 
x   x  n3 
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R   is  the  Rydberg  constant  for  the  species  X  and  Z   is 
x  x 

its  effective  charge  number.   Radio  recombination  line 
frequencies  for  hydrogen  and  helium  have  been  tabulated 
by  Lil ley  and  Palmer  (1968)  . 

The  recombination  line  strength  depends  on  the 
number  of  atoms  in  the  region  which  are  excited  to  the 
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upper  level  of  the  transition.   Since  the  gas  may  not  be 
in  LTE,  it  is  necessary  to  calculate  the  LTE  populations 
and  the  departure  parameters,  which  indicate  the  degree  to 
which  the  gas  is  out  of  LTE.   These  matters  will  be  dis- 
cussed in  the  following  subsections. 

Thermodynamic  Temperatures  and  Level  Populations 

Radio  recombination  line  emission  regions  contain 

excited  neutral  atoms,  ions,  and  electrons.   The  number 

_  3 
densities  (in  units  of  cm   )  for  hydrogen  are  N„0,  NH+, 

and  N   respectively.   N„  is  the  total  hydrogen  density. 
In  LTE  there  is  an  equilibrium  situation  where  recombin- 
ations are  balanced  by  ionizations.   What  is  desired  is 
a  statistical  relation  between  the  kinetic  temperature  of 

the  qas,  N..+  ,  N  ,  and  N  ,  the  number  density  of  atoms  in 
hr   e       n 

the  n    level  of  excitation. 

For  a  gas  in  LTE,  the  population  number  density  of 
ions  in  the  n   excited  level  and  the  i    ionization  state 
is  given  by  the  Boltzmann  equation  (cf.  Aller,  1963,  p. 
112): 


N  .  =  N  .  —  exp(-  x  •  /kT   )  . 
m    oi  g    Y  v  A l n   ex 
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This  equation  is  based  on  the  condition  that  the 
assembly  has  a  Maxwellian  velocity  distribution.   In 
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Equation  11-15,  N  .  is  the  LTE  population  of  the  ground 

state;  g   and  g   are  the  statistical  weights  of  the  ex- 
an     3o 

cited  and  ground  states  respectively;  x-„  "is  the  excitation 
energy  within  the  i    state  of  ionization;  and  Tgx  is  the 
excitation  temperature.   The  statistical  weights  are  the 
number  of  sublevels  in  the  n    level;  g  =  ( 2 J  +  1).   J 
represents  the  total  angular  momentum  —  spin  plus  orbital  — 
of  the  level.   Under  the  conditions  of  LTE,  the  excita- 
tion temperature,  T   ,  is  equal  to  the  kinetic  temperature 
r  ex 

of  the  gas. 

The  radio  recombination  lines, which  have  been 
observed  involve  highly  excited  neutral  atoms,  for  ex- 
ample H  I,  He  I,  and  C  I.   So, the  desired  population  is 

N   expressed  in  terms  of  the  electron  density  and  the 

no    v 

respective  number  of  the  singly  ionized  species.   For 
the  spec  i  es  X , 


N>)  -  4.143  x  lO"16  N*(X)  Ne  T^5  ^  exp(-ein(X)/k  Tex) 
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The  details  of  the  derivation  of  this  equation  are  given 
in  the  Appendix;  it  results  from  placing  the  values  of 
the  fundamental  constants  into  Equation  A-5.   The  par- 
tition function,  u,(T),  is  given  by  Equation  A-2.   ein(x) 
is  the  threshold  energy  which  is  (Equation  A-6) 
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X    RAY 


Figure  11-2.      The  Free-Free  Process.   The  elec 
tron,  e,  approaches  the  central 
ion.   The  impact  parameter,  p , 
is  shown  for  two  cases. 
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eln(x)=-(Xl-X1n)   -  157900     (Z^/n2)   k.  11-17 


The  kinetic  energy  of  a  particle  with  mass  M  is 


eK  =  1/2M  v"  =  3/2  k  TK. 


If  the  gas  has  more  than  one  constituent,  each  with  a 
different  mass,  the  kinetic  temperature  of  each  constitu 
ent  will  be  the  same  although,  because  of  the  different 
masses,  the  RMS  velocity  of  each  constituent  will  be 
different.   The  RMS  velocity  of  the  component  with  mass 
M  i  s 


2    1/2         ,3k    TKJ/2  TT    , 

<V     >  =     ( jjj — )  .  H-1 


If    the    gas    particle    velocity    distribution    should 
suddenly    become    non-Maxwel 1 ian ,    collisions    involving    a 
deflection    of    90°    will    thermalize    them    to    an    equilibrium 
temperature    unless    the    system    is    thermodynamical ly    unstable 
(Spitzer,    1968,    p.    90).       In    a    situation    where    there    is    a 
nonequi  partition    of    energy,    the    electrons    will    have    the 
largest    kinetic    temperature    because    of    their    low   mass. 
The    energies    of    the    electrons    and    ions    will    equalize    first 
because   of    the    long-range    interactions    between    them. 
Following    this,    there    will    be    elastic  collisions     between  the 
ions  and  neutral   atoms.     Collisions  between  electrons  and  atoms  are  not 
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as  an  effective  means  of  equalizing  the  kinetic  energies 
as  are  collisions  between  ions  and  the  electrons  (Spitzer, 
1968,  p.  140). 

Population  and  Depopulation  Mechanisms  of  High  n-Levels 

The  high  n-levels  are  populated  and  depopulated  by 
several  mechanisms  which  are  discussed  by  Brockel hurst 
(1970),  Dupree  and  Goldberg  (1970)  and  M^Daniel  (1964, 
p.  588).   The  parameters,  such  as  cross-sections,  proba- 
bilities, and  coefficients  are  discussed  in  general  by 
these  authors  as  well  as  by  Spitzer  (1968),  Kaplan  and 
Pikelner  (1970),  and  by  the  workers  listed  in  Table  II-l 
who  have  calculated  high  n-level  populations.   The  mecha- 
nisms involved  are  as  follows. 


1.  Electron  recombination  followed  by  a  cascade 

2.  Dielectric  recombination  (in  atoms  with  more 
than  one  electron). 

3.  Excitation  and  de-excitation  induced  by 
el ectron  collisions. 

4.  Collisional  ionization  followed  by  three- 
body  recombination  leaving  an  excited  atom 
X':  H  +  e  ■>  X+  +  e  +  e  and  X   +  e  +  e  ■+ 

X'  +  e. 

5.  Stimulated  and  spontaneous  emission. 
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6.   Fine  structure  excitation;  i.e.  redistribution 
of  angular  momentum  with  a  level:   n  £  -*■  n  &  '  . 

Excitation  to  high  n-levels  from  the  ground  level  by  radi- 
ative excitation  is  not  important  in  H  II  regions  (Dyson, 
1969;  Dupree  and  Goldberg,  1970). 

To  account  for  the  effects  of  a  departure  from  LTE 
with  respect  to  level  populations,  the  departure  coeffi- 
cient, b  ,  is  employed; 


b   =  N  /N 
n    n   n 
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The  departure  coefficient  is  a  function  of  the  electron 
temperature  and  densi ty  and  the  intensity  of  the  radiation 
field.   When  b  is  calculated  for  a  type  of  region,  the 
radiation  temperature  is  chosen  on  the  basis  of  specific 
environmental   conditions  of  the  region  and  bn  is  tabu- 
lated or  plotted  as  a  function  of  n  for  specific  values  of 

T   and  N  .   For  example,  Dupree  (1972a)  has  done  this  for 
e      e 

departure  coefficients  in  a  cool  region.   The  coefficients 

were  first  calculated  for  an  isolated  region  where  IQ(v) 

4 
=  0  and  then  the  effects  of  an  adjacent  10   °K  H  II  region 

on  the  b   values  were  demonstrated, 
n 

Statistical  equilibrium  requires  that  there  are 
as  many  transitions  into  a  level  as  there  are  transitions 
out  of  it.   The  equation  of  statistical  equilibrium  for 
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the  nth  level  is  of  the  form  (Seaton,  1964;  Dupree,  1972a, 

b) 


n-1 


n-1  n+p' 


t"fo 


?      Vj   +    T   .Cn,j  +  ^+1Cn,1+    £-.,    Ln,m       r    ,       n,n,, 
j=n         'J     J=n-p         °    j=n+l  m=n+l  j-  I 


I(v 


=  y     (m         C  +  N         C  )   +   (  y       N     B 

*-    [  n-p     n-p,n       'n+p     n+pV       vmf^+1     m     m,n 

p=l 


+   L      N-   B     n)    I(v)+    £       N     A         +NeNi    (a^an  +  C.jn) 

j=n       J     J'  j=n+l 
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The  symbols  used  in  Equation  11-20,  which  have  not  yet 
been  defined,  are  as  follows: 


A    =  Einstein  spontaneous  emission  probability 

m ,  n 


n  ,m 


m ,  n 


m  >  n  ; 
=  Einstein  absorption  probability,  m  >  n; 
=  Einstein  stimulated  emission  probability, 

m  >  n ; 


C.  ,  =  collision  probability  for  a  transition 
J  >k 

between  the  j    and  k    levels; 

C   .  =  collisional  ionization  probablity  from 
n  ,  i 

1 evel  n ; 
C.    =  three-body  recombination  probability; 


and 
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a,r  =  rate  of  radiative  recombination  into  level  n; 

n 

a      =  rate  of  di  electronic  recombination  into  level 
an 

n  (  =  0  for  hydrogen )  ; 
m '  =  the  maximum  level  considered  in  cascades  from 

hi  gher  1 evel s  ; 
p'  =  the  maximum  p  for  bound  collisional  transitions, 

p  =  | An | ; 

n   =  the  lowest  level  of  the  atom  considered  for 
o 


the  computation; 


N.  =  the  ion  number  density 


Collisional  Transition  Probabilities  and  Cross-sections 


In  order  to  compute  the  departure  coefficients, 
means  of  simplifying  Equation  11-20,  limiting  the  number  of 
equations  or  limiting  the  range  of  the  summations  should 
be  considered.   At  some  level  n  =  m ' ,  the  electron  is  close 
enough  to  the  continuum  that  the  collisional  ionization 
probability,  C   .,  is  much  greater  than  the  probability  for 


a  spontaneous  transition  to  lower  levels,  A 


m ,  n 


The  levels 


beyond  the  m '  cut-off  level  are  in  thermodynamic  equilibrium 

and,  b   =  1  for  n  >m'. 
n 

Sejnowski  and  Hjellming  (1969)  find  that  inclusion 
of  collisional  transitions  of  p  =  | A  n  |  >  1  will  lower  the 
value  of  b   by  only  a  few  percent  from  that  considering 
only  transitions  of  p  =  1 .   Consequently,  many  workers 
limit  their  calculations  to  p '  =  1  . 
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Collisions  involving  a  change  of  I-    state  within 

a  level,  n,£  ■+    n ,  I   +  1,  occur  frequently.   In  cold  hydrogen 

qas  with  N  >  10   cm   ,  fine  structure  collisions  by  pro- 
g- 
tons  as  well  as  el ectrons,  occur  more  frequently  within 

levels  for  which  n  >  100  than  do  radiative  processes 

2 
(Pengelly  and  Seaton,  1964).   Because  of  this,  all  n 

£-states  are  equally  populated.   Also,  at  low  temperatures, 
fine  structure  excitation  by  electron  collisions  provide 
the  main  cooling  mechanism  of  the  gas  (Dalgarno  and  McCray, 
1972).   By  this  process,  electron  kinetic  energy  is  con- 
verted to  low  frequency  radiation  which  leaves  the  gas, 
thus  removing  energy  from  the  thermal  pool. 

The  total  population  of  level  n  is  divided  into 
Jl-states  according  to 


n,£ 


(2&  +  1  ) 


(Dupree,  1969).   The  total  population  is  N   =  E  N   . 
For  hydrogen,  the  levels  are  degenerate  and  b 


n,A 


b  ; 


b   is  independent  of  I. 
n        r 


The  collisional  transition  rates,  C.  ,  for 


=  j  +  p,  depend  on  cross-sections  Q.  , 


by 


Cj,k  ■  Ne  -  «j.k  (v)> 
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The  average  is  obtained  by  integrating  over  a  Maxwellian 
velocity  distribution. 

The  bound-bound  cross -secti ons  have  been  computed 
by  means  of  one  of  three  approximations  or  methods: 
binary  encounter  method,  impact  parameter  method,  and  a 
method  based  on  the  correspondence  principle.   The  binary 
encounter  method  uses  a  classical  approximation,  which  is 
valid  for  electron-atom  collisions  if  the  time  of  inter- 
action between  the  incident  and  bound  electrons  is  short 
compared  to  the  orbital  period  of  the  bound  electron  and 
the  collision  takes  place  in  a  region  smaller  than  the  size 
of  the  atom.   This  time  is  sufficiently  short  when  the 
energy  transferred  between  the  two  electrons  is  large. 
With  this  condition,  the  energy  transferred  is  assumed  to 
be  equivalent  to  that  transferred  between  two  free  elec- 
trons.  Because  of  this  requirement,  the  binary  encounter 
method  gives  best  results  for  ionizing  collisions;  the 
results  are  poor  for  p  =  1.   In  calculating  cross-sections, 
Brockelhurst  (1970)  uses  the  binary  encounter  method  for 
p  >  6.   The  application  of  the  theory  involved  is  discussed 
at  length  by  Burgess  and  Percival  (1968). 

The  impact  parameter  method  for  electron-atom 
collisions  uses  the  assumption  that  the  incident  electron 
follows  a  straight-line  trajectory  past  the  target  atom 
with  an  impact  parameter  p  (cf.  Figure  1 1  -2 ) .   Seaton 
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(1962)  computes  the  collision  cross-section  for  a  bound 
transition  (n,£  •+  n',  £')  for  both  weak  and  strong  cou- 
pling.  For  weak  coupling, 


'j.k 


Pk-(p)  2ttp  dp, 


where  j  =  n,£,  k  =  n',  I ' ,  and  p   is  the  cut-off  value  of 
p>which  is  approximately  equal  to  the  smaller  of  the  two 
atomic  radii  of  the  j  and  k  states.   P^Cp)  is  the  transi- 

K  J 

tion  probability  as  a  function  of  impact  parameter. 

For  strong  coupling,  Seaton  (1962)  defines  p, 

such  that  P..(p-.)  =  1/2.   For  p<p,  ,  Pj.  oscillates  about  a 

2 
mean  value  of  1/2  so  that  in  this  range,  Q.  .  =  1/2-rrp.  . 

J  ,  K  1 

For  a  transition  of  p  =  1,  most  of  Q.  .  comes  from  p>p,  , 
and 


/ 


Qj,k  =  1/2^pl   +  J         P.  .(,  )  2,r    dp 
P 1    "■  J 


Seaton  (1962)  calculates  the  cross-sections  for 
the  weak-  and  strong-coupled  cases  and  accepts  the  smaller 
of  the  two  values.   Pu-(p)  is  calculated  with  the  sym- 
metrized,  time-dependent  perturbation  theory.   Brockelhurst 
(1970)  claims  that  the  impact  parameter  method  will  not 
work  for  p  >  1.   He  assumes  the  case  of  strong  coupling  and 
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finds  that  most  of  Q  comes  from  p  <  p-i  for  which  the  per- 
turbation theory  is  not  valid. 

The  classical  method  only  works  for  large  p;  the 
semi -cl assical  impact  parameter  method  is  valid  only  for 
p  a  1.   The  correspondence  principle  method  can  be  used  to 
obtain  cross-sections  for  intermediate  transitions.   All 
variations  of  the  correspondence  principle  state  that  for 
large  scale  systems, where  the  uncertainty  is  much  larger 
than  Planck's  constant,  h,  the  result  of  a  quantum  mechani- 
cal approach  is  identical  to  that  using  a  classical  ap- 
proach.  Percival  and  Richards  (1970a)  show  that  the  Bohr 
correspondence  principle  can  be  applied  to  the  energy 
transfer  of  a  collisional  excitation  for  the  case  of  a 
weak  interaction  and  if  the  quantum  and  classical  first- 
order  perturbation  theories  are  both  valid.   Also,  the 
collisional  probability  must  be  much  larger  than  that  of 
a  radiative  transition,  a  condition  which  is  satisfied  for 
n  ~  100.   With  the  the  Bohr  principle  for  collisions,  the 
mean  net  energy  transferred  during  the  collision  is  equiva- 
lent to  the  energy  transferred  from  independent  Fourier 
components  of  order  +p  and  -p  of  the  classical  atomic 
orbital  motion.   P . . ( p)  is  computed  from  these  Fourier 
components . 

For  the  case  of  a  strong  interaction,  Percival  and 
Richards  (1970b)  show  that  the  Heisenberg  correspondence 
principle  relates  the  quantum  mechnical  matrix  elements 
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with  the  Fourier  components  of  classical  atomic  orbital 
motion.   The  results  for  H  II  regions  are  published  in  a 
separate  paper  (Percival  and  Richards,  1970c). 

For  collisional  ionization,  either  the  binary  en- 
counter method  (Burgess  and  Percival,  1968)  or  the  classi 
cal  Thompson  theory  (Dyson,  1969;  Seaton,  1964)  can  be 
used  to  obtain  the  collision  rates.   From  the  Thompson 
theory 


C   .  =  3.46  x  10"15  n2  N   T~1/2 
n  ,  i  e 


The  inverse  of  collisional  ionization  is  three-body  re- 
combination.  Using  the  equation  of  detailed  balancing 


N   C   .  =  N.  N   C.  n 
n   n ,  i     l   e   i  ,  n 


Radiative  and  Dielectronic  Recombination  Coefficients 


The  radiative  recombination  coefficient  for  cap- 
ture into  level  n  is  represented  by  a  .  The  partial  re- 
combination coefficient,  a,.v,  "> s  that  for  captures  into 
all  levels  above  and  including  the  k   ;  n  >  k.   It  is 

given  by  a,.\    =   Z°°a  .  .   Usually  k  =  2.   The  total  re- 
IKj    j=k  J 

combination  coefficient  ,  a,,\,  includes  radiative  cap- 
tures into  the  ground  level.  Seaton  (1959)  has  derived 
formulations  for  a,,x  and  a,^  using  a  symptotic  expansion 
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of  the  Kramers-Gaunt  factor.   With  a  notation  similar  to 
that  of  Spitzer  (1968,  p.  117)  , 


a(k)    =    2.065    x    10"11    Z2    T"1/2    Ck    (TR 


11-22 


The  function  £k(T)  involves  the  asymptotic  expansion. 
Spitzer  (1968,  p.  117)  has  tabulated  £,  and  C2-   Tney  are 
pi otted  as  Figure  1 1 -3 . 

Dielectronic  recombination  occurs  in  atoms,  which 
have  two  excitable  electrons.   This  process  is  very  im- 
portant in  establishing  the  high  level  populations  of 
carbon  atoms  and  is,  in  fact,  responsible  for  these  levels 
being  over  populated  with  respect  to  the  LTE  populations. 
It  takes  place  in  two  steps,  the  first  of  which  is  a 
capture  of  an  electron  into  a  high  level,  n.   The  energy, 
rather  than  being  radiated,  is  transferred  to  a  bound 
electron, which  rises  to  an  excited  level;  n'  <<  n.   This 
is  an  unstable  condition  and  the  second  step  is  to  dis- 
pose of  the  energy  of  the  lower  excited  level  by  means 
of  a  downward  transition  to  a  stable  level.   The  energy 
difference  is  either  radiated  —  leaving  the  atom  with  a 
single,  highly  excited  level  —  or  it  is  transferred  to  the 
highly  excited  electron  which  leaves  the  atom  by  auto- 
ionization,  a  process  called  the  Auger  effect  (McDaniel, 
1964,  p.  590;  Shore  and  Menzel  ,  1968,  p.  93).   The  radi- 
ative process  dominates  (Dupree,  1969).   The  dielectronic 
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recombination  mechanism  is  more  effective  at  populating 

highly  excited  states  than  is  the  radiative  recombination 

d     r  4 

mechanism;  i.e.,  a_  >  a  ,  at  least  for  T  ^  10   °K  and 
n     n  e 

N   >  1  cm"3  (Dupree,  1969). 
e  v   r 


Departure  Parameters,  b   and  6 
— t 2 — n mn 


It  will  be  shown  later  that  the  spectral  line 

strengths  depend  not  only  on  the  departure  coefficient, 

b  ,  but  —  since  the  transition  involves  two  b  's  —  it  also 
n  n 

depends  on  the  derivative  of  b   with  respect  to  n.   For 
r  n 

that  reason,  when  b   is  computed,  a  parameter  related 
to  the  derivative  is  also  computed; 


11-23 


where  m  >  n. 

Departure  parameters  have  been  published  by  sev- 
eral workers.   Information  about  this  work  is  summarized 
in  Table  II-l.   The  first  column  gives  the  reference  to 
the  work.   The  parameters  calculated  by  Dupree  (1969) 
and  Hoang-Binh  and  Wal msley  ( 1 974  )  are  for  carbon.   The 
second  and  third  columns  give  the  range  of  n  over  which 
the  parameters  are  tabulated  or  plotted.   Columns  4,  5 
and  6  give  the  range  parameters  used  in  the  calculations: 
n  ,  m',  and  p'  (cf.  Equation  11-20).   The  seventh  column 
indicates  the  method  by  which  the  collision  cross-sections 
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were  calculated.   The  reference  to  the  source  of  the  cross- 
sections  is  given  in  the  parentheses  following  the  designa- 
tion.  Since  the  departure  parameters  are  functions  of  the 
electron  temperature  and  density,  the  last  column  of  the 
table  gives  the  range  of  these  quantities  for  which  the 
parameters  were  computed. 

All,  except  Hoang-Binh  and  Walmsley  (1974),  do  not 
consider  incident  radiation  in  the  calculations.   Dupree 
(1972)  does,  however,  show  the  effect  of  an  adjacent  H  II 
region  on  a  line-emitting  region  with  a  temperature  of 
20  °K.   The  radiation  field  from  the  H  II  region  increases 

b   si i  ghtly  for  n  >  80. 
n     ?        j 

The  physical  significance  of  departure  from  LTE 

can  be  illustrated  by  showing  the  relation  which  n,  T  , 

and  N   have  on  b   and  6   .   This  is  done  with  Figure  1 1  - 4 
e  n      mn 

4  3 

in  which  are  drawn  curves  for  T   =  10   and  7.5  x  10   °K 

e 

computed  by  Sejnowski  and  Hjellming  (1969)  and  curves  for 

T   =  100  °K  computed  by  Dupree  (1972a).   For  each  temper- 

3-3 
ature,  there  are  separate  curves  for  N   =  10   and  10  cm 
r  e 

As  n  decreases,  collisions  become  less  important 
in  determining  level  populations.   For  small  n,  radiative 

processes  dominate  and  the  curves  show  an  asymptotic  be- 

4  3 

havior.   The  heavy  portion  of  the  10   and  7.5  x  10   °K 

curves  in  Figure  1 1  -  4  is  the  radiative  asymptote.   De- 
population by  stimulated  emission  is  much  less  frequent 
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260 


Figure  II-4.   Departure  Parameters  for  Hydrogen  as 

a  Function  of  n,  Te,  and  Ne.   The  curves 
with  a  solid  line  are  for  Te  =  104oK  and 
with  a  dashed  line  are  for  Te  =  7.5  x  1 0Jo K 
(Sejnowski  and  Hjellming,  1969).   The 
dotted  curve  is  for  Te  =  100°K  from  Du- 
pree  (1972a). 
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than  by  spontaneous  emission  (Gordon,  1974),  except  when  an 

intense  radiation  field  is  present. 

As  N   decreases,  collisions  become  less  frequent 
e 

and  the  radiative  asymptote  extends  to  larger  n.   Then,  for 
hydrogen,  there  is  a  decrease  in  the  rate  of  collisional 
excitation  of  lower  n-levels  which  causes  a  decrease  in 
level  population  with  respect  to  the  LTE  values  as  n  de- 
creases.  In  the  case  of  carbon,  however,  dielectronic  re- 
combination is  efficient  enough  in  populating  upper  levels 
that  the  population  will  exceed  the  LTE  values  for  n  >  20 
and  N   <  10   cm"  .   In  any  case,  the  populations  approach 
their  LTE  values  as  n  ■+   °°  (Dupree,  1969). 

It  should  be  noted  that  for  hydrogen,  b   and  6mn 
are  both  positive  in  the  radio  regime.   For  any  two  levels 
m  and  n,  where  m  >  n,  b   is  closer  to  one  than  is  b  .   In 
that  sense,  the  upper  level  is  "overpopulated"  with  respect 
to  the  lower  level.   If  the  slope  of  the  bn  curve  is  large 
enough,  there  is  an  enhancement  of  the  line  intensity  over 
that  obtained  from  a  region  in  LTE  even  though  the  level 
population  is  less  than  its  LTE  value.   This  occurs  be- 
cause of  stimulated  emission  induced  by  an  exterior  source 
or  the  region's  own  thermal  radiation  field.   An  examin- 
ation of  the  radio  recombination  line  absorption  coef- 
ficient and  non-LTE  solutions  of  the  transfer  equation 
will  show  how  this  happens. 


39 


Radio  Recombination  Line  Absorption  Coefficient 

The  rate,  R(v)  at  which  a  bound-bound  absorption 

takes  place  between  enrgy  levels  n  and  m  (m  >  n)  is  equal 

to  the  product  of  the  number  of  absorbers,  N  ,  the  Einstein 

absorption  probability,  B   ,  and  the  energy  flux,  I(v) 
r  nm 

dn/4iT.   The  degree  of  absorption  is  modulated  by  the  spec- 
tral line  profile  function,  tj> .  ( v ) .   Then 


dQ. 


R(v)  =  Nn  ^(v)  Bnm  I(v)  ^ 


The  energy  per  photon  is  hv   .   The  absorption  coefficient, 

i 

k,  (v),  is  defined  such  that  the  energy  absorbed  along  an 
increment  of  path  length  ds  is 


f  HO 

deA=-.L(v)  I(v)  dfl  -  -Nn  *L(v)  BnmI(v)  hvnm^  ds. 

However,  energy  may  be  emitted  by  stimulated  emission  of 
the  amount 


dp, 


de,  -  jj(v)  du  =  NM  *L(v)  Bmn  I(v)  hvmn  ^ds 


and  by  spontaneous  emission  by  the  amount 


dss  =  js(v)  d^  -  Nm  *L(v)  Amn  hvmn  ^  ds 
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\h       (v)  is  the  emission  line  profile,  which  is  usually 
r  mn 

identical  to  the  absorption  profile.   Inserting  deA,  de j 
and  de ^  in  the  transfer  equation  (Equation  1 1  -  2 )  gives 


Let  <$>.   (  v)  =  <jj,  (v)  and  define  a  total  absorption  coef- 
ficient, k  .  (v  ) ,  by 

k.  (v)    =    (N      B         -    N  B      )    <J>.  (v)    5j 

L  n      nm  m  mn         L  4tt 


B         N 
hv    , /    \    R         n      mn      m 
-r-  cf>  (v  )    B         U  - 


4tt 


nm    V      Bnm    Nn 


Usinq    q      B        =    q      B         (Mihalas,    1970,    p.    83),    the    Boltzmann 
y    yn      nm        3m      mn 

equation,    and    the   departure    coefficients,    the    absorption 
coefficient    becomes 


M=%~    *L<v)  Brm  bn  Nn  [1   --b-exP(^/kTex^ 
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For    LTE,    b      =    b      =    1    and    the    LTE    absorption    coefficient    is 


hv 


kl   (v)   =  f    <frL(v)   Bpm  Nn   [1    -  exp(-hv/k  Tex)]. 


11-25 
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The  Einstein  probability  can  be  written  in  terms 
of  the  oscillator  strength,  f 


nm 


nm 


a      2   2 

4  tt   e 

hv  Mc 


nm 
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In  the  above  equation,  e  i s  the  el ec tron  '  s  charge;  M  is 
its  mass.   Oscillator  strengths  for  hydrogen  radio  re- 
combination lines  have  been  computed  by  Goldwire  (1968) 
and  Menzel  (1969)  . 

Substituting  Equation  11-26  into  Equation  11-25 
and  assuming  that  hv  <<  kT   ,  the  absorption  equation  be 
comes 


*      h  v   it  e        * 
ex 


m 


N   can  be  obtained  from  the  Saha-Bol tzmann  equation 
n 

(Equation  A-5)  and  v   can  be  obtained  from  an  approxi- 

2     3 
mat  ion  to  the  Rydberg  equation;  vmn  =  2c  RH  Z  An/n  . 

Using  these  and  the  appropriate  constants  (R.,  =  1.097 
x  105,  Z  =  1,  u(T)  =  1,  and  gp  =  2n2),  the  LTE  absorp- 
tion coefficient  becomes 


k*(v)   =   1-071    x  107  *(v)  ffm  N.   Ne  r2-5  exp(en/k  Te 


/k  T  )     11-27 

"X 


where  e  /k  =  157890/n  .   The  LTE  optical  depth  is 
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t*(v)   =  1.071    x  107<D(v)^fnmr2-5  exp(En/kTex)   E,_.  11-28 


E.     is    the    line    emission    measure   defined    by 


N.N      ds 

^      e 


and  N.  is  the  number  density  of  the  recombining  ion. 

Following  the  notation  of  Dupree  and  Goldberg 
(1970) ,  0   is  defined  by 


lm  =   [1  -  J*  exp(-hv/k  Tex)]  /  [1  -  exp(-hv/k  Tex)]. 
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Then  kl(v)  =  bn  6mn  ^(v)  and  T|_(v)  =  \    3mn  tl(v). 

In  most  cases  for  radio  recombination  lines,  hv/ 

k  T   is  small  enough  that  the  exponentials  involving  them 
ex 

are  equal  to  1  -  hv/k  T   .  Then, 


k  T 


ex 


mn         hv 


1   '  b       lkT       ; 
n  ex 


k  T  b     -  b    -i 

i  ex     /  m         n% 

1   "     hv         [     b         ' 


Usinq  6   ,  which  was  defined  in  Equation  11-23, 
s   mn 


1         k  Tex 
3mn  =  {T^)    (1  "  "nT^  6mn 
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If  Ab  =  b   -  b   is  small,  Ab  ■>  db  and 
m    n 
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mn 


d_b  _  d(1n  b) 
b„,     dn 


A  n  , 


m  =  n 


which  is  a  notation  often  used. 

Gordon  (1971)  and  Hjellming  and  Gordon  (1971) 
define  a  different  3  than  that  defined  in  Equation  11-29; 


k  T 


mn 


(1 


ex 


)  -  1 


k  T 

Ht\T 


ex   d  1 n  b 
dn 


An 


This  is  given  for  reference  only  since  the  definition  of 

B    in  Equation  11-29  will  be  used  in  the  present  work. 
mn 

Li  ne  Prof i 1 es 


There  are  three  line  profile  f uncti ons, whi ch  are 
relevant  to  radio  recombination  line  analysis.   They  are 
the  Gaussian,  Lorentzian,  and  the  convolution  of  these 
two,  the  Voigt  function.   The  Gaussian  function  arises 
from  the  Doppler  shift  resulting  from  thermal,  turbulent, 
expansion,  and  contraction  motions  of  the  line  emitter 
relative  to  the  observer.   A  Lorentzian  profile  for  a 
radio  recombination  line  results  from  pressure  broadening. 
Natural  broadening,  or  radiation  damping,  also  produces 
a  Lorentzian  profile,  but  this  is  such  a  small  effect  that 
it  is  impossible  to  detect  in  the  observed  line  profiles. 
If  pressure  broadening  occurs,  the  effect  would  be  blended 
with  the  Doppler  broadening  and  a  Voigt  profile  would  result 
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Doppler  broadening.   Mainly,  thermal  broadening 
will  be  considered.   The  thermal  motion  of  particles  in  a 
gas  with  a  kinetic  temperature  T.,  has  a  Maxwellian  velocity 
distribution;  so,  the  number  of  particles,  out  of  a  total 
of  N,  which  have  a  radial  velocity  between  v  and  v  +  dv 
i  s 


d  N(v 


0/ 


2   2 
exp[-(v-vc)  /v  J  dv 


In  this  equation,  v   is  the  most  probable  velocity  and 

v   is  the  velocity  of  the  ensemble  as  a  unit.   v   = 
c  o 

1  /  ? 
(2k  T„/M)    where  M  is  the  mass  of  the  particles.   The 
K 

intensity  of  the  spectral  line  at  a  frequency  v^  is  pro 
portional  to  the  number  of  emitters  Doppl er-shi f ted  to 
that  frequency.   The  line  profile  is  a  Gaussian; 


1  2   2 

(v)  = exp[-(v-vc)  /v  ]  dv 

o  ,r^~ 
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The  profile  function  is  a  probability  distribution 


<Mv)  dv  =  1  . 


The  power  spectrum  is  usually  displayed  in  terms 
of  antenna  temperature  as  a  function  of  velocity.   Velo- 
city is  physically  related  to  the  emitter  and  the  line 
half-power  width  is  then  invariant  to  the  line  rest 
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frequency.   The  peak  line  antenna  temperature,  T'(c), 
is  that  for  which  v  =  v  .   The  temperature  profile  is 


T'(v)  =  T'(c)  exp[-(v-vc)2/v^]. 
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The  integrated  line  power,  which  is  the  area  under 
the  profile,  is 


P  =   I   T1  (v)  dv  =  T'  (c)  vQ  JH 
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Since  v   is  not  known,  it  is  desired  to  substitute  it  with 
o 

a  measurable  quantity.   When  the  line  power  drops  to  one 
half  of  its  peak  intensity,  v  =  v,  and 


T'(Vl)  =  1/2  ^-V-  -  v-^p  exp[-(Vvc)2/v2]. 


O  V  IT      0  V  TT 


The  Doppler  half-width  is  defined  by  AvD  =  2(v,-v  );  then 

8k  T, 


AVD  =■  2vo  T^TT  -  {r-jf 


D   In  2)1/2 
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Tn  is  the  Doppler  temperature  which  is  the  effective 
temperature  of  the  gas  if  the  line  width  were  the  result 
of  thermal  motion  only. 

The  integrated  line  power  in  units  of  °K  km/s  is 


P(°K    km/s) 


<4-rir2»1/2   t,'c»«d 
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However,  P  is  usually  given  in  units  of  °K  kHz;  so, 


(°K  kHz)  =  (^-^)1/2  (*£)  T'(c)  Av, 


4  In  2 


1 .0645  [-£■)    T'  (c)  AvQ. 
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In  Equation  11-35,  v   is  the  rest  frequency  measured  in 
kHz;  the  velocity  of  light,  c,  and  the  Doppler  width  are 
both  measured  in  km/s.   Likewise,  the  integrated  optical 
depth ,  t ,  ,  is 


V 

tl(v)  dv  =  1.0645  t(vc)(^)  Avd_   n_36 


The  Doppler  temperature  is  given  by 


M     ,    n2 
D  "  8k  In  2  ^AV 
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Doppler  broadening  includes  the  effects  of  thermal  motion 

by  an  amount  AvT,  turbulence  of  a  scale  smaller  than  the 

antenna  beam  (Av,   ,  ),  and  expansion  or  contraction  of  the 
turb 

emitting  region  (av  ).   These  velocity  components  combine 
to  give  a  total  Doppler  width  of 


2  2 

Avn    =     (AvT    +    Av 


a    2x1/2 

i      +    Av 
turb  x 


11-31 
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Pressure  broadening  and  the  Voigt  profile.   When 
the  electron  density   is  high  enough  that  the  time  be- 
tween collisions  is  shorter  than  the  life-time  of  the 
excited  atomic  level,  the  energy  distribution  of  the 
radiative  transition  is  altered.  The  wavetrain  is 
abruptly  terminated  by  the  passing  electron.   The  line 
profile  emitted  by  the  assembly  of  atoms  is  a  Lorentz 
function  (Kuhn,  1962,  p.  392).   Peach  (1972)  has  shown, 
however,  that  the  profile  is  not  an  exact  Lorentz  func- 
tion, but  it  is  close  enough  for  an  accuracy  of  about  10 
percent.   The  main  effect  of  the  Lorentzian  is  to  redis- 
tribute the  emitted  energy  to  the  wings  of  the  line. 
Griem  (1967)  has  shown  that  proton-atom  collisions  are 
not  important. 

Brockelhurst  and  Seaton  (1972)  have  determined 
that  the  pressure  broadened  half-power  width  for  a  m  •*■    n 
transition  is 


Avp  =  9-4  W       {^r-] 


Pressure  broadening  will  occur  in  addition  to  Doppler 
broadening.   If  the  gas  has  an  electron  temperature  T 
and  a  Doppler  temperature  Tn  (Equation  11-37),  the  ratio 
of  pressure  to  Doppler  broadening  is,  by  Brockelhurst 
and  Seaton  (1972), 
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Av 

Av~ 


0.142 


— ) 

100; 


7.4 


io4 


2xloV/2  JO4,0"1 
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To  obtain  appreciable  pressure  broadening  of  radio  re- 
combination lines,  TD  >  2  x  10  ,  and  Tg  >  10  ,  and  Ng 
>  10  .   So  far,  the  only  possible  radio  recombination  line 
sources  known  to  have  such  parameters  are  dense  knots  in 
emission  nebulae  such  as  M42 ,  W51 ,  and  W3  (Kjellming  and 
Davies,  1970)  or  dense  H  II  regions  associated  with  type 
I  OH  masers  (Habing  et  al . ,  1974).   No  pressure  broadening 
has  been  seen,  however  (Gordon,  1974). 

When  Doppler  and  pressure  broadening  occur  to- 
gether, the  resulting  profile  is  the  convolution  of  a 
Lorentz  with  a  Gaussian  which  is  a  Voigt  profile,  4> .,  ( v  ) . 
The  convolution  integral  cannot  be  evaluated  exactly.  So, 
the  function  is  expressed  in  the  form 


<j>w(v) 


1 


\T  '    Av 


D  v/T 


H(a,oj)  . 
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where  H(a,w)  is  the  Voigt  function 

2 


H(a,co)  = 


•y 


(or-y)2  +  a2 


dy 


and 


y  = 


and  a 


Av^ 

a 
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The  frequency  dispersion,  o,  is  given  by  o  =  vD  2 ( 1 n  2 
=  vD  (— )  (In  2)~1/2.  Posener  (1959)  has  tabulated  the 
Voigt  function  for  values  of  astronomical  interest. 
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Non-LTE  Solutions  of  the  Transfer  Equation  and  Line 
Enhancement 


Two  solutions  of  the  transfer  equation  will  be 
developed  here  for  use  in  the  analysis  of  the  spectra. 
They  will  be  derived  starting  with  Equation  1 1 -4.   Since 
the  frequency  of  the  spectral  line  is  high,  the  Rayleigh- 
Jeans  approximation  to  the  Planck  function  can  be  used  to 
convert  specific  intensity  to  brightness  temperature. 
From  Equation  1 1  -  4 ,  the  brightness  temperature  is 

L, 


TB(v)  =  TL(v)  +  Tc(v)  -  A(v)  Te   (l-e"L)  +  f  (1-e   L)  TQe 
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where     t  =    t    (  v)    +    tt(v  )  . 


and    jl(v)    -    bm 


Since,  from  Equation  1 1  -  6  ,  A(v) 

* 
L 


S(v)/B(v,Te) 


v,T  ), 
e 


A(v 


KT  +  b    K. 

!     m   L 


1  +  bm  (KL/KT) 


:T  +  b   3    k.    l  +  b   3    (<  /<  )  ' 
T     n   mn   L        n   mn    L   T 
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Optically  thin  H  I  region  adjacent  to  a  strong 
source  of  radiation.   Assuming  that  «,>>  kj  at  the  line 
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center      (i.e.    a  t  v  =  v    )  , 


a(v; 


bm  <  J  V 


b     B       (k.  /  kt)     b     3 
n     inn   v   L'      J         n     mn 


Since    the    region    is    optically    thin,    Equation    11-41    re- 
duces   to 


Tn(vl 


bn  emn 


Te   (tt^l)   +Tn  +  Tq   (1-tt-tl). 


When    the    frequency    is    off    the    line    (i.e.v>>vc)      only   the 
continuum    is    observed; 


TB<v>-Tc(v)-rT-T.TT  +  VTo<lJtT> 

n     mn 


Subtracting  this  equation  from  the  preceding  one  gives  only 
the  line, 


TL(v)  -TB(v)  -TC(v)  =  b-f-TeTL-ToxL. 

n  mn 


Since    iL=bn    Bmn    t*, 


Tl(v)=tl(v)      (bmTe-bnBmnTo) 
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This  is  the  line  brightness  temperature  from  an  optically 
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thin  H  I  region  for  which  <L   >>  <T-   The  second  term  in 
Equation  11-43  will  produce  line  enhancement  if  3mn  is 
negative.   To  convert  the  brightness  temperature  to  an- 
tenna temperature,  the  right  side  of  the  equation  must 
be  multiplied  by  the  beam  dilution  factor,  pg,  and  the 
main  antenna  beam  efficiency,  n g .   These  are  defined  by 
PB  e  sis/QM  and  nB  e  QM/^A,  where  9^    is  the  solid  angle  of 
the  line-emitting  region;  qm  is  the  solid  angle  of  the 
antenna  main  beam;  and  fift  is  the  total  solid  angle  of 
the  antenna  beam.   The  antenna  temperatures  are  denoted 
with  a  prime;  so,  Equation  11-43  becomes 


TL(v)    - 


Tin      Tl/  (  v)      (  b         T 

r|B      Kv    ;    v    m      e 


b       3 


T    ) 
mn      o 
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Isolated  region.   This  solution  of  the  transfer 
equation  will  be  for  the  case  of  an  isolated  emitting 
region  for  which  the  incident  radiation  is  weak  and  can 
be  ignored.   Then  T   =  TM  =  0  and  Equation  11-41  reduces 

J  O      IN 


to 


-(t,+t  ) 
TD(v)  -  A(v)  T   [1-e      T  ]. 


Off  the  line,  A(v)  will  be  unity;  so,  the  continuum  tem- 
perature is 


Tc(v)  =  Te(l-e   T). 
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Taking  A(v)  from  Equation  11-42  and  converting  Tg(v) 
to  antenna  temperature,  TL(v)  =  TB(v)  -  T"c(v)  and 


TL(v)  =  PBnBTe 


J m L__ 

tt  +   b     3       T, 
T         n     nm     L 


-(t.+t    )  -tt 

[1-e       L     T  ]   -   (1-e     '; 
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Line  enhancement .   Equations  11-44  and  -45  will 
be  used  for  the  analysis  of  this  study's  observations. 
The  equations  indicate  how  non-LTE  conditions  affect  the 
radio  recombination  line  temperature.   For  LTE  bn  =  bm 


0;  and  Equation  11-44  becomes  T.  (v 


=  t,  (T  -T 
L    e   o 


Fiqure  1 1  -  4  shows  that  as  N   decreases,  collisions 

°  e 

become  less  frequent;  so,  the  level  populations  and  b^ 
decreases.   As  the  departure  from  LTE  increases,  6mn  be- 
comes larger.   When  6   >  (hv/kT),  6mn  is  negative  and 
amplification  of  the  line  takes  place.   This  process  is 
often  referred  to  as  the  "partial  maser"  process.   The 
populations   of  the  transition  levels  are  not  actually 
inverted  as  is  apparently  true  for,  say,  H^O  masers 
(Sullivan,  1973). 

The  line  enhancement  can  be  seen  from  Equations  II 
44  and  -45  in  which  tl  =  bn  3mn  tl  is  negative.   Usually, 

•k 

t.  <<   tt  because  the  path  length  involved  with  t  T 
is  much  larger  than  that  involved  in 
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t.  .   Also,  b   is  always  positive.   Therefore,  if   b 
L         m  in 

3   t, /tt   <  1,  A(v)  increases  and  the  total  absorption, 
mn   L   T I 

t,  +  t-t,  decreases . 

The  line  enhancement  will  be  lessened  somewhat 

because  b   <  1.   However,  for  a  decrease  in  electron 
m 

3        -3 
density  from  10  to  10  cm  ',5    increases  by  more  than  a 

J  mn 

factor  of  ten  while  b   decreases  by  about  14  percent. 

The  net  result  is  a  substantial  increase  in  line  power. 

Equation  11-44  also  indicates  the  possibility 

that  the  recombination  lines  may  be  seen  in  absorption 

against  a  strong  radio  source  if  the  relative  magnitudes 

of  b   and  6   are  such  to  make  3    positive.   However  even 
n      mn  mn 

for  the  case  of  a  dense  H  II  region  (E.  ^  10   pc  cm   ) 
seen  through  a  tenuous  H  II  region,  Brockelhurst  and 
Seaton  (1972)  have  shown  graphically  that  absorption 
cannot  be  greater  than  8  percent.   So  far,  radio  re- 
combination lines  have  not  been  seen  in  absorption. 


The  Radiation  Field 


The  radiation  field  intensity  entering  a  recom- 
bination line-emission  region,  and  used  in  Equation  11-20 
to  determine  the  level  popul ations,  i s  computed  from  IQ(v) 

=  W(v)  B(v,T).   Usually  the  radiation  source  is  a  nearby 

4 
nebula  whose  blackbody  temperature  is  10   °K.   The  dilu- 
tion factor,  W(v),  is  a  product  of  three  factors:, 
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W(v)  =  WQ  WA(v)  A(v 


11-46 


W.(v)  represents  absorption  of  the  radiation  by  the  in- 
tervening medium  between  the  source  and  the  line  emission 
region.   If  t  is  the  intervening  optical  depth,  w/\(v) 
=  e"T.   If  the  nebula  is  adjacent  to  the  line-emission 
region,  W. (v )  =  1 . 

A(v)  represents  the  departure  of  the  radiation 
source  from  a  blackbody.   It  is  defined  in  Equation  II-6. 

The  geometrical  dilution  factor,  Wg,  can  be 

found  as  follows.   Consider  a  radiation  source  which  is 

roughly  spherical  in  shape  with  a  diameter  dg  and  located 

at  a  distance  rR  from  the  line-emitter.   The  geometry 

is  shown  in  Figure  1 1  -5a .   The  line-emitting  region  sees 

the  radiation  source  as  a  disc  of  angular  diameter 

2 


and  filling  a  solid  angle  Q 


it  eR/4.      Then 


tt  d^/4 
4    7T    rl 


16    r 


IK. 

4tt     ' 


11-47 


If  the  two  regions  are  contiguous,  rR  -  1/2  dg 
and  W~  =  1/4.   If,  however,  the  radiation  source  is  much 

b 

larger  than  is  the  line  emitter,  Figure  1 1  -  5  b  applies  and 

W„  =  1/2.   Finally,  if  the  source  of  radiation  surrounds 
G 

the  line  emitter   or  if  the  line  is  emitted  by  a  dense 
H  II  region,  W g  =  1 . 
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Figure  1 1  - 5 .   Determination  of  the  Geometrical  Dilution 
Factor,  Wq.   Part  a  shows  the  line  emit- 
ting region  and  the  continuum  source  to 
be  about  the  same  size  and  rR  >  ds.   Part 
b  shows  the  line  emitting  region  adjacent 
to  a  large  continuum  source. 
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A  numerical  example  will  illustrate  the  case 

where  rD  >  d  /2.   Consider  the  radiation  of  temperature 
K     S 

TR  entering  a  small  line  emitting  region  located  100  pc 
from  an  H  II  region,  such  as  NGC  7000  (the  North  American 
Nebula),  for  which  T  =  10   °K  and  r   =  25  pc  (Kaplan  and 
Pikelner,  1970,  p.  152).   For  this  case  A(v)  =  1.   Assume 
that  W.(v)  =  1  also.   Then,  using  Equation  11-47,  it  is 
found  that  TD  =  131  °K.   This  indicates  that,  unless  the 
line-emitting  region  is  directly  adjacent  to  a  large  H  II 
region,  the  line-emitting  region  can  be  considered  as 
being  isolated.   In  the  case  where  the  line  emitter  is 
itself  an  H  II  region  (T   =  104  °K),  Dyson  (1969)  shows 
that  the  effect  of  the  free-free  continuum  on  the  departure 
parameters  is  very  small.   Dupree  (1972)  and  Hoang-Binh 
and  Walmsey  (1974)  demonstrate  the  effect  of  a  large  H  II 
region  adjacent  to  the  1 i ne- emi tti ng  region. 

Ionization  Rates 

Knowing  the  radiative  recombination  coefficient, 
a(T),  the  required  rate  of  ionization  per  atom  of  species 
X,  £v»  can  be  determined  from  the  condition  of  detailed 

A 

balance.   For  hydrogen, 


NH  cH  ■  Ne  NH+  a(T) 
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The  radiative  recombination  coefficient  may  be  obtained 
from  Equation  11-22.   Habing  (1970)  has  a  formula  for 
computing  c H  which  uses  the  assumption  that  a(T)  a  T 


CH  =  1.27  x  TO"15  v^1  T"°-3  tt(vt)/  JT21(V)  dv 
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In  this  formula,  tt(vt)  is  the  continuum  optical  depth 
at  a  frequency  vy  and  t2,(v)  is  the  21  cm  line  optical 
depth  discussed  in  the  next  section. 

Ionization  rates  can  be  calculated  theoretically 
from  a  given  source  of  ionization  such  as  low  energy  cos- 
mic rays  (Dalgarno  and  McCray,  1972)  and  OB  stars.   These 
rates  can  be  compared  to  the  ionization  rate  required  by 
a  model  for  a  line  emitting  region. 

As  an  example  the  ionization  rate  within  the 

Stromgren  sphere  surrounding  an  OB  star  will  be  considered 

Let  S   be  the  total  number  of  ultraviolet  photons  with 
u 

energy  greater  than  hv,  emitted  from  the  central  star, 
v,  is  the  frequency  of  the  Lyman  limit.   Also  let  r<.  be 
the  Stromgren  radius  which  is  the  radial  distance  from  the 
star  at  which  the  ionizing  radiation  flux  is  zero.   Then 


Su  *  J*  rS  Ne  V  a(T) 
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(Spitzer,  1968,  p.  116).   From  Equation  11-48,  a(T) 

=  N..  Cm/M   N..+  and  the  ionization  rate  of  an  OB  star  is 
H   H'  e   HT 


3S 


H    4  tt  r3  N 
4  *  rS   H 
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The  ultraviolet  flux,  S  ,  from  stars  of  spectral  class  05 
through  Bl  are  tabulated  by  Spitzer  (1968,  p.  117). 


The  21  cm  Hydrogen  Line 


The  well  known  21  cm  hydrogen  line  is  emitted  by  a 
magnetic  dipole  transition  between  the  hyperfine  levels 
of  the  ground  state.   This  is  the  result  of  the  difference 
of  energy  levels  between  the  proton  and  electron  spins 
being  parallel  (upper  state,  F  =  1)  and  antiparallel 
(lower  state,  F  =  0).   The  frequency  of  the  transition 

is  1420.405  MHz.   The  upper  state  is  metastable;  the 

-1  5 
spontaneous  emission  probability  is  A,„  =  2.85  x  10 

sec"1  (Kaplan  and  Pikelner,  1970,  p.  410).   Usually  the 
metastable  state  is  depopulated  by  collisional  excitation 
more  frequently  than  it  is  by  the  21  cm  radiative  transi- 
tion.  The  F  =  1  state  is  populated  by  atom-atom  col- 
lisions (Kaplan  and  Pikelner,  1970,  p.  420). 

Designating  the  lower  and  upper  states  by  0  and 
1  respectively,  the  absorption  coefficient  is 
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hv 


10 


Kio(^  =  TT^(v)  [No  Boi  "  Ni  Bio] 


The  parameters  are  analagous  to  those  of  the  recombination 
line.   Since  collisions  play  an  important  part  in  popu- 
lating the  upper  state,  it  will  be  assumed  for  the  moment 
that  the  populations  are  LTE  values  as  determined  by  the 
Boltzmann  equation.   The  statistical  weights  of  the  states 
are  g,  =  3  and  gQ  =  1 .   Then,  since  gQ  BQ1  -  g]  B1Q  and 
hv-,  Q/kT  <<  1  , 


10 


!v) 


hVl0    *t»\    hVl°    N      B 

(v)  FT"  No  Boi 

s 


4  TT 


T   is  the  excitation  temperature  which,  for  the  21  cm 
s 

transition,  is  called  the  spin  temperature  because  the 
transition  results  from  a  change  of  electron  spin  orien- 
tation.  The  absorption  coefficient  is  expressed  in 
terms  of  the  spontaneous  emission  coefficient,  A,q5 

3 


2h  v 


10 


'10 


10 


Then 


3N   h  c' 
o 


10 


(v)  =  o- "  t  -  -  Aln  <(.(v 


tt  k  T  v  iri  '10 
s  10 


11-51 


The  absorption  coefficient  is  used  to  obtain  the 
neutral  hydrogen  column  density  along  a  line  of  sight. 
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For  a  cold  gas,  the  total  neutral  hydrogen  density,  NHO 

is  Nn  +  N,  .   From  the  Boltzmann  equation,  the  number  den 
sity  at  some  point  s  in  the  line  of  sight  is 


NH(s)    =    NQ(s)    [1+3(1 -hv1Q/k    Ts)]    -    4NQ(s 


since  hv,n/k  =  0.07.   Substituting  in  values  for  the 
constants  and  integrating  over  the  line  profile  and  along 
the  line  of  sight  gives 


k1q(v,s)  ds  dv  =/  x21(v)  dv  =  2.58  x  10"ID   I  -^ ds: 
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where  the  path  length,  L,  is  in  cm  and  v , Q  is  in  Hz. 

If  the  1 i ne  prof i 1 e  is  Gaussian  with  a  half-power 


w 


idthv   and  a  peak  optical  depth  t 9 ,  ,  Equation  11-36  can 


be  used  to  obtain 
rL  Nu(s) 


21 


,14 


£ ds  =  4.126  x  10'"  t21  AvD 
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It  may  be  necesary  to  determine  whether  or  not  the 
spin  and  gas  kinetic  temperatures  are  equal.   Starting 
with  the  equation  of  detailed  balance, 
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Nl    NH    C10(V    +    Nl    A10    +    N1     B10    ^lO'V 


l0    NH    CQ1(TK)    +    NQ    BQ1     I(v10fTR) 


Cm ( T„ )  and  C,n(T„)  are  the  collisional  excitation  and  de- 
(J  I    K         I  U   In 

excitation  probabilities,  respectively,  whose  formulation 
is  given  by  Kaplan  and  Pikelner  (1970,  p.  420).   Solving 
the  detailed  balance  equation  for  T   in  terms  of  T..  and 


the  radiation  temperature,  TR,  gives 


557.3  T°-5  NH  +  TR 
557.3  T~0,5  NH  +  1 
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T  =  T.,  if  the  hydrogen  density  is  large  enough  that 
S      K. 

557.3  T°'5  NH  >>  TR.   For  H  I  regions,  TK  >  64  °K  and 

O 

N   >  1  cm   is  not  unreasonable  (Baker  and  Burton,  1975); 
e 

so  557.3  T?'5  Nu  >  4.4  x  103  °K  which  is  much  larger 
K     n 

than  one  would  expect  TR  to  be.   Therefore,  unless  NH 

<<  1  cm"3,  T   -  T„. 
s      l\ 

Absorption  spectra  of  the  21  cm  line  have  been 
obtained  by  three  techniques.   Hagen,  Lilley,  and  McClain 
(1955)  obtained  absorption  spectra  in  the  direction  of 
discrete  sources  by  averaging  the  spectra  from  off-source 
beams  and  subtracting  this  from  the  on-source  beam  spectrum 
The  accuracy  to  which  this  method  yields  the  true  absorp- 
tion spectrum  depends  on  how  uniform  the  hydrogen  gas  is 
over  the  grid  of  beams. 
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Radhakri  shnan  ej^  aj_.  (1972a)  and  Hughes,  Thompson, 
and  Colvin  (1971)  have  observed  absorption  spectra  in  the 
direction  of  radio  sources  (e.g.  3C391  )  with  two-element 
interferometers.   The  front-ends  and  synchronus  detectors 
are  switched  in  such  a  way  that  the  system  is  insensitive 
to  all  but  the  absorption  spectrum  in  the  direction  of  the 
source . 

A  third  method  has  been  used  to  obtain  absorption 
spectra  in  the  direction  of  pulsars.   Some  of  this  work  has 
been  done  by  Gordon  e_t  aJL  (1969),  Manchester  et  a_l_.  (1969): 
Gordon  and  Gordon  (1970),  Guelin  ejt  a]_.  (1969),  and 
Encrenaz  and  Guelin  (1970).   The  pulsar  is  observed  with 
a  gated,  spectral  line  receiver  which  alternately  stores 
and  integrates  the  spectrum  with  the  pulsar  on  and  with 
it  off.   At  the  end  of  the  integration  period,  the  pulsar- 
off  spectrum  is  subtracted  from  the  pulsar-on  spectrum. 
The  difference  of  the  two  is  the  absorption  spectrum.   The 
pulsar  is  a  particularly  useful  beacon  because  its  dis- 
persion measure  provides  a  value  of  the  integrated  electron 
density  along  the  line  of  sight.   Therefore,  two  pieces  of 
information  about  the  intervening  interstellar  medium  are 
known . 

Pulsar  Dispersion  Measure 


The  arrival  time  of  the  pulsed  radio  emission  from 
pulsars  is  frequency  dependent.   The  pulse  travels  at  the 
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group  velocity.   The  group  index  of  refraction  is  a  func- 
tion of  the  electron  density  of  the  medium  through  which 
the  pulse  travels.   The  difference  of  arrival  time,  t,  -  t? 
of  the  respective  frequencies  v ,  and  v  „  is,  then,  the  in- 
tegrated electron  density  along  the  path  length  L.   This  is 
the  dispersion  measure,  DM,  which  is  calculated  from  the 
equation  (cf.  Terzian,  1972) 


DM  e    f  Ne  ds  =  2.4  x  10-4  ( t -,  - 1 2 )  (v~2-v~2)-1 


11-55 


From  this  equation,  the  unit  of  DM  is  pc  cm 


-3 


CHAPTER  III 

THE  TELESCOPE  SENSITIVITY  AND 
THE  OBSERVING  SCHEME 


The  observations,  which  will  be  discussed  in  the 
next  chapter,  were  made  in  such  a  way  as  to  make  the  most 
efficient  use  of  the  telescope.   This  chapter  will  discuss 
the  basis  for  the  observational  scheme.   The  important  fac- 
tor was  the  system  sensitivity  which  is  affected  by  the 
microwave  frontend  of  the  receiver  and  the  one-bit  auto- 
correlation spectrometer. 

The  general  theory  involved  with  the  operation  of 
a  one-bit  autocorrelation  spectrometer  for  astronomical  use 
was  discussed  in  detail  by  Weinreb  (1963).   The  operation 
of  the  instrument  used  for  the  observations  is  discussed 
in  two  reports.   The  IF  filter  system,  which  precedes  the 
correlator  and  determines  the  width  of  the  spectral  window, 
is  described  by  Mauzy  (1968).   The  operation  of  the  entire 
system,  with  emphasis  on  the  digital  autocorrel ator ,  is 
described  by  Shalloway  et  al .  (1968).   A  discussion  of 
system  sensitivity  for  the  different  methods  of  observing 
radio  spectra  is  presented  by  Ball  (1972). 
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System  Sensitivity 


The  nominal  system  sensitivity  is  the  RMS  noise 
temperature,  AT,  given  by 


Y  T 


AT  = 


sys 


/Af  t. 


III-l 


The  parameters  are  as  follows: 

Y  =  a  constant  which  depends  on  the  system  and  the 
method  of  its  operation; 

T    =  the  system  noise  temperature; 
sys 

Af  =  the  frequency  resolution  bandwith;  and 

t .  =  the  integration  time. 
1 

T     is  the  combination  of  the  noise  temperatures 
sys 

received  by  the  antenna  as  well  as  the  power  generated  by 
the  antenna-transmission  line  system  and  the  receiver. 

The  total  bandwith,  B,  observed  by  the  spectrom- 
eter  is  divided  up  into  N  channels.   Both  B  and  N  are 
selectable.   Each  channel  has  the  effect  of  measuring  the 
spectral  power  within  a  filter  whose  bandpass  is  of  the 
form  (sin  v)/v  having  a  half-power  width  of 


Af  =  1.21  J 


1 1 1  -  2 


(Shalloway  e_t  aj_.  ,  1968).   The  spacing  between  channel 
centers  is  B / N . 
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The  factor  y    in  Equation  III-l  is  greater  than  one 
and  is  the  degradation  of  the  system  sensitivity  due  to  the 
operation  of  the  autocorrelator.   It  is  a  one-bit  digital 
machine  which  computes  the  autocorrelation  function,  p(nAt), 
of  the  undetected  IF  power;  nAt  is  the  time  delay  at  the 
n   channel.   The  first  step  in  obtaining  the  autocorrela- 
tion function  is  to  clip  the  signal,  a  random  Gaussian  volt- 
age, v(t).   The  clipper  generates  a  function,  y  ( t )  ;  y(t)  = 
+Y  if  v(t)  >  0  and  y(t)  =  -Y  if  v(t)  <  0.   The  function 
y(t)  is  then  sampled  at  a  rate  equal  to  2B  =  1/At.   The  out- 
put of  the  sampler  is  fed  into  an  N-bit  shift  register  and 
N  one-bit  multipliers.   The  shift  register  is  updated  at 
the  time  intervals  of  At.   The  output  of  each  of  the  multi- 
pliers updates  a  counter.   At  the  end  of  the  integration 
time,  each  counter  contains  the  sum  which  represents  one 
point  of  the  autocorrelation  function.   The  function  is  in 
the  form  of  an  array  which  is  transferred  to  the  on-line 
computer  which  computes  the  power  spectrum  by  taking  the 
Fourier  transform  of  it. 

The  effect  of  using  one-bit  quantization  is  to 
increase  AT  above  that  which  would  result  from  a  many-bit 
quantization  (Weinreb,  1963,  p.  32).   This  is  accounted 
for  in  the  constant  y    of  Equation  III-l.   This  constant 
is  actually  a  product  of  two  factors;  y    =  a$.   The  ef- 
fect of  one-bit  quantization  is  that  B  =  1.39  (Weinreb, 
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1963,  p.  49).    The  value  of  the  other  factor  depends  on 
the  weighting  function  used  in  computing  the  correlation 
function.   A  uniform  weighting  function  was  used  for  which 
a  =  1.099  (Weinreb,  1963,  p.  22).   Therefore,  y    =  1.528 
and  Equation  III-l  becomes 


1 .528  T 


AT 


sys 


/Af  t. 


1 1 1  -  3 


The  Observing  Scheme 

There  are  several  observing  schemes  which  can  be 
used  to  observe  spectral  lines  (cf.  Ball,  1972).   The  ob- 
serving scheme  must  provide  a  reference  spectrum  which  is 
used  to  remove  the  effects  of  the  receiver  passband  from 
the  des  i  red  spectra  . 

The  "Total  Power"  scheme  can  provide  an  efficient 
means  of  obtaining  the  reference  spectrum.   Before  observ- 
ing a  source  on  the  program,  a  spectrum  is  obtained  while 
the  antenna  is  pointed  in  a  direction  close  to  the  program 


Burns  and  Yao  (1969)  have  determined  that  the  clip- 
ping factor,  3,  is  a  function  of  the  sampling  rate.   For 
a  sampling  rate  of  2B,  3  =  tt / 2  =  1.571.   If  this  is  cor- 
rect, y  =  1-726  and  the  RMS  noise  of  the  data  discussed  in 
Chapter  IV  should  be  multiplied  by  a  factor  of  1.13.   How- 
ever, the  data  reduction  routines  continue  to  use  3  =  1-39 
(Cram,  1973,  p.  133). 
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source  and  where  the  spectral  line  is  not  detected.   This 
spectrum  is  the  receiver  passband  function.   The  nomen- 
clature of  the  Total  Power  method  is  to  designate  the  ref- 
erence spectrum  by  "OFF"  and  a  program  spectrum  by  "ON." 
If  the  receiver  passband  remains  unchanged  while  the  OFF 
and  ON  spectra  are  observed,  subtracting  the  OFF  spectrum 
from  the  ON  spectrum  removes  the  effect  of  the  passband 
from  the  desired  spectrum.   The  efficiency  of  this  observ- 
ing scheme  is  realized  if  one  OFF  spectrum  can  be  used  with 
more  than  one  ON  spectrum  if  there  are  several  sources  in 
the  observing  program. 

Let  the  observing  time  on  an  OFF  source  be  tr  and 
that  on  an  ON  be  tM.   If  an  OFF  source  is  followed  by  n 

IN 

ON  sources,  each  of  which  is  observed  for  a  time  t^,  the 

total  observing  time,  t  ,  will  be  tQ  =  tp  +  n  t[N|.   The 

times  tF  and  tN  and  the  ratio  n  must  be  chosen  so  that 

the  receiver  passband  remains  unchanged  during  the  time 

t  .   The  most  efficient   use  of  observing  time  for  a  given 
o 

number  of  program  sources  will  be  obtained  if  the  ratio 

a  =  tr/t.,  is  chosen  to  minimize  AT. 
r   IN 

When  an  OFF  spectrum  with  an  RMS  noise  (AT)  is 
subtracted  from  an  ON  spectrum  with  an  RMS  noise  (AT)N  the 
RMS  noise  of  the  resulting  spectrum  will  be 


T  =  [(AT)jj  +  (AT)*] 


1/2 


69 


For  spectral  line  observations,  T    is  not  expected  to 

b  y  o 

change  more  than  a  few  percent  between  the  OFF  and  ON  po- 
sitions; then  an  average  system  temperature,  T  ,  can  be 
used  to  get 


/Af    vrN    V/ 


Referring  to  Equation  III-l,  the  integration  time  is 


t . 
i 


tN  tF 


1 1 1-4 


Since  t  F  =  a  t N  and  t   =  (a  +  n ) t N , 


AT 


ave 


1/2 


/Af  t 


a  +  1  +  n  +  n/a 


The  value  of  a  is  determined  when   T  is  a  minimum,  for  which 


dun  =  o 

da 


?  -1/2 

Then  (1 -n/a  )  (a  +  1+n+n/a)   '   =  0.   The  only  reasonable 

solution  is  a  =  /7T  .   Therefore,  the  most  efficient  Total 

Power  observing  pattern  exists  for 


tF  =  /n  tN 


1 1 1  -  5 


CHAPTER  IV 
THE  RADIO  RECOMBINATION  LINE  OBSERVATIONS 

Introduction 

The  H157a  radio  recombination  line  spectra  for 
this  study  were  obtained  between  18  and  24  April,  1973 
with  the  43  meter  radio  telescope  of  the  National  Radio 
Astronomy  Observatory  (NRAO)   located  in  Green  Bank, 
West  Virginia.   The  rest  frequency  of  the  line  is  1.683200 
GHz.   The  H157a  line  was  chosen  for  two  reasons.   It  is 
conveniently  located  within  the  passband  of  the  NRAO  18 
cm  cooled  receiver  which  has  good  sensitivity  and  a  rep- 
utation of  being  very  stable.   Also,  the  frequency  of  this 
line  is  such  that  the  line  strength  is  relatively  strong. 
The  total  power  observing  scheme  was  used  to  obtain  the 
spectra . 

The  sources  used  for  the  observations  are  the  two 
supernova  remnants  (SNR)  3C391  and  3C396  and  the  pulsar 
1858+03.   These  were  chosen  on  the  basis  of  the  objectives 
discussed  in  Chapter  I;  they  are  located  in  the  galactic 


NRAO  is  operated  by  Associated  Universities,  Inc 


under  contact  with  the  National  Science  Foundation. 
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ridge  and  independent  data  from  the  intervening  inter- 
stellar medi urn  are  avai 1 abl e  .   A  radio  galaxy,  3C295  (Gold- 
stein 1962),  was  also  observed,  but  its  galactic  lati- 
tude is  61°  and  no  radio  recombination  line  was  observed 
in  its  direction.   The  purpose  of  observing  it  was  to  dem- 
onstrate that  the  recombination  lines  came  from  the  ridge 
gas  and  were  not  affected  by  instrumental  effects. 

Along  with  3C391  and  3C396,  a  region  adjacent  to 
each  of  them  was  also  observed.   These  will  be  referred 
to  by  the  designation  COMP.   The  COMP  spectra  were  observed 
for  two  reasons.   The  primary  reason  was  to  be  able  to 
subtract  out  emission  from  the  region  immediately  surround- 
ing each  of  the  SNR's  and,  thereby,  make  these  ob- 
servations compatible  with  those  made  with  narrow  beams. 
The  second  reason  for  obtaining  the  COMP  spectra  was  to 
determine  the  angular  scale  over  which  the  medium  changes. 

Table  I V  - 1  lists  the  center  position  of  each  beam. 
The  right  ascension  and  declination  (1950.0  epoch)  and 
the  galactic  longitude  and  latitude  are  given  for  each  of 
the  regions  studied.   The  OFF  beams  are  the  reference 
beams  used  for  the  Total  Power  observing  scheme.   Their 
selection  and  use  will  be  discussed  in  the  following  sec- 
ti  on . 

Figure  I V  - 1  shows  the  continuumradiation  environ- 
ment of  the  regions  observed.   The  cont i nuum  maps  are  adapted 
from  the  1414  MHz  maps  of  Altenhoff  et  aj_-  (1970).   The 
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size  of  the  circle  at  each  region  is  about  equal  to  the  half- 
power  width  of  the  antenna  beam,  18  arc  minutes  at  the 
H157a  frequency. 

Observing  Scheme  and  Data  Reduction 


The  receiver  front-end  was  mounted  at  the  prime 
focus  of  the  43  meter  paraboloid.   A  scaler  horn  was  lo- 
cated at  the  focal  point  of  the  reflector.   The  signals 
of  linear  orthogonal  polarizations  were  separated  behind 
the  horn  and  fed  into  separate  and  identical  receivers. 
The  box  containing  the  front-end  equ i pment  was  rotated  so 
that  the  E-plane  of  the  horn  was  aligned  normal  to  the 
celestial  equater.   Each  receiver  was  load-switched  and 
contained  a  two-stage  parametric  amplifier.   The  first 
stage  of  the  parametric  amplifier  and  the  switched  load 
were  refrigerated. 

The  local  oscillator  frequency  was  generated  by 
a  frequency  synthesizer  followed  by  a  6x  frequency  multi- 
plier.  The  intermediate  frequency  (IF)  sent  to  the  con- 
trol room  from  the  front-end  box  was  150  MHz.   The  synthe- 
sizer frequency  was  chosen  to  compensate  for  the  Doppler 
effects  of  all  motions  of  the  telescope  with  respect  to 
the  recombination  line  emitter  and  so  that  the  velocities 
of  the  spectra  would  be  with  respect  to  the  local  standard 
of  rest  (LSR).   Since  it  was  expected  that  the  recombina- 
tion line  emission  regions  and  galactic  H  1  regions  are 
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ki nemati cal ly  related,  the  21  cm  emission  velocity  lon- 
gitude maps  of  Westerhout  (1969)  were  consulted  to  deter- 
mine the  velocity  offset   needed  to  positon  the  antici- 
pated H157a  and  C157a  spectral  features  in  the  2.5  MHz 
spectrometer  passband.   The  velocity  offset   for  the  pul- 
sar was  +50  km/s.   All  other  program  regions  involved  an 
off  set  of  -10  km/s.   The  synthesizer  frequencies  were  cal- 
culated by  the  NRAO  computer  program  DOPSET  before  any  ob- 
servations were  made. 

The  150  MHz  IF  from  each  receiver  was  monitored  by 
a  continuum  reciever.   The  remaining  signal  was  fed  to 
the  NRAO  Model  II  autocorrelation  spectrometer.   The  spec- 
trometer creates  413  frequency  channels  which  were  divided 
into  three  sections.   One  section  of  192  channels  was  con- 
nected to  each  of  the  front-end  recei  vers  .   The  remaining  29- 
channel  section  of  the  spectrometer  was  connected  to  one 
of  the  recievers  to  serve  as  a  monitor.   It  was  not  used 
for  the  acquisition  of  the  data. 

Each  section  of  the  spectrometer  had  a  passband 
of  2.5  MHz.   This  allowed  a  channel  separation  of  13.0 
kHz.   Each  channel  has  a  sin  v/v  profile  with  a  half-pow- 
er width  of  15.8  kHz.   At  a  frequency  of  1.6832  Ghz,  the 
velocity  resolution  of  each  channel  was  2.82  km/s. 

The  operation  of  the  spectrometer  system  is  de- 
scribed by  Shalloway  e_t  aj_.  (1968)  and  the  Total  Power  ob- 
serving scheme  is  described  in  Chapter  III.   The  autocorrelator 
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periodically  transfers  the  autocorrelation  function  array 
into  the  on-line  computer  which  performs  the  Fourier  trans- 
form and  loads  the  resulting  spectrum  array  onto  a  seven- 
track  magnetic  tape.   This  array,  along  with  source  coor- 
dinates and  control  information,  constitutes  a  record.   The 
control  information  contains  a  number  designation  of  the 
spectrum  and  whether  it  is  a  spectrum  of  a  program  (ON) 
source  or  a  reference  (OFF)  source,   The  identity  of  its 
OFF  spectrum  is  also  specified. 

Each  source  is  observed  for  a  length  of  time  ( t  F 
or  t N  as  used  in  Chapter  III)  during  which  a  sequence  of 
records  is  placed  on  the  tape.   This  sequence  of  records 
is  treated  as  a  unit  during  off-line  analysis   and  is 
called  a  "scan."   At  the  beginning  of  the  off-line  pro- 
cessing, the  records  of  each  scan  are  averaged  together 
and  the  spectrum  of  each  scan  is  loaded  on  another  mag- 
netic tape.   Then,  for  each  ON  spectrum,  its  OFF  spectrum 
is  located  and  the  antenna  temperature  spectrum  is  com- 
puted from 

T.m  _  /S(i)  -  R(ih  T 
TL(i)  -  (    ^ )  Tsys  , 

where  S(i)  and  R(i)  are  the  source  (ON)  and  reference  (OFF) 
spectra,  respectively,  as  a  function  of  the  channel,  i. 

The  system  temperature,  T    ,  is  computed  by  the  off- 
line  computer  on  the  basis  of  a  calibration  noise  source 
whose  power  is  periodically  injected  into  each  receiver. 


The  accuracy  to  which  the  system  temperature  is  determined 
depends  on  how  well  the  calibration  temperature  is  known. 
During  the  observations  the  values  determined  during  bench 
testing  of  the  receiver  were  entered  into  the  DDP-116  pro- 
gram.  The  calibration  signal  was  injected  into  each  re- 
ceiver by  a  loop  "antenna"  placed  in  the  waveguide  at  a 
point  directly  following  the  coupler  which  seperates  the 
two  polarizations.   During  the  bench  testing  the  horn  was 
replaced  by  a  termination.   In  spite  of  the  careful  match- 
ing of  the  termination  to  the  system,  its  characteristics 
are  somewhat  different  from  those  of  the  horn  observing 
the  sky.   Consequently  the  amount  of  calibration  noise 
power  entering  the  system  is  not  exactly  the  same  in  the 
two  situations  (Turner,  1973). 

The  calibration  temperature  discrepancy  was  cor- 
rected by  daily  calibrations  using  the  radio  source  Virgo 
A  whose  continum  spectrum  is  well  known.   Once  a  more  ac- 
curate calibration  temperature  was  known,  a  correction  fac- 
tor was  later  applied  in  the  off-line  processing.   The  con- 
tinum flux  of  Virgo  A  at  1.6832  GHz  is  182  Jy  which  was 
determined  by  plotting  the  fluxes  published  by  Kellerman 
e_t  aj_.  ( 1  969  )  and  interpolating  to  this  frequency.   The  con- 
tinum receiver  connected  to  the  IF  was  used  to  measure  the 
relative  lev  Is  of  Virgo  A  and  the  cold  sky  with  and  with- 
out the  calibration  signal.   Table  IV-2  shows  the  results 
of  the  daily  calibration  by  which  the  system  noise  temperature, 
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T    »  was  also  determined.   The  first  calibration  on  22 
sys 

April  gave  unexpl ai nabl e  bad  results;  so,  the  calibration 
was  repeated  later.   The  spectra  obtained  in  the  meantime 
showed  no  indication  of  a  problem.   The  values  obtained 
during  this  bad  calibration  were  not  used  to  compute  the 
averages  from  which  the  correction  factors  were  computed. 

Before  the  telescope  was  pointed  toward  any  source, 
the  receiver  frontends  were  tuned  so  that  the  spectro- 
meter passband  was  located  in  a  portion  of  the  reciever 
passband  (25  MHz  wide  at  the  3  db  points)  which  was  rea- 
sonably flat.   The  daily  observing  procedure  began  by  peak- 
ing up  on  Virgo  A  in  order  to  check  the  pointing  accuracy 
of  the  telescope  and  obtain  a  correction  factor  for  the 
calibration  temperature.   No  problems  with  the  pointing 
were  found . 

Before  observing  the  first  program  region,  a  five- 
minute  OFF  scan  followed  by  a  five-minute  ON  scan  of  Ml  7 
was   made  to  check  the  operation  of  the  spectrometer.   The 
OFF  regions  for  this  and  the  program  regions  were  at  least 
four  degrees  above  or  below  the  galactic  plane  where  it 
was  certain  that  the  H157a  line  was  not  present.   The  H175a 
spectrum  of  M17,  an  emission  nebula,  is  very  strong. 
A  five-minute  scan  gave  adequate  indication  that  the  spec- 
trometer was  operating  properly.   The  daily  M17  spectra 
were  compared  to  look  for  degradation  of  the  system.   None 
was  noticed. 
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Table  I  V- 2 
Calibration  and  System  Noise  Temperatures 


Date  T    D 
sys  B 

Apri  1 

18  62.2 

19  66.4 

20  64.5 

21  61 .9 
(22a  176.9 

22  60.4 
2  3          84.5 

Average        66.6 
"Bench" val ues 
Correction  factor: 


cal  A 

3.12 
3.45 
3.  35 
3.02 
8.57 
3.21 
4.40 
3.42 
3.23 
0.975 


T 


cal  B 


58 

5 

59 

4 

60 

5 

59 

6 

187 

0 

58 

0 

84 

4 

63.5 


3 

17 

3 

98 

3 

55 

3 

75 

7 

83) 

3 

21 

3 

59 

3.54 
3.63 
0.972 


These  are  the  results  of  the  bad  calibration  and 
are  shown  for  reference  only.  They  were  not  used  to  com- 
pute the  averages . 


Most  of  the  observing  program  for  3C391  ,  its  COMP, 
3C396,  and  its  COMP  was  performed  with  the  Total  Power 
scheme  using  a  sixty-minute  cycle;  a  twenty-minute  OFF 
scan  was  followed  by  a  ten-minute  scan  on  each  of  the 
four  program  regions  in  succession.   The  ratio  of  OFF- 
time  to  ON-time  was  determined  on  the  basis  of  Equation 
1 1 1-5.   To  make  the  most  efficient  use  of  observing  time, 
it  is  desirable  to  make  the  scans  as  long  as  possible  so 
that  the  total  time  moving  between  regions  is  minimized. 
The  maximum  length  of  the  cycle  was  limited  by  the  length  of 
time  during  which  the  receiver  passband  remained  stable. 
During  the  observations  it  was  found  that  the  quotient  spec- 
tra of  the  last  few  records  of  the  sixty-minute  cycle  be- 
gan to  show  a  curvature  of  the  baseline.   This  means  that 
the  receiver  passband  had  changed  enough  since  the  OFF  scan 
had  been  recorded  that  this  OFF  spectrum  no  longer  repre- 
sented the  current  passband  well  enough  to  produce  a  flat 
continuum  on  which  the  spectral  features  are  enhanced.   The 
curvature  only  affected  the  last  few  3C396  COMP  spectra; 
so,  the  cycle  was  maintained  at  sixty  minutes. 

The  observations  of  these  four  regions  continued 
until  the  RMS  noise,  expected  after  combining  all  scans  of 
each  region,  was  about  0.002°K.   Most  of  the  last  two  days 
were  used  to  observe  the  pulsar  and  3C295.   Each  was  ob- 
served individually,  usually  with  a  fifty-minute  cycle 


82 


consisting  of  a  twenty-five  minute  OFF  followed  by  a  twen- 
ty-five minute  ON.   The  3C295  spectra  showed  no  spectral 
features  or  abnormal  baselines.   Therefore,  these  data  will 
not  be  discussed  further  in  this  study. 

After  each  day's  observations  had  been  completed, 
the  seven-track  "telescope"  tape  was  taken  to  Charlottes- 
ville, Virginia  where  an  IBM  360/55  computer  was  used  to 
average  the  records  of  each  scan  and  load  the  resulting 
spectra  on  a  nine-track  "user"  tape  for  further  analysis. 
At  the  same  time  a  summary  of  records  and  the  quotient  spec- 
trum of  each  scan  was  printed.   These  tasks  were  performed 
by  the  first  phase  of  the  TPOWER  data  reduction  program 
(Cram, 1973).   The  printout  was  sent  to  the  observatory  the 
next  day.   With  these  "quick-look"  summaries,  it  was  possi- 
ble to  determine  which  scans  had  problems  which  would  make 
them  unusable.   On  the  basis  of  these  reports  some  scans 
were  deleted  and  reobserved  on  the  last  day  of  the  program. 

After  all  observations  had  been  made,  the  final 
spectra  were  produced  using  the  second  phase  of  TPOWER,  a 
system  of  routines  controlled  by  a  command  language.   It 
was  during  this  process  that  the  correction  factor,  de- 
termined from  the  daily  noise  tube  calibration,  was  applied 
to  adjust  the  antenna  temperature  scale.   The  scans  were 
summed  and  weighted  by  the  effective  integrating  time  of 
each  scan.   Then  the  spectra  of  the  two  receivers  were  av- 
eraged together,  thereby  doubling  the  integrating  time.   The 
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resulting  spectra  were  convolved  with  a  five-channel  "box- 
car" function  and  plotted.   The  plotted  spectra  were  stud- 
ied in  order  to  estimate  the  structure  of  the  baseline. 
Curvature  of  the  baseline  of  each  spectrum  was  removed  with 
a  TPOWER  routine  which  first  computes  a  specified  order 
Chebyshev  polynomial  by  means  of  a  least  squares  fit  to 
a  specified  portion  of  the  baseline.   The  routine  then 
subtracts  the  polynomial  baseline  from  each  spectrum  and  the 
spectra  are    then  plotted  in  their  final  form. 

A  final  task  was  performed  in  order  to  facilitate 
the  analysis  of  the  spectra.   Gaussian  functions  were  fitted 
to  the  spectral  features.   The  TPOWER  routine  begins  with 
initial  estimates  of  the  center  channel  and  half- intensity 
width  of  each  feature  and,  from  that,  determines  a  best- 
fit  center,  half- intensity  width,  and  central  intensity  of 
the  features . 


The  Spectra 


The  basic  spectra  for  the  program  regions  are  shown 
in  Figures  I V  -  2 ,  -3,  and  -4.    Preceding  the  figures  is 
Table  I V- 3  which  lists  the  total  effective  integration  time 
(t.)  and  the  order  of  the  Chebyshev   polynomial  base  line 
which  was  removed.   The  effective  integration  time  was  com- 
puted with  Equation  III-4.   It  includes  the  doubling  ob- 
tained by  averaging  the  spectra  of  the  two  receivers. 
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The  RMS  noise,  listed   in  the  third  column,  is  com- 
puted from  Equation  1 1 1  -  3  ; 

1 .528  T 


AT  = 


All 


3ti 


T     is  the  average  system  temperature  of  both  receivers 
sys 

and  is  65.0  °K.   The  effective  bandwidth,  8,  is  five  times 
the  channel  separation  of  13.0  kHz  because  the  spectra  are 
displayed  after  being  convolved  with  a  five  channel  box- 
car function.   Table  I V  -  3  also  provides  a  table  of  contents 
which  includes  a  list  of  the  figure  numbers  for  all  spectra 
including  those  showing  the  Gaussian  fits  to  the  spectral 
features . 

The  lower  portions  of  Figures  IV- 5  through  10  show 
the  results  of  fitting  Gauss  in  functions  through  the  spec- 
tral features.   The  Gaussian  functions  of  the  individual 
features  are  added  together  and  the  combined  function  is 
plotted  through  the  original  spectrum.   The  upper  portions 
of  these  figures  show  the  plots  of  the  residuals  which  are 
obtained  by  subtracting  the  combined  Gaussian  function  from 
the  spectrum. 

The  two  possible  features  at  -118.8  and  58.2  km/s 
on  the  COMP  3C396  spectrum  are  only  about  four  times  the 
RMS.   Therefore,  it  was  decided  not  to  present  Gaussian 
features  fit  to  them. 

The  center  velocity,  half-power  width,  peak  antenna 
temperature  and  integrated  line  power  of  each  feature  in  the 
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spectra  are  listed  in  Table  I V- 4  which  follows  the  spectra 
Also  listed  for  each  parameter  is  its  standard  error  (stan- 
dard deviation  of  the  mean)  which,  except  for  the  standard 
error  of  the  integrated  line  power,  is   computed  by  the 
Gaussian  fitting  routine.   The  standard  deviation  of  the 
line  power,  sp,  is  computed  from  the  following  equation 
(Parratt  1961  ,  p.  117)  with  Equation  11-35: 


sp  =  [(gy)   sT   + 


3P 

9Av 


*l   ],/2 

AvJ 


=  1.0645  [(Av)2  s2  +  T2s2v]1/2. 

Here,  sT  is  the  standard  error  of  the  peak  line  temperature 

and  s,   is  the  standard  error  of  the  line  half-power  width. 

Av 

The  discussion  of  the  spectra  will  be  a  part  of  the 
analysis  of  the  emitting  region.   Chapter  V  will  present 
the  analysis  of  the  spectrum  from  the  directions  of  3C391 
and  its  COMP  beam.   The  model  derived  there  will  be  applied 
to  the  analysis  of  the  medium  in  the  direction  of  3C396, 
its  COMP  beam,  and  the  pulsar  1858+03. 


Figure  IV-2.  H157a  Spectra  from  the  Direction 
of  3C391  and  Its  COMP  Beam.  The 
velocity  resolution  is  11.6  km/s 
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Figure  IV-3.   H157a  Spectra  from  the  Direction  of 
3C396  and  Its  COMP  Beam.   The  veloc- 
ity resolution  is  11.6  km/s. 
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Figure  IV-5.   Gaussian  Profiles  Fitted  to  the 

3C391  Spectrum  with  a  Single  Gaus- 
sian Through  the  Hi gh- Vel oci ty 
Feature.   The  lower  curve  is  the 
fitted  spectrum.   The  upper  curve 
is  thatof  the  residuals. 
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Figure  IV-6.   Gaussian  Profiles  Fitted  to  the 

3C391  Spectrum  with  Two  Gaussians 
Through  the  High- Vel oci ty  Feature. 
The  lower  curve  is  the  fitted  spec- 
trum.  The  upper  curve  is  that  of 
the  residuals. 


0.06 


0.04 


0.0  2 


0.00 


0.06 


96 


I  < 1 1 1 1 1  I H 


0.04 


0.02 


0.00 


-170      -130      -90       -50       -10 


30         70  110        150 

v        km/s 


Figure  IV-7.   Gaussian  Profiles  Fitted  to  the 

3C391  COMP  Spectrum  with  a  Single 
Gaussian  Through  the  Hi gh- Vel oci ty 
Feature.   The  lower  curve  is  the 
fitted  spectrum.   The  lower  curve 
is  that  of  the  residuals. 
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Figure  IV-8.   Gaussian  Profiles  Fitted  to  the  3C391 
COMP  Spectrum  with  Two  Gaussians 
Through  the  High-Velocity  Feature. 
The  lower  curve  is  the  fitted  spec- 
trum.  The  upper  curve  is  that  of 
the  residuals. 
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Figure  I V-9 .   Gaussian  Profiles  Fitted  to  the  3C396 
Spectrum.   The  lower  curve  is  the 
fitted  spectrum.   The  upper  curve  is 
that  of  the  residuals. 
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Figure  IV-10.   Two  Gaussian  Profiles  Fitted  to  the 
Pulsar  1858+03  Spectrum.   The  lower 
curve  is  the  fitted  spectrum.   The 
upper  curve  is  that  of  the  residuals 
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CHAPTER  V 


MODELS  FOR  THE  INTERSTELLAR  MEDIUM  IN 
THE  DIRECTION  OF  3C391 


3C391  :   A  Summary  of  Previous  Observations 

The  radio  source  3C391  has  the  nonthermal  spectrum 
of  a  supernova  remnant  (SNR).   Caswell  et  aj_.  (1971) 
give  its  spectral  index  between  750  MHz  and  5  GHz  to  be 
-0.57.   However,  there  is  a  turnover  of  the  spectrum  and 
the  spectral  index  between  80  and  408  MHz  is  about  0.13. 
Low  frequency  optical  depths  derived  from  radio  recombina- 
tion lines  and  also  the  continuum  spectrum  itself  strongly 
indicate  that  the  low  frequency  decrease  in  the  flux  is 
due  to  absorption  by  the  interstellar  gas  rather  than  due 
to  a  process  within  the  SNR  itself  (Bridle  and  Kesteven, 
1971;  Gordon,  1972;  Dulk  and  Slee,  1972). 

Assuming  that  the  absorption  is  due  to  the  inter- 
stellar gas,  the  optical  depth  at  80  MHz  may  be  determined 
by  assuming  that  the  flux  from  the  SNR  follows  a  power  law 
in  the  frequency  range  between  80  MHz  and  5  GHz.   The  spec- 
tral index  is  constant  over  this  range  and  was  found  to  be 
-0.50  using  the  fluxes  measured  at  408  MHz  5  GHz.   The 
observed  flux  at  80  MHz,  S(80)  was  compared  to  that  expected 
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from  the  power  law,  S  (80).   For  3C391  , 


t(80)  -  In  |°^1  =  0.84  (±0.3) 


The  antenna  beamwidth  used  was  (axfi)  =  3.7  x  4.8  arc 
minutes  which  is  the  size  of  the  angle  subtended  by  3C391 . 

Caswell  et  aj_.  (1971)  compute  this  optical  depth  to 
be  1.3  assuming  a  spectral  index  of  -0.57  between  80  MHz 
and  5  GHz.   Goss  (1972)  makes  use  of  the  H109a  recombina- 
tion line  power  and  the  5  GHz  brightness  temperature  com- 
puted by  Jackson  and  Kerr  (1971)  to  compute  t(80)  for 
the  direction  1  =  31?2,  b  =  0°.   This  position  is  0?7 
away  from  3C391  ,  but  the  optical  depth  was  assumed  to  be 
valid  for  the  direction  of  3C391 .   For  an  electron  tempera- 
ture range  of  10?  to  104,  Goss  finds  that  x(80)  is  close 
to  one. 

From  the  observations  tabulated  by  Caswell  et  al  . 
(1971)  it  is  evident  that  the  size  of  3C391  does  not  change 
appreciably  with  frequency.   Their  Mills  Cross  observations 
at  408  MHz  give  its  size  as  (axfi)  =  3.2  x  3.4  arc  minutes. 
Their  drift  scan  map  shows  the  source  to  be  "an  intense 
compact  feature  within  a  weaker  more  extended  structure." 
Their  angular  resolution  is  (ax5)  =  2.06  x  3.47  arc 
minutes.   Bridle  and  Kesteven  (1971)  detect  the  same 
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structure  with  a  10.63  GHz  map.   Goss  (1972)  determined  the 
size  at  31.4  GHz  to  be  3.2  arc  minutes  in  diameter;  his 
resolution  was  2.8  arc  minutes.   A  map  made  by  Chaisson 
(1974a)  at  15.5  GHz  with  a  resolution  of  2.1  arc  minutes 
shows  an  arc  of  structure  which  is  generally  surrounded  by 
the  lowest  intensity  isophote  of  the  10.63  GHz  map  of 
Bridle  and  Kesteven  (1971)  which  they  claim  to  be  of  low 
r  e  1  i  a  b  i  1  i  ty  . 

The  distance  to  3C391  can  be  estimated  by  two 
methods,  from  an  empirical  relation  between  radio  surface 
brightness  and  from  21  cm  absorption  profiles.   The  21  cm 
absorption  profile  obtained  by  Radhakri  shnan  et  aj_.  (1972) 
(shown  in  Figure  V-l)  has  an  absorption  peak  with  a  radial 
velocity  of  107  km/s  with  respect  to  the  local  standard  of 
rest  (LSR).   This  is  only  slightly  higher  than  the  maximum 
velocity  of  105  km/s  predicted  by  the  Schmidt  rotation  model 
for  the  galaxy  (Schmidt,  1965)  and  observed  by  Burton 
(1970)  at  a  galactic  longitude  of  31?9.   Assuming  pure 
circular  galactic  rotation  of  the  gas,  the  maximum  velocity 
corresponds  to  emission  from  the  tangent  point  where  the 
line  of  sight  at  31?9  intersects  the  5.3  kpc  (radius)  orbit. 
Also,  it  is  assumed  that  the  sun  is  10  kpc  from  the  galac- 
tic center.   The  distance  from  the  sun  to  this  point  is 
8.5  kpc.   Since  the  21  cm  line  is  seen  in  absorption  against 
the  SNR  at  this  velocity,  it  can  be  concluded  that  the  SNR 
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lies  no  nearer  than  8.5  kpc.   Because  of  the  distance 
ambiguity  of  the  radial  velocity-distance  relation  (cf. 
Mihalas  and  Routly,  1968,  p.  126),  it  cannot  be  determined 
with  certainty  whether  the  other  absorption  features  lie  in 
front  of  the  tangent  point  or  behind  it.   Therefore,  all 
that  can  be  said  is  that  the  distance  to  3C391  is  greater 
than  8.5  kpc  and  may  be  as  close  as  13.4  kpc.   Ilovaisky 
and  Lequeux  (1972)  derive  a  distance-luminosity  function  for 
galactic  supernova  remnants.   They  use  3C391  as  a  calibra- 
tor for  which  the  distance  was  assumed  to  be  11  ±  2  kpc. 

There  is  considerable  disagreement  as  to  the  nature 
of  the  interstellar  medium  in  the  direction  of  3C391  .   The 
intent  of  this  study  is  to  provide  another  data  point  which 
can  be  used  with  previous  studies  to  generate  a  consistent 
model.   Of  particular  interest  are  the  recombination  line 
observations  of  Cesarsky  and  Cesarsky  (1973a, b)  who  ob- 
served wide  (^30  km/s)  H92a,  H157a,  and  H166a  lines  in  the 
direction  of  3C391 .   Using  their  H92a  line  and  x(80) 
measured  by  Dulk  and  Slee  (1972),  they  employ  a  unique 
technique  to  derive  an  electron  temperature  of  20°K  with 
the  gas  in  non-LTE.   Chaisson  (1974a)  takes  issue  with  this 
solution  and  presents  an  argument  in  support  of  the  low 
frequency  absorption  resulting  mainly  from  an  H  II  region 
imbedded  in  the  two-phase  medium  and  the  remaining  absorption 
from  the  partially  ionized  cool  component  of  the  medium. 
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The  wide  recombination  lines  are  indicative  of  thermal 
broadening  of  a  hot  emitter.   For  that  reason  Chaisson 
suggests  that  the  recombination  lines  are  emitted  by  the 
H  II  regions  from  which  most  of  the  absorption  results. 
Downes  and  Wilson  (1974)  also  conclude  that  the  line 
emission  can  come  from  H  II  regions  rather  than  from  H  I 
regions.   The  structure  seen  in  Chaisson's  map  could  pos- 
sibly be  evidence  for  there  being  one  or  more  H  II  regions 
in  the  intervening  medium.   However,  it  is  not  certain  that 
the  arc  of  sources  which  surrounds  the  east,  south,  and 
west  of  3C391  exist  in  the  medium.   Since  they  are  included 
within  an  isophote  of  Bridle  and  Kesteven's  (1971)  map, 
the  sources   form   a    structure  within  3C391  itself. 
Also,  both  Chaisson  (1974a)  and  Bridle  and  Kesteven  (1971) 
find  an  increase  i n  the  continuum  flux  at  high  frequencies 
above  10  GHz  over  that  expected  from  the  spectral  index  at 
lower  frequencies.   This  excess  flux,  they  suggest,  could 
be  due  to  at  least  one  H  II  region  imbedded  in  the  inter- 
vening medium.   Goss  (1972)  does  not,  however,  see  any 
evidence  for  an  excess  at  31.4  GHz  and  suggests  that  the 
low  frequency  absorption  and  radio  recombination  lines  are 
both  due  to  a  diffuse  medium  with  a  temperature  ranging 
from  102  to  104  °K. 

Figure  V-l  shows  the  21  cm  absorption  and  emission 
spectra  obtained  by  Radhakri  shnan  e_t  al_.  (1972)  placed 
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above  this  study's  H157a  spectrum.   Figure  V-2  shows  the 
recombination  line  spectra  which  the  Cesarskys  (1973a, b) 
observe  placed  above  that  of  this  study's.   The  Cesarskys' 
H147a  and  H166a  spectra  show  a  possible  feature  at  a 
velocity  of  -45  km/s  which  they  suggest  may  be  real  and  may 
be  the  C157a  line.   The  velocities  in  both  figures  are  with 
respect  to  the  local  standard  of  rest.   On  this  study's 
spectrum,  the  horizontal  bar  labeled  with  a  C  indicates  the 
range  of  velocities  in  which  the  C157a  line  would  be  if 
observed,  assuming  that  it  is  colocated  with  one  or  more  of 
the  hydrogen  line  emission  regions.   The  vertical  arrows 
labeled  with  CO  indicate  where  the  CO  features,  seen  by  the 
Dickels  (1974),  are  located. 

Downes  and  Wilson  (1974)  looked  for  but  did  not  de- 
tect H134a  and  C134a  recombination  lines.   Their  signal- 
to-noise  ratio  is  such  that  the  H134a  peak  temperature  is 
less  than  0.02  °K  and  that  of  the  C134a  is  less  than  0.03 
°K.   Since  T|_a   1/n,  these  limits  are  close  to  the  expected 
H134a  peak  temperature  and  probably  above  the  expected 
C134a  peak  temperature. 

The  work  which  has  been  discussed  is  summarized  in 
Table  V-l.  The  conclusion  which  one  might  draw  from  this 
discussion  is  that  there  are  three  possible  models  of  the 
interstellar  medium  —  at  least  in  the  direction  of  3C391  — 
which  must  be  considered:  a  group  of  cold,  partially 
ionized  H  I  regions,  one  or  more  H  II  nebulae,  or  an 


Figure  V-2.   Comparison  of  Radio  Recombination  Line 

Spectra  from  the  Direction  of  3C391.   The 
upper  two  spectra,  H166a  and  H157a,  are 
from  Cesarsky  and  Cesarsky  (1973a).   The 
H92a  spectrum  is  from  Cesarsky  and  Cesarsky 
(1973b).  The  bottom  spectrum  is  the  H157a 
spectrum  of  the  present  study.   All  velo- 
cities are  with  respect  to  the  local  stan- 
dard of  rest . 
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extended  diffuse  H  II  region.   The  remainder  of  this 
chapter  will  discuss  this  study's  observations  in  terms 
of  parameters  consistent  with  these  models. 

The  Models 

The  D-Method  Models 

This  is  the  method  used  by  the  Cesarskys  (1973b) 
in  obtaining  a  model  from  their  H92a  recombination  line 
spectrum.   They  began  by  assuming  a  transfer  equation 
solution  suitable  for  an  H  I  region  observed  in  front  of 
a  strong  continuum  source.   From  Equation  11-43,  which 
assumes  that  i^  and  tj<<1  ,  the  line  temperature  is 


T, (v)  =  x*(v)  (b  T  -b  6   T  ) . 
L v  '  L v  '  v  m  e   n  mn  R; 


The  radiation  temperature  is  the  sum  of  several  terms: 


TR  =  TT+TN+TSNR+2-7 


T$nr  1S  the  radiation  temperature  of  3C391  as  experienced 
by  the  line-emitting  region  and  Tj  and  TR  are  the  respec- 
tive temperatures  of  the  thermal  and  nonthermal  ambient 
radiation  fields  excl udi ng  that  of  3C391  .   The  2.7  °K 
term  represents  the  universal  thermal  background  radiation 
The  line  power  emitted  along  a  unit  pathlength,  As,  is 
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AP(v)  =  kl  (bmTe-bn3nmTR)As 


V-l 


The  line  absorption  coefficient  (Equation  11-34)  is 


k*   =   fK*(v)dv    =    fl.071xl07<f>(v)Ajl  fnm   Tg2-5    exp(en/kT)Ni    Nedv, 


Line 


Line 


where  N-  is  the  number  density  of  the  recombining  ion 
(hydrogen).   The  thermal  continuum  absorption  coefficient 
at  80  MHz  along  As  is  found  by  combining  Equations  1 1  -  7 
with  the  defintion  of  optical  depth; 


At(80)  =  4.710T"1-5  1 n (0. 6194  Tj-5)<NeEZ-  Nj>As.     V-2 


Dividing  Equations  11-28  and  V-l  by  t(80)  and  Ar(80) 
respectively  provides  two  functions,  D  and  DQ,  defined  as 
fol lows : 

nDP, 


o  "  tt(80 


so , 


a      a  (b  T     -  b  S     TD)   exp(e n/kT 

D     =  2.274xl04nRMf       1-n.e         n  mn  R_ .  i^_ZL_ 

o  B  n      mn 


)  E. 


ln(0.6194  T^-5) 


V-3 


and 


so , 


At(80) 
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a         a             (b  T     -  b  3     Tn)   exp(c  /kT)       <N  Nu>As. 
D  =  2.274xl04nR^f       iiSLS n  mn   R'       ^   n       L.       *"       L     -  ,V-4 


ln(0.6194  TgJ) 


<N  £Z.N.>AST 


The  antenna  main  beam  efficiency,  nB >  is  the  ratio  of  main 
beam  solid  angle  to  total  antenna  solid  angle  (Kraus, 
1966,  p.  73).   Neglecting  the  contribution  of  carbon  to  the 
electron  pool,  the  Cesarskys  use 


¥ j  Nj  ■  HH+  +  NHe  + 


1  .1  N 


H+  ' 


so,  ET  =  1.1  EL  if  the  absorption  and  line  emission  media 

are  colocated  and  their  respective  path  lengths,  LT  and 

L,,  are  equal.   If  C+  is  present,  the  abundance  is  NC=NC+  - 

4xl0"4  NH  (Spitzer,  1968,  p.  122).   Then  ZZ?  Nj«  1.1NH  is 

a  good  approximation  only  if  the  fractional  ionization, 

x  =  N  /N  ,  is  high.   Nu  is  the  total  hydrogen  density; 
e   H  H 


'H+ 


NH  -  NH0  H  N, 

This  technique  computes  D  for  models  specified  by 
choices  of  TR,  Ne,  and  Te.    For  each  model,  curves  of  D(Te) 
for  specific  values  of  N   and  TR  will  be  computed  and  com- 
pared to  the  observed  quantity,  DQ,  in  such  a  way  that 
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N   and  T  can  be  specified.   The  estimate  of  Tn  will  remain 
e      e         K  R 

the  same  for  all  models.   There  are  two  problems,  however. 
First,  there  is  one  "known,"  DQ ,   and  three  "unknowns," 
TR,  N  ,  and  T  .   This  is  complicated  by  the  fact  that  bn 

and  8   are  functions  of  these  variables  and,  particularly 

mn  r 

at  low  temperatures  (T  < 1 00  °K) ,  this  functional  relation- 
ship is  not  well  known.   There  is  controversy  regarding 
the  bound-bound  collision  cross-sections  at  these  tempera- 
tures which  correspond  to  kinetic  energies  near  the  reac- 
tion threshold  (Percival  and  Seaton,  1972;  Dupree,  1972b). 
Ne  is  not  completely  eliminated  in  Equation  V-4  since  bn 
and  3mn  are  nonanalytic  functions  of  it. 

The  second  problem  involves  the  rel ation  between  DQ 
and  D.   D   is  a  ratio  of  integrated  quantities;  D  is  a 
ratio  of  differentials.   The  Cesarskys  (1975b)  show  that, 
for  the  H92oc  line,  two  extreme  conditions  provide  limits 
for  D.   D(H92a)  will  be  largest  for  an  LTE  solution  for 
which  Te =  20  °K.   Since  for  this  case  bn  =  Bmn  =  l  ,  Ne  does 
not  have  to  be  specified  and  does  not  enter  into  the  solu- 
tion.  The  conditions  for  which  D  will  be  a  minimum  are  a 
non-LTE  situation  where  Te  =0  0  °K  and  Ne<10"3  cm-3.   At 
these  electron  densities  collisions  need  not  be  considered 
and  the  departure  parameters  are  better  known.   These  con- 
ditions dictate  that  6<DQ(H92a)  <  60  °K  kHz. 
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The  Cesarskys'  spectral  line  gives  DQ=6.8  °K  kHz 
from  which  they  conclude  that  the  most  probable  temperature 
of  the  emitter  is  20  °K  with  a  possible  range  of  tempera- 
ture of  0<T  <1000  °K  due  to  the  observational  error  of  70 
percent.   Substitution  of  Te=20  °K  into  Equation  11-51 
gives  F_,=<Ne  NH+>L=4  pc  cm"6. 

There  is  another  matter  which  should  be  considered. 
In  a  two-phase  medium,  the  line  emission  from  each  phase 
must  be  considered.   In  an  earlier  paper,  the  Cesarskys 
(1971)  showed  that  for  a  two-phase  medium  like  those 
proposed  by  Field,  Goldsmith,  and  Habing  (1969)  and 
Hjellming,  Gordon,  and  Gordon  (1969),  the  line  emission 
will  come  predominantly  from  the  cold  regions.   The  frac- 
tion of  the  line  power  from  the  hot  region  increases  as  Tg 
decreases  because  TL<*  T"2-5  N2.   However,  emission  from 
the  hot  component  will  generate  a  weak  and  wide  (^80  km/s) 
line  because  of  the  temperature  dependence  and  the  fact 
that  this  gas  is  widely  distributed  so  that  there  is  line 
broadening  resulting  from  differential  galactic  rotation 
along  the  line  of  sight.   Consequently  the  hot  region  spec- 
tral feature  will  be  difficult  to  identify  because  it  will 
be  immersed  in  the  baseline  noise  or  be  hidden  by  one  or 
more  features  from  the  cool  regions.   For  these  reasons, 
the  Cesarskys  (1975b)  consider  emission  only  from  a  cool 
region . 
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It  should  be  noted  that  the  antenna  beam  used  to 
obtain  this  study's  spectra  is  18  arc  minutes  in  diameter 
which  is  considerably  larger  than  Cesarskys'  beam  (2.3 
arc  minutes)  and  3C391  (3.5  arc  minutes).   It  will  be 
assumed  that  t(80)  is  constant  within  the  18  arc  minute 
beam  and,  therefore,  the  D-method  can  be  used  for  this 
study.   The  comparison  of  the  high-velocity  features  of 
the  3C391  and  its  COMP  spectra  (Figure  IV-2)  shows  a  vari- 
ation of  line  power  by  a  factor  of  2.9  over  a  distance  of 
17  arc  minutes.   The  42  km/s  feature  varies  oppositely  by 
a  factor  of  1.2  and  is  not  seen  in  the  H92a  beam  which  is 
2.3  arc  minutes  wide. 

Since  a  high-velocity  feature  appears  on  both  the 
3C391  and  its  COMP  beams,  it  seems  to  originate  from  an 
extended  region.   Using  this  feature  only,  a  solution  con- 
sistent with  the  H92a  observations  may  be  possible  using 
the  D-method.   However,  considering  that  the  solution  is 
expected  to  yield  a  low  temperature,  it  is  relevant  to  con- 
sider the  21  cm  hi gh- veloci ty  absorption  features.   There 
are  two  prominent  features,  one  at  about  80  km/s  and  another 
at  107  km/s.   This  raises  the  question  as  to  whether  the 
H157ct  high-velocity  feature  is  single,  centered  on  99.82 
km/s  and  having  a  half-power  width  of  35.7  km/s  or  is 
double  with  components  at  80.66  km/s  (Av  =  18.4  km/s)  and 
102.8  km/s  (Av  =  29.2  km/s).   It  is  not  certain  whether  or 
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not  this  feature  is  double  on  the  H92a  spectrum.   There 
may  be  a  component  at  about  80  km/s,  but  its  peak  is  no 
more  than  two  times  the  RMS  noise.   The  Cesarskys1  (1973a) 
H166a  and  H157a  spectra  show  some  structure  other  than  the 
103  km/s  feature  which  has  an  intensity  roughly  three 
times  the  RMS.   Taking  each  spectrum  by  itself,  it  is  not 
clear  that  there  should  be  a  feature  around  80  km/s.   Taken 
together,  the  spectra  indicate  a  strong  possibility  of 
there  being  a  feature  there. 

To  see  what  effect  the  choice  of  situation  has  on 
the  solution,  both  possibilities  will  be  considered.   As- 
suming that  the  feature  is  single  and  using  t(80)  from 
Dulk  and  Slee  (1972)  , 

Do(H157a)  =  o.ll  UO.If  =  1M4  (±2'89)  °K  kHz- 

Taking  only  the  103  km/s  feature  if  a  double  feature 
exists, 


Do(H157a)  =  HI   (±o:3)}  =  9-37  (±2'37>  °K  kHz' 


Assuming  that  the  uncertainty  quoted  for  t(80)  is  one 
standard  deviation,  the  standard  deviation  of  DQ  is  com- 
puted by 
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'D 


P2    s2    +    t(80)2      s2 

T  P 


1/2 


(Parratt,    1961  ,    p.    117) . 

For   the    H157a   line,    Equation    V-4    becomes 

D  =  3.83X105   (bmTe  -   bn3mnTR)   exp(6.4059/Te)   ^    E± 
T     In   (0.6194  T^5)  ET 


V-5 


T  and  Tn  are  brightness  temperatures.   D  and  D   are  in 

e  K  J  r  Q 

terms  of  antenna  temperature;  therefore,  the  constant 
multiplier  in  Equation  V-5  contains  the  factor  hb  =  °-8> 
the  beam  efficiency  of  the  43-meter  NRAO  telescope  at  the 
H157a  frequency  (Gottesman  and  Gordon,  1971). 

TR  may  be  computed  from  the  transfer  equation  solved 
for  the  situation  shown  in  Figure  V-3.   If  the  gas  between 
3C391  and  the  line  emitter  has  a  kinetic  temperature  Te 
and  continuum  optical  depth  tj,  then 


TR  =  Te  (1-e-  T)  +  W  TSNR  +  TN  +  2.7 


V-6 


The  dilution  factor,  W,  is  given  by 

l     d2 
U  -  A(v)  -i|4  exP(-TT) 

RR 


V-7 


(cf.  Equation  11-46) 
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LINE    EMITTING 
REGION 


OBSERVER 


Figure  V-3.   Geometry  for  Determining  TR.   This  is  a 
schematic  drawing  and  is  not  to  scale. 
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Ilovaisky  and  Lequeux  (1972)  determine  the  diameter 
of  3C391  to  be  11  pc.   For  this  discussion,  it  can  be 
assumed  that  TqNR  =  104  °k  and  the  gas  between  3C391  and  the 
line  emitter  has  the  characteristics  of  the  intercloud 
medium  of  the  Fi eld,  Gol dsmi th ,  and  Habing  (1969)  two-phase 
model;  then,  Te  -  7500  °K  and  Ne  =  0.012  cm"3.   Also  it 
will  be  assumed  that  Ne  =  NH+,  ET  =  1.1<N2>RR  =  1 . 6x1 0"4  RR. 
Since  3C391  is  between  8.5  and  13.4  kpc  from  the  sun 
(Caswell  et  a]_.  ,  1971),  the  range  on  RR  is  0<RR<4.9  kpc. 

W  can  be  calculated  directly  from  a  flux  measurement. 
The  flux,  S  ,  measured  at  1667  MHz  is  16.4  (±1.5)  Jy  (Cas- 
well et  al_.  ,  1971).   The  COMP  beam  is  in  the  ridge;  so, 
subtracting  its  observed  continuum  temperature  from  that 
observed  in  the  3C391  beam  should  provide  a  good  approxi- 
mation to  the  radiation  temperature  of  3C391  itself.   From 
Table  IV-3,  T^(3C391 )  =    10.7  °K  and  T^(COMP)  =  6.1  °k. 
Their  difference  is  4.6  °K. 

The  flux  measurements  give 


^S   c2 


C    B  nM  2kv 


5.7  °K, 


where  ^M  and  Q$  are  the  respective  solid  angles  of  the 
antenna  main  beam  and  3C391  as  seen  at  the  telescope. 
Therefore  the  flux  at  1667  MHz  provides  a  reasonable 
estimate  of  WTSNR  at  1683  MHz. 
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At  the  line  emitter  3C391  subtends  a  solid  angle  of 

2     2 
n+   which  is  equal  to  ird  /4  RR.   Then 


WT 


2  sc  n, 

c     S   ' 


SNR 


2kv2  QS  4^ 


The  following  table  shows  that  TR  ^  TN  +  2-7  except  when  the 
line  emitter  is  adjacent  to  the  SNR.   The  optical  depth  is 
computed  from  Equation  II-8.   Since  the  optical  depth  is  so 
small,  TD  is  not  attenuated  and  is  equal  to  the  observed 


*R 
[pc] 


Table  V-2 
Calculation  of  Radiation  Temperature 


LT 
Lpc  cm   J 


5.5 

0.009 

2.9x10 

550 

0.088 

1 .4x10 

-10 


W  T 


SNR 
[°K] 


60 
6x10' 


VTN 
[°K] 


17.7 
2.7 


continuum  brightness  temperature  and  Equation  V-6  reduces  to 


TN  +  2.7 


V-8 


Then  for  3C391,  TR  =  10.7/nB  =  13.4  °K.   This  makes  the  ther- 
mal component  20  percent  of  the  total  continuum  if  only  the 
2.7  °K  background  is  present.   This  is  midway  between  the 
percentages  computed  by  Gottesman  and  Gordon  (1970)  who 
argue  for  10  percent  and  Matthews,  Pedlar,  and  Davies  (1973) 
who  obtain  30  percent. 
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Figure  V-4  shows  the  plot  of  Equation  V-5  with  all  of 
the  parameters  which  have  been  drived  inthe  discussion.   Two 
extreme  cases  are  plotted.   The  solid  curve  is  the  LTE 


function  for  which  b>n  =  bm 


mn 


1 .   The  dashed  curve  i  s 


-4    -3 
the  col  1 i si onless  non-LTE  case  where  Ne<10   cm   .   The  non- 

LTE  departure  parameters  were  computed  by  Dupree  (1972). 

Their  values  are  listed  in  Table  V-3. 


Table  V-3 

Non-LTE  Departure  Parameters  with  N e< 1 0 
[from  Dupree  (1972)] 


Te  [oK] 


b  3 
n  mn 


10 

20 

100 

1000 


0.0583 
0.  1150 
0.2720 
0. 5350 


-0.0438 
-0. 1442 
-1.0414 
-8.7328 


The  values  of  D  ,  9.37±2.37  and  11.44±2.89,  are 
indicated  on  the  plot  by  labeled  horizontal  lines  with 
error  bars.   Considering  both  values  of  D  ,  it  is  seen  that 
the  solution  is  an  LTE  gas  with  a  temperature  of  1 4 . 3±0 . 3  °K. 

The  solution  disagrees  markedly  in  one  respect  from 
that  of  the  Cesarskys1  (1973b)  which  is  a  non-LTE  solution 
with  a  most-probable  temperature  of  20  CK.   To  obtain  this 
solution  for  this  study's  spectrum,  DQ  would  have  to  be 
25.3  °K  kHz,  a  value  well  outside  of  the  error  quoted.   The 
value  of  DQ  (H157a)  which  the  Cesarskys  predict  is  about 
25  °K  kHz. 
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If  t(80)  =  1.3,  as  derived  by  Caswell  et  al_.  (1971), 
the  H157a  solution  would  be  closer  to  TR,  but  the  H92a 
solution  would  yield  a  temperature  of  less  than  10  °K. 
This  is  evident  from  Cesarskys'  (1973b)  Figure  3.   Recent 
unpublished  observations  suggest  that  t(80)  may  be  close  to 
0.4  (Chaisson,  1974b).   If  this  is  the  correct  value, 
D  (H157a)  will  be  about  24.2  °K  kHz  which  will  give  a  non- 
LTE  solution  with  a  temperature  of  about  60  °K  if  P  is  7.87 
°K  kHz  and  an  indeterminate  solution  if  P  is  9.61  °K  kHz. 
However,  with  x(80)  =  0.4,  DQ(H92a)  will  be  14.2  °K  kHz 
giving  a  non-LTE  solution  with  a  temperature  of  about  1000  °K. 

Even  though  the  surrounding  hot  medium  will  generate 
a  wide,  shallow  line,  it  is  necessary  to  determine  what 
effect  this  component  will  have  on  the  solution.   For  this 
purpose,  Equation  V-5  can  be  modified  to 


n  -    AP(c)  +  AP(h)  v  q 

U    At (80,  c)  +  At (80,  h)  ' 


where,    letting    x    equal    to    c    for    the    cold    component   and    h    for 
the    hot    component, 


AP(x)   -   1.84X106  r2-5  exp(6.4049/Tx).(bm  T  -b  S^,  Tr)aELx 


and 


At(80,x)   -  4.710  T"1,5     ln(0.6194  T1'5)     AE,. 

X  XIX 
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Because  of  the  different  number  densities  in,  and  vastly 
different  path  lengths  through,  the  two  regions,  it  cannot  be 
assumed  that  their  emission  measures  will  be  equal.   There- 
fore, they  will  not  cancel  and  at  least  two  more  variables 
are  added  to  the  solution. 

To  solve  equation  V-9,  four  combinations  of  emission 
measures  will  be  chosen.   For  each  combination  Ej  =  1.1  E , 
for  both  components.   The  cold  region  non-LTE  departure 
parameters  are  those  computed  by  Dupree  (1972)  and  listed  in 

Table  V-3.   For  the  surrounding  hot  region,  Te  will  be 

_3 
7500  °K  and  N   is  assumed  to  be  near  0.01  cm   for  the 

purpose  of  computing  the  departure  parameters.   The  values 

of  T   and  N   correspond  to  those  of  the  stable,  hot  phase  of 

the  two-phase  model  of  Field,  Goldsmith,  and  Habing  (1969). 

_  3 
Also,  N   =  0.01  cm    is  analogous  to  the  col  1 i s ionless  case 

for  which  Dupree  (1972)  computed  the  "Ne  =  0"  parameters. 

Sejnowski  and  Hjellming  (1969)  have  computed  departure 

parameters  for  T   =  7500  °K,  but  in  order  to  obtain  them  for 

_3 
N   =  0.01  cm   ,  they  must  be  estimated  by  extrapolating  the 

radiative  asymptote  through  n  =  157  (cf.  Figure  II-5). 

With  this  extrapolation,  the  departure  parameters  for  the 

_3 

non-LTE  hot  medium  are  b   =  0.75  and  6    =  10   ,  which  makes 

m  mn 

b  3    =  -68. 95. 
n  mn 

For  each  of  the  models  examined,  D(H157a,  Te(c)) 
will  be  computed  for  each  of  the  possibilities: 
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(a)  cold  region  in  LTE,  hot  region  in  non-LTE 

(b)  cold  region  in  LTE,  hot  region  in  LTE 

(c)  cold  region  in  non-LTE,  hot  region  in  non-LTE 

(d)  cold  region  in  non-LTE,  hot  region  in  LTE. 

In  each  of  these  cases,  the  continuum  emission  measure  is 
assumed  to  be  related  to  the  electron  density  by  ET(c)  = 
Ne(c)2  pL  L  and  ET(h)  -  Ne(h)2  (1  -  pL)  L,  where  L  is  the 
total  line-of-sight  path  length  through  both  media  and   p. 
is  the  line-of-sight  filling  factor  of  the  cold  regions.   L 
will  be  assumed  to  be  equal  to  the  distance  to  3C391  al- 
though this  may  not  be  a  good  assumption  since  the  electron 
density  may  be  appreciably  different  in  the  ridge  gas  than 
it  is  in  the  solar  neighborhood.   From  Figure  I V - 2  it  is 
seen  that  a  portion  of  the  line  of  sight  does  not  lie  within 
the  ridge.   It  is  not  clear  what  the  ridge  boundary  is  like 

and  how  much  NQ  differs  in  and  out  of  the  ridge.   Pulsar 

e  3 

data  seem  to  indicate  that,  in  the  solar  neighborhood,  N   = 

_3 
0.03  cm    (Guelin,  1974).   The  error  in  choosing  L  =  8.5  kpc 

should  have  little  effect  on  the  solution  if  it  is  shown  that 

the  hot  region  does  not  contribute  much  to  the  recombination 

line  power. 

Model  1  is  essentially  the  same  as  that  previously 
examined;  it  is  based  on  the  Cesarskys'  (1973b)  model  for  a 
cold  region  emitter:   E>  (c)  =  4  pc  cm"   and  Te  =  20  °K. 

Model  2  is  based  on  Model  1  but  with  E.(c)  reduced 
by  a  factor  of  10  to  see  what  the  effect  of  doing  so  would 
be.   The  electron  density  remains  the  same  and   the  path 
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length  through  the  cold  region  becomes  85  pc.   If  the  region 
is  spherical,  it  would  subtend  an  angle  of  34  arc  minutes. 
This  is  a  reasonable  size  because  the  line  strength  observed 
by  the  COMP  beam  is  about  one-third  of  that  observed  by  the 
3C391  beam.   This  suggests  that  the  size  of  the  emitting 
region  may  be  about  34  arc  minutes  in  diameter. 

Model  3  is  based  entirely  on  the  results  of  Field, 
Goldsmith,  and  Habing  (1969).   A  filling  factor  of  P[_  =  0.006 
is  used  for  this  analysis  which  falls  into  the  range  they 
suggest.   The  cold  region  emission  measure  is  based  on  the 
information  given  in  their  Figures  1  and  2. 

Model  4  was  inspired  by  the  results  of  Chaisson's 
(1972)  recombination  line  observations  in  the  direction  of 
W3A,  and  H  II  region  located  at  1  =  134°  and  b  =  0°.   His 
spectrum  indicates  the  presence  of  hydrogen,  helium,  carbon, 
and  a  blend  of  heavy  element  recombination  lines.   The 
residuals  he  obtains  by  removing  Gaussian  features  from  his 
spectrum  suggest  the  possibility  of  an  additional  weak, 
narrow  component  of  the  prominent  hydrogen  feature.   He 
concludes  that  the  weak,  narrow  feature  originates  from  a 
cold  region  where  the  heavy  elements  (Mg,  Si,  S,  and  Fe) 
are  located.   The  prominent  feature,  then,  originates  from 
the  H  II  region.   Using  the  H  I  and  carbon  lines  and  assuming 
the  abundances  given  by  Spitzer  (1968,  p.  122),  he  obtains 

Te  =  100  °K,  NR  =  10  cm"3,  NR+  =  1.8  (±1.4)xl0~4  cm"3  and 

-3    -3 
N   =  7  ( ±  5 )  x 1 0    cm    for  the  cold  region.   This  gives  most 
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-5 


probable  values  of  N„NUi  =  1.26x10   cm   and  N   I  Z 
K  e  H+  e  l 


'e"H 
5    -6 


e  j   j   J 


1.1N   NUi  =  1.39x10   cm   .   Chaisson  does  not  suggest  a 
e   H  + 

value  for  p. ,  but  a  value  of  0.01  will  be  used  here  because 
of  the  argument  used  for  Model  2 . 

Figures  V-5,  6,  7,  and  8  show  the  D-curves  plotted 
for  each  model  and  combination  of  thermodynamic  condition. 
The  parameters  used  for  these  are  summarized  in  Table  V-4. 
D  was  computed  with  a  WATFIVE  computer  program. 

The  LTE  cold  curves,  labeled  a  and  b  in  each  figure, 
are  asymptotic  at  T   =  TR  =  13.4  °K.   As  the  plots  show, 
neither  DQ  =  9 . 37 (±2 . 37 ) °K  kHz  nor  DQ  -  11.44  (±2.89)  °K 
kHz  is  large  enough  to  obtain  a  non-LTE  solution  for  the 
cold  region.   The  result,  therefore,  remains  the  same  as 
that  for  Figure  V-4.   The  hot  region  has  no  effect  on  the 
solution  and  the  error  in  assuming  that  L  =  8.5  kpc  is  of 
no  consequence.   The  RMS  electron  densities  shown  in  Table 
V-3  are  computed  from 


<Ne(x)2>  ^2 


ET(x)/sx 


1/2 


where  sx  is  either  sc  =  p.L  or  sn  =  C  ~  PL)L-   Table  V-3 
also  shows  the  fraction,  in  percent,  of  the  line  power  which 
originates  from  each  of  the  two  phases.   It  is  apparent  that 
the  hot  medium  has  little  effect  on  the  line  power  except 
when  E  (c)  <<  E(h) . 
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Two  matters  must  be  considered.   First,  there  is  the 

question  of  the  validity  of  the  cold  region  departure 

<    -4 
parameters.   They  were  chosen  for  the  case  where  N   ^  10 

Because  of  the  proximity  of  the  neighboring  curves,  there  is 

some  tolerance;  that  is,  Ne  can  be  somewhat  larger  without 

affecting  the  departure  from  LTE.   The  amount  of  tolerance 

depends  on  n  and  T  .   At  n  =  157  and  Tg  =  10  °K,  the  electron 

-3    -3 
density  of  the  medium  should  not  exceed  10   cm   .   If  T  = 

_3 
1000  °K,  the  electron  density  could  be  as  large  as  0.1  cm 

The  error  from  using  col  1 i s i onl ess  values  for  a  medium  where 

_3 
N   is  actually  about  0.1  cm   — as  was  done  for  Models  1,  2, 

and  3  —  is  probably  not  large.   Percival  and  Seaton  (1972) 

claim  that  there  are  large  uncertainties  in  the  bound-bound 

collision  cross-sections.   For  that  reason  the  Cesarskys 

(1973b)  estimate  that  the  col  1 i s ionl es s  condition  may  be 

valid  at  electron  densities  which  are  60  times  that  expected 

from  Dupree's  (1972)  curves.    If  the  D-curve  for  Ng  =  0.1 

_3 
cm   were  used,  it  would  fall  between  curves  a  and  c  because 

collisions  will  be  more  frequent  and  bring  the  level  popula- 
tions closer  to  their  LTE  values.   Consequently,  the  D  -lines 
will  not  intersect  the  non-LTE  curves  and  the  solution  is 
not  affected. 

Finally,  there  is  the  question  as  to  the  use  of  the 
D-method  for  the  H157a  spectrum  since  the  antenna  beam  is 
much  larger  than  3C391 .   As  mentioned  before,  Figure  1  of 
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Cesarsky  and  Cesarsky  (1973b)  predicts  that  D   obtained 
with  an  antenna  beam  2.3  arc  minutes  in  diameter  should  be 
25  °K  kHz.   The  fact  that  the  18  arc  minute  beam  yields  a 
much  smaller  D  may  indicate  clumping  or  small  scale  struc- 
ture within  the  emitting  region  seen  by  the  antenna  beam. 
Because  of  beam  dilution  the  antenna  power  is  less  than  it 
would  be  if  the  line  emitting  region  filled  the  beam. 
Since  the  D-method  gives  uncertain  results,  another  approach 
will  be  taken . 

Isolated  Cold  Region  Models 

Since  the  background  radiation  is  quite  small,  a 
solution  to  the  transfer  equation  can  be  used  which  is  ap- 
plicable to  an  isolated  region.   The  solution  will  be  left 
in  terms  of  line  power  in  order  to  avoid  the  use  of  a 
parameter  with  a  large  uncertainty  such  as  x(80).   A  solu- 
tion which  is  consistent  between  the  H92a  and  the  H157a 
spectra  will  be  found  by  adjusting  emission  measures  until 
both  observed  line  powers  are  duplicated.   By  doing  this, 
both  T   and  E,  can  be  determined  if  there  is  a  consistent 


solution 


From  Equation  11-45,  the  peak  line  antenna  tempera 


t u r e  is 


TL(vc)  =  PB  nB  Te 


A(vJ  (1  -  e"(TL+TT))  =  (i 


e-TT) 


V-10 
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where , 


A(v, 


VVl 

VTT 


1+b  t. /tt 

1+b  B  nT,/TT 

n  m  n  L   T 


Both  optical  depths  are  peak  values  (i.e.,  at  v  -  vc)  and 

*  * 

t,  =  b  3  T.  .   The  optical  depths  t,  and  tt  are  computed 
L    n  mn  L        r  L      ' 

from  Equations  11-28  and  1 1  -  8  respectively.   Let  «L>  «M , 
and  n,    be  the  solid  angles  of  the  line  emitting  region, 
antenna  main  beam,  and  total  antenna  beam  respectively. 

Then  pb  =  ^J%    and  nB  =  VV    For  these  models  1t  is 
assumed  that  the  line  emitting  region  fills  the  H157a  and 
H92a  beams.  So,  in  both  cases  pB  =  1.   It  is  also  assumed 
that  ET  =  1.1  E..   The  integrated  line  power,  P,  in  °k  kHz, 
is  computed  by  the  use  of  Equation  V-10  and  Equation  11-35; 


P  -  1.0645  TL(vc)~  AvD 


In  computing  the  Hl57a  line  power,  a  single  high- 
velocity  Gaussian  feature  with  full  half-power  width  of 
35.73  km/s  (AvD  -  200.6  kHz)  is  used.   The  expected  H157a 
and  H92a  line  powers  are  computed  for  a  range  of  emission 
measures.   Both  LTE  and  col  1 i s i onl es s  non-LTE  functions  are 
computed.   The  departure  parameters  used  are  those  computed 
by  Dupree  (1972),  the  same  as  those  used  for  the  D-method 
and  tabulated  in  Table  V-3. 
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A  consistent  solution  between  this  study's  spectrum 
and  the  H92a  spectrum  is  found  graphically.   Figure  V-9 
shows  non-LTE  curves  and  Figure  V-10  shows  the  LTE  curves. 
Each  curve  is  the  integrated  line  power,  P,  plotted  as  a 
function  of  electron  temperature  and  is  labeled  by  the  value 
of  the  emission  measure  used  in  computing  it.   In  both 
figures,  the  H157acurves  are  placed  above  those  for  H92a  and 
the  horizontal  lines  represent  the  respective  observed  line 
powers:   PQ(H92a)  =  5.7  (±2.4)  °K  kHz  and  P  (H157a)  =  9.61 
(+0.26,  -1.27)  °K  kHz.   The  error  with  PQ(157a)  is  deter- 
mined by  considering  the  possibility  of  the  feature  being 
double.   The  error  bar  for  both  P  -lines  is  represented  in 
the  figures  by  shaded  strips. 

The  intersection  of  the  P  -lines  with  each  of  the 

o 

respective  curves  in  Figures  V-9  and  10  is  plotted  in  Figure 
V-ll.   Each  of  the  four  curves  in  Figure  V-ll  is,  then,  the 
locus  of  points  which  relate  the  temperature  and  emission 
measure  required  to  produce  the  observed  line  power.   Error 
strips  are  not  shown.   An  intersection  of  the  respective 
LTE  or  non-LTE  loci  defines  a  consistent  solution. 

The  LTE  curves  do  not  intersect;  therefore,  the 
emitting  medium  is  not  in  LTE.   The  non-LTE  loci  do,  however, 
intersect.   The  angle  between  them  is  small,  making  the 
solution  difficult  to  identify.   The  most  probable  values 
are  Te  =  53  (±6)  °K  and  EL  =  9.9  (±0.6)  pc  cm"6.   This 
solution  does  not  involve  the  errors  in  either  locus.   The 


Figure  V-9.   Non-LTE  Integrated  Line  Power  as  a  Function 
of  Te  for  an  Isolated  Cold  Region.   The 
upper  curves  are  for  the  H157a  line;  the 
lower  curves  are  for  the  H 9 2 a  line.   Each 
curve  is  labeled  with  the  emission  measure 
(in  pc  cm"6)  used  to  compute  it.   The  ob- 
served line  power  is  represented  by  the 
line  labeled  with  P0  and  the  shading  repre- 
sents the  error  bars. 


151 


100 


■    '  '  '    I 


HI57a 


* (— i — i — i  i  i  i  |  i 1 — * — • — i  >  ■  * 


50 


H92a 


\20  50 


1000 


Figure  V-10.   LTE  Integrated  Line  Power  as  a  Function  of 
Te  for  an  Isolated  Cold  Region.   The  upper 
curves  are  for  the  H157a  line;  the  lower 
curves  are  for  the  H92a  line.   Each  curve 
is  labeled  with  the  emission  measure  (in 
pc  cm -6)  used  to  compute  it.   The  observed 
line  power  is  represented  by  the  line 
labeled  with  P0  and  the  shading  represents 
the  error  bars . 


153 


1000 


154 


1000 


100 


10 


1 

■ ■ ' ' — ' 

. 

///' 

• 

/S                H920/X 

//-"HI  5  7a 

;               H92a         yS 

y^NON-LTE                  ; 

■ 

i 

• 

10 


20 


50  _ 
pc  cm 


100 


Figure  V - 11 .   Isolated  Cold  Region  Model  Solutions 
The  consistent  solution  is  indicated 
by  the  intersection  of  the  lines. 
There  is  no  LTE  solution.   The  non- 
LTE  solution  is  Te  =  53°K  and  El  =  9 
pc  cm-6. 
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error  strip  for  each  locus  would  encompass  both  lines  and 
render  an  uncertainty  to  the  solution  equal  to  the  range  of 
the  loci.   Hence  the  solution  presented  is  only  a  best 
estimate. 

This  solution  temperature  does  not  seem  to  agree 
very  well  with  the  most  probable  temperature  derived  by  the 
Cesarskys  (1973b)  but  is  closer  to  the  21  cm  spin  temperature 
of  60  °K  derived  by  Radha kr i shnan  and  Goss  (1972).   But 
considering  the  uncertainties,  the  agreement  is  all  right. 
This  result,  along  with  the  alignment  of  the  high-velocity 
recombination  line  features  with  the  84  and  105  km/s  21  cm 
absorption  features,  tends  to  indicate  that  their  regions  of 
origin  are  spatially  related. 

The  beam  dilution  was  taken  to  be  unity  for  both 
lines  because  the  variation  of  the  H157a  line  power  between 
its  two  beams  seems  to  indicate  that  the  emitting  region 
fills  both  of  them.   This  suggests  a  H157a  emission  region 
size  of  about  85  pc .   This,  however,  is  quite  a  bit  larger 
than  the  average  size  of  the  21  cm  absorption  regions  (5  pc) 
derived  by  Baker  and  Burton  (1975).   If  85  pc  is  equal  to 
the  total  path  length  through  the  recombination  line  emis- 
sion region,  the  emission  measure  of  9.9  pc  cm   derived 

_3 
here  results  in  an  RMS  electron  density  of  0.34  cm   .   It 

is  possible  that  the  recombination  lines  may  come  from  a 

_3 
cluster  of  5  pc  clouds.   Since  Nu>0.34  cm   ,  Equation  11-54 

H 

shows  that  the  spin  and  electron  kinetic  temperatures  are 
equal  . 
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It  is  difficult  to  explain  the  width  of  the  recombi- 
nation line  features  if  the  emitter  has  a  temperature  of 
53  °K.   Since  53  °K  is  the  gas  kinetic  temperature,  the  ther- 
mal line  broadening  should  be  1.6  km/s  (Equation  11-34). 
This  is  obviously  much  smaller  than  the  observed  line  width 
which  is  between  30  and  36  km/s.   Turbulence  within  the  cold 
emission  region  will  cause  some  broadening,  but  the  average 
total  half-intensity  width  of  21  cm  absorption  features  seen 
in  the  galactic  plane  is  about  3  km/s;  this  includes  all 
broadening  effects  (Radhakrishnan  and  Goss,  1972).   The 
widths  of  the  80  and  107  km/s  features  observed  in  the  direc- 
tion of  3C391  are  10  and  3  km/s  respectively.   The  80  km/s 
feature  may  be  a  blend  of  two  or  more  features  (Radhakrishnan 
and  Goss,  1972)  . 

The  discrepancy  between  the  21  cm  absorption  and 
radio  recombination  line  widths  suggests  considering  other 
model s . 


The  Hot  Spot  Model s 


The  observational  results  of  Matthews,  Pedlar,  and 
Davies  (1973),  Gottesman  and  Gordon  (1970),  and  Chaisson 
(1974a)  suggest  that  the  recombination  line  emission  regions 
must  be  hotter  than  1000  °K.   One  possibility  is  a  distribu- 
tion of  very  small,  hot  H  II  regions.   Altenhoff  et  al . 
(1973)  have  observed  10  discrete  continuum  sources  within  36 
arc  minutes  of  the  position  (1=  24°,  b  =  0°)  at  which 
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Gottesman  and  Gordon  (1970)  and  Gordon  and  Gottesman  (1971) 
observed  recombination  lines.   This  is  a  region  in  which 
Altenhoff  et  al .  (1970)  indicate  there  are  no  resolved  con- 
tinuum sources.   The  10  discrete  sources  were  observed  at 
a  wavelength  of  11  cm  and  with  a  spatial  resolution  of  4.4 
arc  minutes.   According  to  Guelin  (1974),  three  or  more  of 
these  sources  produce  25  km/s-wide  recombination  lines.   These 
are  about  as  wide  as  those  observed  in  the  direction  of  3C391. 
Weak  and  broad  (Av  =  80  km/s)  lines  were  observed  between 
these  sources. 

Assuming  that  these  "Hot  Spot"  sources  are  widely 
distributed  throughout  the  galactic  ridge,  expected  line 
powers  will  be  calculated;  the  calculations  will  be  based 
on  as  much  evidence  as  is  available.   The  expected  line 
powers  of  H157a  and  H92a  emission  will  be  compared  to  those 
observed . 

Guelin  (1974)  provides  no  more  information  on  these 
sources  than  has  been  mentioned  except  that  it  was  implied 
that  they  are  separated  by  a  minimum  of  4.5  arc  minutes. 
With  this  spacing,  13  of  them  can  be  distributed  within  the 
18  arc  minute  H157a  antenna  beam.   The  maps  of  the  region 
surrounding  3C391  made  by  Caswell  et  al .  (1971),  Dulk  and 
Slee  (1972),  Bridle  and  Kesteven  (1971),  and  Chaisson 
(1974a)  show  no  collection  of  discrete  sources  within  18 
arc  minutes  of  3C391.   Since  Guelin  does  not  mention  the 
intensity  of  the  sources,  one  cannot  be  certain  of  their 
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absence  since,  in  all  cases,  the  maps  were  presented  with 
the  background  removed.   The  structure  seen  by  Chaisson 
(1974a)  with  a  2.1  arc  minute  beam  at  15.5  GHz  is  enclosed 
by  the  weakest  isophotes  of  all  other  maps  and  is  probably, 
therefore,  structure  within  the  SNR  rather  than  of  the 
intervening  medium. 

Braes  and  Miley  (1972)  have  produced  an  aperture 
synthesis  map  of  a  l°xl°  region  in  which  Cygnus  X-3  is 
located.   This  is  in  the  Cygnus  X  region,  a  distribution  of 
H  II  regions  appearing  roughly  circular  in  shape  with  a 
diameter  of  about  six  degrees.   There  are  about  100 
known  thermal  continuum  sources  in  this  region  (Wendker, 
1970).   Radio  recombination  lines  have  been  observed  in 
several  places  between  these  continuum  sources  (Gottesman 
and  Gordon,  1970;  Seacord  and  Gottesman,  unpublished  obser- 
vations).  The  map  of  Braes  and  Miley  (1972)  does  show  some 
structure.   There  are  a  few  extended  sources,  one  of  which 
coincides  with  an  OB  association  (Cygnus  0B2  or  VI  Cygni) 
and  the  others  coincide  with  known  H  II  regions.  There  are 
about  15  discrete  sources  having  a  flux  density  considerably 
less  than  that  of  Cyg  X-3,  0.021  Jy.   However,  no  more  than 
five  of  these  sources  exist  in  any  18x18  arc  minute  area 
over  the  entire  l°xl°  map  and  some  of  these  sources  are 
closer  together  than  4.5  arc  minutes. 

There  are  two  cases  for  which  discrete  H  II  regions 
have  been  studied.   Habing  et  aj_.  (1973)  have  observed  H  II 
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regions  which  are  related  to  Type  I  OH  maser  sources. 

One  of  these,  ON  - 1 ,  had  previously  been  studied  by  Winnberg 

et  al_.  (1973).   ON- 1  is  located  at  1  =  69?5   and  b  =  -1?0 

o         _  c 

and  has  an  emission  measure  of  about  2.7x10   pc  cm   , 

5    -3 
an  electron  density  of  about  10   cm   ,  and  a  diameter  of 

0.03  pc  or  less.   Winnberg  ejt  aj_.  scanned  a  field  of  24x24 
arc  minutes  around  ON  - 1  with  a  resolution  of  1.3  arc  minutes 
and  found  no  other  sources.   If  the  sources  observed  by 
Altenhoff  et  aj_.   (1973)  are  like  that  associated  with  ON- 1  , 
it  is  apparent  that  they  do  not  occur  frequently  out  of  the 
galactic  ridge.   This  agrees  with  the  finding  that  radio 
recombination  lines  are  not  observed  in  the  general  inter- 
stellar medium  outside  the  ridge. 

Another  case  on  which  to  base  hot  spot  models  is 
the  hot,  dense  condensations,  or  "knots,"  within  emission 
nebulae.   An  example  of  this  kind  of  H  II  region  is  in  the 
core  of  M42.   For  the  knots  there,  Hjellming  and  Davies 

(1970)  find  that  the  average  temperature,  electron  density, 

4  4-3        7 

and  emission  measure  are  10   °K,   2x10   cm   ,  and  10   pc 

cm"   respectively.   These  averages  are  made  over  an  effective 

path  length  of  0.02  pc  which  is  assumed  to  be  the  size  of 

the  sources . 

A    series    ofmodels  for  a  distribution  of  0.02  pc  H   II   regions 

will    be  considered.     They  are  assumed  to  be  distributed  across  and  along 

an  18  arc    minute    beam    such    that    13    of    them   are    observed    in 


160 


it  and  there  is  one  of  them  in  the  2.3  arc  minute  H9  2a  beam. 
The  distribution  has  a  centroid  located  at  a  distance  of 
8.5  kpc  from  the  sun.   The  angle  subtended  by  one  of  these 
sources  is  about  0.48  arc  seconds.   Therefore,  the  beam 
dilution  factor,  pB,  is  1.65x10"   for  the  H157a  beam  and 
6.80xl0"6  for  the  H92a  beam. 

The  integrated  line  power,  P,  is  computed  with 
Equation  V-10.   LTE  and  non-LTE  cases  are  computed  for  which 
N  =  1  ,  10,  102,  and  10   cm"  .   The  temperatures  and  range 
of  emission  measures  are  based  on  the  results  of  Hejllming 
and  Davies  (1970)  and  Habing  et  al_.  (1974);  Tg  =  7.5xl03 
and  104  °K  with  104  <  EL  <  5x10   pc  cm"  .   The  non-LTE 
departure  parameters  were  computed  by  Sejnowski  and  Hjellming 
(1969)  and  are  tabulated  in  Table  V-5.   The  slope  of  the  bn- 

curve,  6mn  =  Cvbn)/bm'  is  alS°  shown- 

Figures  V-12  and  13  show  P(H157a)  and  P(H92a) 
plotted  as  a  function  of  emission  measure.   Figure  V-12 
is  for  Te  -  7500  °K  and  Figure  V-13  is  for  Tg  =  10000  °K. 
In  both  figures  the  H157a  curves  are  plotted  above  those  of 
H92a  and  the  vertical  line  labeled  with  PQ  represents  the 
observed  line  power:   PQ(H157a)  =  9.61  (+0.26,  -1.27)  °K  kHz 
and  P  (H92a)  =  5.7  (±2.4)  °K  kHz.   The  shaded  strip  repre- 
sents the  error  bars.   The  non-LTE  curves  are  labeled  with 
the  electron  density  for  which  the  departure  parameters  were 
computed.   The  H92a  LTE  curves  are  drawn  dashed.   For  both 
temperatures,  the  H157a  LTE  line  powers  are  less  than 
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Table  V-5 


Departure  Parameter  Used  in  the  Hot  Spot  Models 
[Computed  by  Sejnowski  and  Hjellming  (1969)] 


N  [cm-3] 
eL    J 

b 
m 

6 
mn 

mn 

H157a 

Te 

■  104 

°K 

10 

0.9024 

2. 71 OxlO-3 

-335.38 

102 

0.9739 

8.900xl0"4 

-109.27 

3 
10J 

0.9951 

1 .085xl0"4 

-  12.433 

io4 

1 .0000 

3. 550x1 0"5 

-   3.3947 

H157a 

Te 

=  7.5x 

IO3  °K 

10 

0.8858 

3.200xl0"3 

-297.05 

io2 

0.9701 

1 .lOOxlO"3 

-101 .24 

io3 

0.9922 

2.230xl0"4 

-  19.708 

io4 

1 .0000 

4.31xlO"5 

-   3.057 

H92a 

Tp 

=  IO4 

10 

0.7340 

1  .  428x1 0"3 

-  34.858 

,o2 

0.8056 

5. 750x1 0"3 

-144.02 

io3 

0.9265 

3.740xlO"3 

-  93. 133 

io4 

0.9810 

1 .OOOxlO-3 

-  24. 100 

H92a 

Tp 

-  7.5x 

IO3 

10 

0.6952 

1 . 670x1 0~3 

-  30.459 

!02 

0.7743 

6.850xl0"3 

-128.71 

,o3 

0.9138 

4. 1 70x1 0~3 

-  77.749 

!04 

0.9762 

1 . 200x1 0"3 

-  21.594 

Figure  V-12.   Integrated  Line  Power  as  a  Function  of 
Emission  Measure  for  a  7500°K  Hot  Spot. 
The  upper  curves  are  for  the  H 1  5  7  a  line. 
The  lower  curves  are  for  the  H92a  line. 
Each  non-LTE  curve  is  labeled  with  the 
electron  density  for  which  it  is  calcu- 
lated.  The  observed  line  power  is  repre- 
sented by  the  line  labeled  with  P0  and 
the  shading  represents  the  error  bars. 
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Figure  V-13.   Integrated  Line  Power  as  a  Function  of 

Emission  Measure  for  a  10  000°K  Hot  Spot. 
The  upper  curves  are  for  the  H157a  line. 
The  lower  curves  ar  for  the  H  9  2  a  line. 
Each  non-LTE  curve  is  labeled  with  the 
electron  density  for  which  it  is  calcu- 
lated.  The  observed  line  power  is  repre- 
sented by  the  line  labeled  with  P0  and 
the  shading  represents  the  error  bars. 
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0.01  °K  kHz  and,  therefore,  are  not  shown  in  the 
figures. 

Figure  V-12  indicates  a  consistent  solution  for 

Te  -  7500  °K  because  the  PQ(H157a)  and  PQ(H92a)  lines 

-3  ^ 

intersect  the  respective  Ng  =  100  cm   curve  at  EL  ^ 

2.6x10   pc  cm"  .   The  intersection  is  well  within  the  error 


limits 


The  situation  for  T   =  10000  °K  is  not  as  simple. 


The  solution  occurs  between  the  curves  for  Ne  =  10  and  100 
cm   .   Simple  interpolation  is  impossible  because  the  de- 
pendence of  the  departure  parameters  on  the  electron  density 

inverts  the  E.-P  relationship  of  the  HI 57a  and  H92a  curves 

_3 
between  N   =  10  and  100  cm   .   Consequently,  in  order  to 

obtain  a  consistent  solution,  the  functions  were  replotted 
as  P  being  a  function  of  N  ,  using  EL  as  a  parameter.   This 
is  shown  in  Figure  V-14.   The  locus  of  the  points  of  inter- 
section between  the  P   lines  and  the  constant  E.  curves 

is  shown  in  Figure  V-15.   The  H157a  and  H92a  loci  intersect 

7  -3 

at  the  values  EL  =  4x10   and  <KQ>      =  81  cm   .   The  inter- 
section occurs  within  the  range  dictated  by  Figure  V-13. 

Table  V-6  compares  the  derived  models  with  the  proto- 
types on  which  they  were  based.   The  agreement  between  the 
emission  measures  is  not  good  and  there  is  no  agreement 
between  the  electron  densities.   <N  >,  the  electron  density 
averaged  over  the  path  length  L,  is  at  least  two  orders  of 


Figure  V-14.   Integrated  Li ne  non-LTE  Power  as  a  Function 
of  Electron  Density  for  a  10  000°K  Hot  Spot 
The  upper  curves  are  for  the  H157a  line. 
The  lower  curves  are  for  the  H92a  line. 
Each  curve  is  labeled  with  the  emission 
measure  for  which  it  was  calculated.   The 
observed  line  power  is  represented  by  the 
line  labeled  with  a  P0  and  the  shading 
represents  the  error  bars. 
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pc  cm 


Fi  gure  V-l 5 . 


Hot  Spot  Model  Solution  for  T?  =  10  000°K 
The  solution  is  given  by  the  inter- 
section of  the  two  curves;  E[_  =  4  x  10' 
pc  cm"6  and  <  N  e  >  ' '   =81  cm" 


31 
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Table  V-6 
Comparison  of  Models  with  Their  Prototypes 


Emission  Nebula  Knots 


This  Study 


<Ne> 

<nV/2 

e 

2  1/2 
<Ne>/<Ng>'^ 


TT 


10000  °K 
4x10   pc  cm 
81  cm"3 
4.47xl04  cm"3 


0.16 
4.8 


M42 


(Hjellming  and  Davies,  1970) 


10000  °k 

1 . 26x1 0   pc  cm 

2.5xl04  cm"3 

3    -3 
6x10   cm 


M.2 


>1 


1 


H  II  Region  Associated  with  Type  I  OH  Masers 
This  Study 


<Ne> 

<  N  >/<N2>1//2 

Me  '   e 


7500  °K 

2 

6xl07  p 

c 

cm 

100  cm"3 

3 

64xl04 

cm 

-3 

0 

28 

5 

5 

0N-1 

(Winnberg  et  al . ,  1973 


-6 


8000  °K 
o 

2.7x10   pc  cm 


1 . 3x1 0   cm 


>1 


Notes 
*1 


Gordon  (1974 


'A  diameter  of  3660  au  is  assumed 
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magnitude  too  small  for  either  type  of  region.   <N  >  can  be 

compared  with  the  RMS  electron  density  computed  by  dividing 

the  emission  measure  by  the  path  length,  0.02  pc. 

.2  1/2 


Table  V-6  also  shows  the  ratio  <N  >/<N  > 


and  the 
1.683  GHz  continuum  optical  depth.   The  ratio  would  indicate 

the  amount  of  clumping  along  the  line  of  sight  if  <Nfi>  and 

2  1/2 
<N  >    were  obtained  by  independent  measurements.   Winnberg 

et  al  ■  (1973)  made  no  determination  of  <N  >.   Hjellming  and 
Davies  (1970)  solved  Equation  V-10  for  several  a-,  B-,  y-, 
6  -  ,  and  e-lines  to  obtain  values  of  <  T  e  >  >  <^e>'  a  n  ^  *" ^ I  >  ' 
So,  <N  >  and  the  RMS  value  of  it  should  be  independent  except 
for  the  fact  that  they  both  were  obtained  from  radio  re- 
combination line  observations.   This  is  true  in  the  case  of 
this  study  also.   An  independent  determination  of  <N  > 
would  be  through  the  measurement  of  a  pulsar's  dispersion 
measure  (cf.  Equation  11-55).   Since  it  is  probable  that  there 
are  more  electrons  along  the  line  of  sight  than  are  respon- 
sible for  the  radio  recombination  line  emission,  the  <N  > 

does  not  represent  an  average  over  the  line  of  sight  and, 

2  1/2 
therefore,  <N  >/<N  >    does  not  measure  clumping. 

However,  it  is  apparent  that  the  ratio  should  repre- 
sent the  same  quantity  in  both  the  analysis  in  this  study 
and  in  that  of  Hjellming  and  Davies  (1970)  and  would, 
therefore,  indicate  agreement  or  disagreement  with  the  type 
of  source  being  considered.   There  is  clearly  a  lack  of 
agreement,  not  only  with  this  ratio,  but  with  all  parameters 
der i  ved . 
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The  number  of  sources  in  the  H157a  beam  was  chosen 
to  be  the  maximum  number  allowed  with  a  minimum  intersource 
spacing  of  4.5  arc  minutes.   To  determine  what  effect  a 
smaller  number  of  sources  in  the  18  arc  minute  beam  would 
have  on  the  solution,  Figure  V-16  shows  the  line  power- 
emission  measure  curves  for  13,  10,  and  5  sources  in  the 
beam.   The  number  of  sources  enters  into  the  calculation  by 
the  beam  dilution  factor,  pg.   Each  curve  in  Figure  V-16 
is  labeled  by  the  number  of  sources.   The  upper  set  of  curves 

is  for  T   =  7500  °K;  the  lower  curve  is  for  Te  =  10000  °K. 

-3 
In  both  sets,  the  dashed  curves  are  for  Ne  =  100  cm   and 

-3 
the  solid  curves  are  for  Ne  =  1 0  cm 

Reducing  the  number  of  sources  has  the  effect  of 

reducing  the  electron  density  and  increasing  the  emission 

measure.  So,  the  size  of  the  emitting  region  increases 

greatly.   The  recombination  line  emission  seems  to  favor  an 

electron  density  more  like  that  of  the  more  tenuous  gas 

surrounding  the  dense  knots  of  an  emission  nebula.   This  is 

in  agreement  with  the  conclusion  of  Brockelhurst  and  Seaton 

(1972).   Consequently,  the  results  of  this  analysis  do  not 

seem  to  support  the  possibility  that  radio  recombination 

lines  observed  from  the  galactic  ridge  are  emitted  by  a 

collection  of  small  (^0.02  pc),  dense,  hot,  H  II  regions 

like  ON  - 1  or  the  condensations  in  emission  nebulae. 


Figure  V - 1 6 .   The  Effects  of  the  Number  of  Sources  on 
the  Hot  Spot  Model  Curves.   The  non-LTE, 
H157a  integrated  line  power  is  plotted  as 
a  function  of  emission  measure  for  Te  = 
7500°K  (upper  set  of  curves)  and  Te  = 
10000°K  (lower  set  of  curves).   The 
dashed  curves  in  both  sets  are  for  Ne  = 
100  cm-3;  the  solid  curves  are  for  Ne  = 
10  cm-3.   Each  curve  is  identified  by 
the  number  of  sources  in  the  antenna 
beam.   The  observed  line  power  is  iden- 
tified by  P0. 
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H  II  Nebula  Models 


Chaisson  (1974a)  suggests  that  the  H92a  line  ob- 
served by  the  Cesarskys  (1973b)  is  emitted  by  a  hot  (7500  °K), 

_3 
low  density  (N   =  22  cm   )  H  II  region  embedded  in  a  two- 
phase  medium  rather  than  by  the  cold  regions  which  the 
Cesarskys'  results  suggest.   There  are  several  arguments  in 
favor  of  this  which  also  apply  to  this  study's  spectra.   The 
line  widths  (^30  km/s)  are  like  those  observed  from  emission 
nebulae  which  are  much  larger  than  those  of  the  correspond- 
ing 21  cm  absorption  features  (£10  km/s).   Also,  the  80 
MHz  optical  depth  measured  by  Dulk  and  Slee  (1972)  can,  in 
part,  be  caused  by  a  partially  ionized  cool  region  with  the 
remaining  absorption  caused  by  an  H  II  nebula.   In  fact, 
recent  unpublished  observations  indicate  that  t(80)  may  be 
0.4  (Chaisson,  1974b)  which  is  the  value  it  would  be  if  the 
recombination  lines  were  emitted  only  by  the  nebula  (Chaisson, 
1974a).   He  detects  an  excess  continuum  flux  at  15.5  GHz 
which  could  result  f rom  t he  nebula.   However,  it  should  be 
noted  that  Goss  (1972)  points  out  that  his  31.4  GHz  observa- 
tions do  not  indicate  an  excess  flux.   Finally,  the  Cesarskys 
(1973b)  derive  an  ionization  rate  of  Cu^"10   ~  sec   from 

their  cold-region  model.   However,  Hjellming,  Gordon,  and 

-  1  5    - 1 
Gordon  (1969)  derive  Cu^lO    sec    from  pulsar  dispersion 

measurements  and  Hughes,  Thompson,  and  Colvin  (1971)  obtain 

a  range  of  1.9x10    <  u  <  4.0x10     sec    from  21  cm 

absorpti  on  spectra . 
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These  arguments  suggest  considering  an  H  II 
nebula  model.   The  model  curves  were  computed  with  Equation 
V-10  for  LTE  and  non-LTE  cases.   The  temperature,  T  ,  was 

restricted  to  7500  °K  and  the  non-LTE  models  were  computed 

-3 
for  electron  densities  of  10  and  100  cm   using  the  departure 

parameters  of  Sejnowski  and  Hjellming  (1969)  given  in  Table 
V-5.   Because  of  the  variation  of  the  line  intensities 
between  the  3C391  and  its  C0MP  beam,  it  was  assumed  that  the 
nebula  subtends  an  angle  of  about  34  arc  minutes,  giving  it 
a  diameter  of  85  pc  at  a  distance  of  8.5  kpc.   Therefore, 
since  it  completely  fills  the  H157a  and  H92a  beams,  the  beam 
dilution  factor  is  unity  for  both. 

Figure  V-17  shows  the  integrated  line  power  plotted 
as  a  function  of  emission  measure.   The  H157a  curves  are 
drawn  solid;  those  for  the  H92a  line  are  drawn  dashed.   Each 
set  consists  of  three  curves,  the  LTE  curve  and  the  non-LTE 
curves  labeled  by  the  electron  density  used  to  compute  it. 
The  observed  integrated  line  power  for  each  line  is  indi- 
cated by  the  horizontal  lines  labeled  with  P  (line).   As  for 
J  o 

the  previous  models,  PQ(H157a)  =  9.61  (+0.26,  -1.27)  °K  kHz 

and  PQ(H92a)=5.7  (±2.4)  °K  kHz. 

A  constant  solution  is  obtained  when  the  lines  of 

P„(H157a)  and  P  (H92a)  intersect  their  respective  LTE  or 
oo 

non-LTE  curve  at  a  common  value  of  Ei.   Because  of  the  large 
errors  involved,  it  is  difficult  to  ascertain  exactly  what 
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Figure  V - 1 7 


H  II  Nebula  Model  Solution.   The  inte- 
grated line  power  is  plotted  as  a 
function  of  the  line  emission  measure. 
The  H157a  curves  are  drawn  solid; 
those  for  the  H92a  line  are  drawn 
dashed.   The  non-LTE  curves  are  labeled 
with  the  electron  density  used  in 
their  computation.   The  observed  line 
powers  are  identified  by  P0  (line) 
and  the  error  bars  are  indicated  by 
the  shaded  strip. 
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values  apply  to  the  solution.   Perusal  of  the  areas  of  inter- 
section in  Figure  V-17  is  easier  with  a  table.   Thus,  Table 
V-7  indicates  the  value  of  emission  measure  where  the  lower 
error  limit,  P  (line),  and  upper  error  limit  intersect  each 
of  the  curves.   From  the  figure  and  table,  it  appears  likely 
that  the  solution  is  close  to  E^  =  3.8x10   pc  cm   with  Ne  = 
10  cm-3.   Since  PQ(H157a)  =  9.6  °K  kHz  is  closer  to  the 

upper  error  limit,  the  solution  E,  =  3.3x10   pc  cm   with 

_3 
Ne  =  100  cm   is  not  unlikely.   The  table  and  figure  indi- 
cate that  an  LTE  solution  is  unlikely  since  a  common  inter- 
section is  so  close  to  both  error  limits.   There  is  astro- 
physical  evidence  which  also  suggests thi s . 

The  likely  range  of  parameter  values  agrees  well  with 
those  Chaisson  (1974a)  obtains  by  an  independent  method: 
<Ne>  =  22  cm"   and  ET  ^4x10   pc  cm"  ,  giving  EL  =  Ey/1.1  % 

o         _  c 

3.6x10   pc  cm   .   This  solution  also  compares  favorably  with 
the  parameters  applicable  to  visible  H  II  nebulae.   For 
example,  the  North  American  nebula  complex  (W80)  has  a 

_  o 

diameter  of  about  50  pc  and  <N  >  %  10  cm    (Kaplan  and 
Pikelner,  1970,  p.  153)  which  gives  ET^5xlO   pc  cm 

Finally,  it  should  be  noted  that,  in  a  very    recent 
paper,  Jackson  and  Kerr  (1975)  conclude  that  radio  recombina- 
tion lines  originate  from  small,  low- density  H  II  nebulae 
like  this  study's  H  II  nebula  model  except  that  the  electron 
temperature  they  derive  is  close  to  4400  °K.   Their  results 
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will  be  discussed  in  more  detail  in  Chapter  VI  where  the 
ridge  as  an  entity  will  be  discussed.   Also,  as  this  was 
being  written,  Pankonin  (1975)  reported  on  his  H248arecombi - 
nation  line  observation  of  3C391.   His  line  indicates  nothing 
which  has  not  been  observed  before  and  he  favors  a  low 
density  H  II  region  as  the  source  of  the  line. 

Table  V-8  summarizes  the  models  derived  in  this 
chapter.   The  discussions  made  during  the  derivation  of  these 
models  suggest  that  the  most  likely  models  are  the  Isolated 
Cold  Region  and  the  H  II  Nebula  models.   In  the  following 
chapter,  additional  arguments  will  be  considered  from  which 
a  most  probable  model,  including  the  distribution  of  the 
emission  sources,  will  be  considered. 
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CHAPTER  VI 
FURTHER  CONSIDERATIONS  AND  APPLICATIONS  OF  THE  MODELS 

Discussion  of  the  Models 

This  discussion  pertaining  to  the  models  suggests 
that  the  most  plausible  condition  is  that  the  recombination 
line-emitting  region  is  either  an  isolated  cold  (Te~53  °K) 
region  or  a  low  density  emission  nebula  similar  to  W80. 
In  order  to  decide  on  a  final  model,  several  more  parame- 
ters will  be  calculated  and  compared  to  the  results  of  oth- 
er work.   Three  of  these  parameters  are  the  1683  and  80 
MHz  optical  depths  and  the  continum  antenna  temperature  at 
1683  MHz, which  are  expected  to  be  observed  from  each  of  the 
model s . 

The  D-Method  models  i n vol ve  a  two-component  medium 
which  will  be  considered  to  be  a  cold  region  surrounded  by 
hot  gas.   The  continuumantenna  temperature,  T-j- ,  is  given  by 


TT  =  nBUe(h)[l-exp(-Tf(h))]  +  ^  Te(c)exp(-Tf(h))[l 


exp(-i(c))]  +  T  (h)exp(-T(c)-Tf(h))[l-exp(xb(h)]} 
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In  this  equation,  ft   is  the  solid  angle  subtended  by  the 
cold  region.   For  these  models  it  is  greater  than  the  an- 
tenna main  beam  solid  angle,  ft..;  so,  ft/^M  =  1-   All  op- 
tical depths  are  for  the  thermal  continuum.   That  of  the 
hot  gas  in  front  of  the  cold  region  is  x .  ( h )  and  that  of 
the  hot  gas  behind  the  cold  region  is  x.(h).   It  will  be 
assumed  that  the  hot  gas  is  distributed  equally  in  front 
of  and  behind  the  cold  region;  so,  if(h)  =  Tb(h).   For  a 
single  component  recombination  line  emitter  with  a  contin- 
uum optical  depth  xT  and  observed  through  an  optically  thin 
gas  (t^3x10~  ),  Equation  VI-1  reduces  to 


nB  ^  Te[l-exP(-TT)], 


where  ft,  is  the  solid  angle  subtended  by  the  line-emitting 
region.   For  all  except  the  Hot  Spot  models,  ft,/ftw  =  1-   For 
the  Hot  Spot  models  Q,/ftM  =  1.65x10"  . 

The  continuum  optical  depths  at  1683  and  80  MHz 
and  the  1683  MHz  continuum  antenna  temperature  for  each  of 
the  models  are  listed  in  Table  VI-1.   If  the  isolated  cold 
region  were  surrounded  by  a  hot  gas,  the  situation  would 
be  similar  to  that  of  the  two-component  model.   From  the 
table  it  is  evident  that  this  hot  region  would  contribute 
a  neglibible  amount  to  the  antenna  temperature.   If  the 
Hot  Spots  are  surrounded  by  a  low  density  HII  nebula,  the 
observed  contribution  to  the  continuum  temperature  would 
be  about  4  °K. 
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Gordon  and  Gottesman  (1971)  observed  H1973  recom- 
bination lines  in  Cygnus  X  and  at  two  positions  in  the 
ridge,  at  longitudes  of  23?9  and  25?1.   The  ratio  of  H157a 
to  H 1 970  integrated  line  powers  are  3.21,  3.95,  and  4.56 
for  longitudes  23?9  and  25?1  and  Cygnus  X  respectively. 
The  H 1 9 7 3  line  powers  expected  in  the  direction  of  3C391 
are  calculated  for  each  model  using  the  parameters  of  Tab- 
le V-7  and  Equation  V-10  except  for  the  D-Method  models 
where  the  integrated  line  power  is  given  by 


P(H1 973 )  =  1 .0645Avd[0.8t£(H1976)  (T  -Tr)] 


The  power  ratio,  P(H1 57a)/P ( HI  976)  ,  for  each  model  is  listed 
in  the  fifth  column  of  Table  VI- 1.   The  velocity  half-pow- 
er width  of  the  line  used  to  calculate  line  power  was  taken 
to  be  same  as  that  observed  for  the  H157a  line,  35.7  km/s 
or  201.8  kHz.   A  comparison  of  the  power  ratios  in  Table 
VI-1  with  those  found  by  Gordon  and  Gottesman  is  made  with 
the  assumption  that  a  radical  change  from  region  to  region 
within  the  ridge  is  not  to  be  expected.   Acutally,  their 
ratios  did  not  change  much  between  the  ridge  regions  and 
that  in  Cygn  us  X . 

Another  parameter  which  will  be  used  to  evaluate 

the  models  is  the  hydrogen  ionization  rate,  cu  =  H  N,.  + 
J         a  H     e  H 

a?(T)/N„  (Equation  11-48).   The  partial  recombination  co- 
efficient, a?(T),  is  computed  from  Equation  11-22.   Then 
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a2(r 


2.065X10"11  T"1/2c2(T; 


and  ?p(T)  is  obtained  from  Figure  1 1  -  3 .   Because  radia- 
tive recombination  directly  into  the  ground  state  gener- 
ates a  photon  which  immediately  ionizes  a  neighboring  hy- 
drogen atom,  the  ionization  rate  from  outside  sources  is 
related  to  recombination  into  all  levels  expect  the  ground 
level.   Therefore,  the  partial  recombination  coefficient 
is  used  rather  than  the  total  coefficient  (Dalgarno  and 
McCray  1972). 

In  calculations  of  ^H  for  the  Hot  Spot  and  HII  Nebula 
models,  total  ionization  in  the  regions  is  assumed.   Then 
NH  =  NH  +  .   Also,  Ne  =  l.lNj,+  to  account  for  electrons  con- 
tributed by  helium.   Then  <N   Nu  +  >/NL,  reduces  to  <N  >. 

e   n     n  e 

For  the  D-Method  and  Isolated  Cold  Region  models, 
NeHH+  is  assumed  to  be  equal  to  N  .   The  total  hydrogen  den- 
sity is  nearly  equal  to  the  neutral  hydrogen  density,  Nu0 

H 

which  is  related  to  the  21  cm  line  integrated  optical  depth 
by  means  of  Equation  11-52; 

r        ,                              -u  r  ^u°(s  ) 
J  x21(v)dv  =  2.58x10  ID  /  -~ ds. 

line  o     s 

y\2-|(v)  dv  is  obtained  from  Table  2  of  Radha  kri  s  hnan  et  al. 
(1972b).   There  are  three  hi gh- vel oc i ty  features  in  the  ab- 
sorption spectrum  which  are  located  at  80,  101,  and  107  km/s. 
In  calculating  NH  for  the  cool  models,  is  assumed  that 
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N  ^  N  o  and,  since  the  effective  beam  width  of  the  21  cm 
H   H 

measurements  is  no  greater  than  the  18  arc  minute  beam  used 
in  this  study,  all  three  features  are  assumed  to  be  repre- 
sented by  a  blend  in  the  H157a  spectrum.   NH  is,  therefore, 
computed  using  all  of  them.   Then 


<  j  > 
lH>=3.876xl014— j-—  h2]  {' 


<T  > 
3.876xl014-T— [1  .0045  £t01  (i)AvnJ 


L  L'  -uu^  .21  v  '  /uvDi 
l 


18.97— j-^-  . 


The  path  length,  L,  is  in  pc.  <NH>  is  listed  in  the  sixth 
column  of  Table  VI-1.   The  fractional  ionization,  xeN6/Nh, 
is  listed  in  the  seventh  column. 

A  perusal  of  this  table  indicates  problems  with 
most  of  the  models.   For  all  cases  except  the  H  II  Nebula, 
the  recombination  line  emission  region  contributes  a  neg- 
ligible fraction  to  the  observed  continuum  temperature. 
In  these  cases  the  thermal  portion  of  the  total  continuum 
is  the  2.7  °K  universal  background  radiation. 

The    80    MHz    continuum    optical    depth    computed    on    the 
basis    of   model    parameters    shows    considerable    disagreement 
with   measured    values. The  value  oft  (80)  for  D-l    is  equal   to  that 
obtained    by      Caswell    et   al_.    (1970).    That  for  the  Isolated  Cold 
Region  model    is  a  bit  lower  than  the  value  of  0.83  derived  by  Dulk 


and  Slee  (1972).   t(80)  for  the  H  II  Nebula  model  is  equal 
to  that  which  Chaisson  (1974a)  suggests. 

Of  the  eight  models  considered  here,  only  three 
produce  H157a/H197B  ratios  which  are  reasonably  close  to 
those  obtained  by  Gordon  and  Gottesman  (1971),  model  D-l 
and  the  Isolated  Cold  Region  and  H  II  Nebula  models.   The 
extremely  high  ratio  for  model  D-4  results  from  P ( H 1 9  7  3)  a 
t*-El  -  1 .07xl0"3  pc  cm"6. 

If  the  recombination  lines  arise  from  a  cold  region, 

one  would  expect  to  find  a  hydrogen  density  between  10  and 

-  3  - 1 5     - 1 

100  cm   and  an  ionization  rate  near  10    sec   .   These 

values  are  in  agreement  with  the  two- component  model  of 
Field,  Goldsmith,  and  Habing  (1969)  for  which  the  ioniza- 
tion rate  is  4x10     sec    as  expected  from  2  Mev  cosmic 

ray  protons.   The  ionization  fraction  of  their  model  is 

-4 
about  5.5x10   .   None  of  the  D-models  agree  very   well  with 

their  model.   For  a  temperature  near  14.3  °K,  their  model 

-  3 
would  require  that  N H  ^  100  cm   .   Lock man  and  Gordon  (1973' 

obtain  a  low  temperature  ( T< 1 00  ° K )  solution  from  H159a  re- 
combination lines  observed  in  the  galactic  plane  in  regions 
near  the  center.   They  did  not  detect  the  H2003  line.  So, 
P (HI  59a) /P (H200B )  was  greater  than  would  be  expected  from 
a  gas  in  LTE.   This  result,  the  results  of  Cesarsky  and 
Cesarsky  (1973b),  Table  VI-1,  and  the  discussion  immediate- 
ly following  the  derivation  of  the  D-Method  models  suggest 
that  a  14.3  °K,  LTE  solution  is  not  valid. 
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The  Isol  ated Cold  Region  model  fares  only  a  little 
better  with  respect  to  NM,  x,  and  £„.   The  Field,  Goldsmith, 


and  Habing  (1969)  model  suggests  that  NH~30  cm    and  x 


,-3 


1.3x10    if  T   =  53  °K  and  the  carbon  abundance  is  a  f a c - 
e 

tor  of  ten  less  than  what  is  usually  assumed  (Spitzer,  1968, 
p.  122).   The  depletion  of  carbon  (due  to  accretion  on  neg- 
atively-charged grains)  is  supported  by  preliminary  analy- 
sis of  ultraviolet  observations wi th  satellites  (Dupree,  1974) 
and,  perhaps,  by  the  fact  that  the  C157a  line  was  not  ob- 
served in  this  study.   Nevertheless  the  number  densities 
and  fractional  ionization  of  the  model  derived  here  are    too 
small.   The  agreement  is  not  improved  by  introducing  a  time- 
dependent  source  of  ionization  as  suggested  by  Gerola  et  al  . 
(1974)  . 

The  hydrogen  densities  derived  for  the  Hot  Spotmodels 

are  much  too  small  since  the  prototypes  require  that  N   > 

3    -  3 
10   cm   .   It  is  difficult  to  imagine  how  such  an  electron 

density  can  be  obtained  in  the  Hot  Spot  models  where  N,,  ^ 

_  3 
100  cm   .   Evaluation  of  these  models  by  means  of  the  ioni- 
zation rates  is  not  possible  since  the  source  of  energy  is 
not  known. 

Of  all  models  studied,  the  H  II  nebula  produces  the 
most  promi  s  i  ng  set  of  parameters.   The  value  of  N,,  derived 
is  reasonable  when  compared  to  that  of  W80  (Kaplan  and 
Pikelner  1970,  p.  153).   If  the  region  is  ionized  by  a  hot 
star,  total  ionization  is  to  be  expected.   The  ionization 
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rate  due  to  a  star  with  a  flux  S   at  frequencies  greater 

u  3 

than  that  of  the  Lyman  limit  is  (Equation  11-50) 


3S 


'H    4irr3  Nu 
s   H 


S   is  tabulated  by  Spitzer  (1968,  p. 116).   The  Stromgren 

radius,  r  ,  for  this  model  is  42.5  pc.   Then,  for  an  05 

-  1  3     -  1 
star,  Cm  =  3.61x10    sec   .   If  r   is  decreased  to  25  pc, 
H  s  p 

-  1  2     -  1 
the  size  of  W8  0,  cH  =  1.77x10    sec   .   In  either  case,  a 

cooler  star  provides  a  smaller  ionization  rate.   In  order 

to  obtain  the  variation  of  line  strength  between  the  3C391 

and  C0MP  beams,  the  size  of  the  nebula  cannot  be  very    much 

less  than  85  pc.   So,  in  order  to  exist  with  an  ionization 

-  1  2     -  1 
rate  of  3.27x10    sec   or  less,  the  nebula  must  be  a  small 

cluster  of  nebulas  together  occupying  a  region  85  pc  in  dia- 
meter. 

The  arguments  favoring  and  rejecting  the  models 
studied  here  are  summarized  in  Table  VI- 2.   It  is  apparent 
that  the  best  model  for  the  source  of  radio  recombination 
lines  from  the  galactic  ridge,  at  least  in  the  direction  of 
3C391 ,  is  a  low-density  H  II  region.   Other  aspects  of  this 
model  —  particularly  the  relation  of  H  I  absorption  spectra 
to  it — will  now  be  discussed. 

A  recent  paper  by  Jackson  and  Kerr  (1975)  supports 
this  conclusion  except  that  for  their  model  they  derive  a 
temperature  of  440  (+600)  °K.   Also,  they  claim  that 
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consideration  of  non-LTE  will  not  affect  their  results  sig- 
nificantly.  Their  argument  for  assuming  LTE  was  discussed 
in  an  earlier  paper  (Jackson  and  Kerr,  1971)  where  they  show 
that  the  ratio  of  the  line  temperature,  T.  ,  expected  from 
a  non-LTE  gas  to  the  line  temperature,  Tj* ,  expected  from 
an  LTE  gas  is  roughly  equal  to  b  .   They  assume  that  b  =1 
for  the  HllOa  line  (v  =  4874  MHz).   However,  consulting  Fig- 
ure II-4,  it  is  seen  that  b  -,  -,  -,  <  0.7.   The  best  way  of  de- 
termining how  much  the  kinetic  temperature  is  above  4400  °K, 
is  to  observe  higher  order  lines  with  rest  frequencies  close 

to  that  of  the  HllOa  line.   The  line  strength  depends  on  3 

3     v  mn 

as  well  as  b  .   The  gas  in  these  nebulae  may  not  be  far  out 

m 

of  LTE  because  the  solution  shown  in  Figure  V-7  indicates 
only  a  small  difference  between  the  non-LTE  solutions  for 

_  o  o  _  c 

N   =  10  and  100  cm    if  E ,  ~ 3 x 1 0   pc  cm   .   In  this  case  there 
e  L       r 

_  3 
is  no  appreciable  difference  between  the  100  cm   non-LTE 

and  the  LTE  solutions. 

Hjellming  and  Davies  (1970)  determine  that  the  re- 
combination line  emission  regions  in  nebulae  such  as  M42 
and  Ml  7  are  in  non-LTE  even  though  they  are  denser  and  hot- 
ter than  is  the  model  derived  by  Jackson  and  Kerr  (1975). 
It  has  been  pointed  out,  however,  that  the  recombination 
lines  arise  from  the  more  tenuous  gas  which  surrounds  the 
dense  knots.   It  is  difficult  to  explain  why  such  knots 
would  be  out  of  LTE  since,  with  their  temperature  and  elec- 
tron density,  collisional  excitation  must  dominate  for 
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large  n.   So,  the  recombination  lines  must  come  from  the 
low  density,  non-LTE  gas  surrounding  them. 

It  is  interesing  to  note  that  the  model  developed 
in  this  study  is  not  very    much  different  than  the  model  de- 
veloped by  Matthews,  Pedlar,  and  Davies  (1973).   The  aver- 
age temperature  and  emission  measure  which  they  compute  from 
H166a  lines  are  6.2xl03  °K  and  2.7xl03  pc  cm"6  respective- 
ly  as  compared  to  7.5x10   °K  and  3.8x10   pc  cm    derived  by 
this  study.   Like  Jackson  and  Kerr,  they  assume  a  condition 
of  LTE.   Another  difference  between  their  model  and  that 
of  this  study  is  the  size  of  the  region.   The  positions 
which  they  observed  were  spaced  5°  to  10°  apart  and   they 
could  not  observe  the  change  in  the  spectrum  over  a  dis- 
tance as  small  as  that  observed  here.   From  their  observa- 
tions they  conclude  that  the  emission  regions  are  extended 
3°  to  5°  in  the  galactic  plane  which  is  equivalent  to  a 
size  in  order  to  0.5  to  10  kpc.   This  study  suggests  that 
the  structure  is  much  smaller,  being  in  the  order  of  85  pc 
or  smaller.   This  is  in  agreement  with  the  conclusions  of  Jack- 
son and  Kerr  (1975). 

It  has  been  mentioned  that  high  spatial  resolution 
maps  of  the  region  surrounding  3C391  (cf.   Chaisson  1974a) 
do  not  indicate  the  presence  of  discrete  sources  like  the 
ten  such  sources  reportedly  observed  by  Altenhoff  et  al . 
(1973)  near  the  direction  1  =  23?92,  b  =  -0?32.   But,  the 
possibility  of  their  existence  in  the  direction  of  3C391 
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cannot  be  discounted  completely.   If  they  do  exist,  it  is 
not  completely  clear  how  they  fit  into  the  model  proposed 
by  this  study.   Without  knowing  more  about  the  size,  dis- 
tribution, continuum  flux  density,  and  recombination  line 
spectra  of  these  sources,  it  can  only  be  assumed  that  they 
form  a   small-scale  structure  within  the  H  II  nebulae. 
Perhaps  they  are  dense  condensations  associated  with  star 
formation.   If  that  is  so,  it  is  not  clear  why  they  are 
not  seen  in  Chaisson's  (1974a)  map.   The  structure  seen  in 
his  map  is  included  within  the  isophotes  of  the  map  of 
Bridle  and  Kesteven  (1971)  and  it  appears  to  be  struc- 
ture within  3C391  and  not  of  the  intervening  medium. 

There  must  be  an  explanation  for  the  alignment  of 
the  21  cm  absorption  features,  which  are  associated  with 
cold  gas,  with  the  H  II  nebulae.   This  can  be  explained  by 
cool  cloud-like  regions  which  lie  close  to  the  nebulae. 
Being  close  to  each  other,  they  will  both  have  the  same  cen- 
tral velocities.   This  arrangement  is  similar  to  the  struc- 
ture of  the  Orion  nebula  described  by  Balick,  Gammon,  and 
Hjellming  (1974). 

The  distribution  of  the  H  II  nebulae  must  explain 
the  observed  central  velocities  and  relative  intensities  of 
the  recombination  line  features  from  the  region  around  3C391 
including  the  COMP  beam.   If  circular  galactic  rotation  can 
be  assumed,  the  Schmidt  rotation  model  (Schmidt,  1965)  can 
position  the  emission  regions  along  the  line  of  sight.   In 
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the  region  inside  the  sun's  orbit,  there  is  a   twofold  am- 
biguity between  velocity  and  radial  distance.   FiguresVI-1 
and  2  show  the  positions  of  the  emission  regions.   Figure 
VI-1  is  for  the  case  where  the  high-velocity  feature  is 
single;  Figure  V I -  2  is  for  the  case  where  it  is  double.   In 
both  figures,  an  outer  boundary  with  a  radius  of  8  kpc  is 
drawn  in  to  encompass  the  low-velocity  feature. 

The  outer  boundary  of  the  ridge  has  been  generally 
thought  to  be  a  circular  disc  or  an  annulus  with  an  outer 
radius  of  7  kpc  (Matthews,  Pedlar,  and  Davies,  1973).   If  the 
ridge  is  a  disc-shaped  distribution  of  H  II  nebulae  with  a 
nebular  density  which  drops  off  rapidly  beyond  7  kpc,  the 
boundary  may  not  be  well  defined.   At  longitudes  progressive- 
ly near  the  tangent  to  the  boundary,  the  lines  of  sight  will 
pass  through  fewer  nebulae  and  recombination  line  emission 
will  become  patchy.   So,  even  though  recombination  lines 
have  not  been  detected  in  the  direction  of  the  tangent  to 
the  8  kpc  circle  (1  =  53°)(Gordon  and  Cato,  1972),  it  is 
reasonable  to  expect  that  some  nebulae  lie  this  far  from 
the  galactic  nucleus. 

The  Hughes,  Thompson,  and  Colvin  (1971)  21  cm  absorp- 
tion spectrum  shows  an  absorption  feature  at  42  km/s.   There 
is  a  weak  feature  at  this  velocity  in  the  Radhakri shnan  e t  a! , 
(1972b)  spectrum.   From  this  it  can  be  assumed  that  either 
3C391  is  beyond  14  kpc  or  the  H  II  nebula  and  its  neighbor- 
ing cold  cloud  responsible  for  the  42  km/s  features  are  lo- 
cated at  a  distance  of  about  3.1  kpc  from  the  sun. 
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If  these  H  II  nebulae  are  to  explain  the  spectra  of 
other  galactic  ridge  surveys,  it  must  be  hypothesized  that 
the  ridge  phenomena  are  caused  by  a  distribution  of  H  II 
nebulae.   It  is  confined  to  a  disk-  or  ring-shaped  region 
extending  about  70  pc  above  and  below  the  plane  (Gordon, 
Brown,  and  Gottesman,  1972)  and  limited  to  a  gal  actocentri c 
radius  of  about  8  kpc.   In  order  to  test  this  hypothesis, 
the  spectra  from  the  directions  of  the  pulsar,  3C396,  and 
its  COMP  beam  will  be  examined.   Following  this,  the  re- 
sults of  this  study  will  again  be  compared  with  other  sur- 
veys of  the  plane. 

The  Pulsar  1858+03 

Determining  the  distance  to  a  pulsar  is  not  a  triv- 
ial task.   The  distance  to  the  stronger  pulsars  can  some- 
times be  determined  from  the  21  cm  absorption  spectrum  which 
is  obtained  by  subtracting  the  pulse-off  spectrum  from  the 
pulse-on  spectrum  (Gordon  e_t  a_l_.  1969).   An  estimate  of  the 
distance  can  be  calculated  from  the  dispersion  measure, 


DMe  /  N   ds  =  <N  >  L, 
/   e        e 


if  the  mean  electron  density  can  be  determined  independent- 
ly.  The  determination  of  <N  >  is  particularly  difficult 
in  the  galactic  plane.   A  recent  survey  of  interstellar 
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medium  research  (Guelin,  1974)  indicates  that  the  average 

electron  density  within  about  2  kpc  of  the  sun  is  0.03 

_3 
cm   .   <N  >  is  much  larger  in  the  ridge. 

The  reason  for  choosing  the  pulsar  1858  +  03  as  a  di- 
rection in  which  to  observe  radio  recombination  lines  is 

that  it  appears  to  lie  in  the  ridge.   Its  dispersion  measure 

-  3 
is  402  (+2)  pc  cm   (Manchester  and  Taylor,  1972)  which  is, 

by  far,  the  largest  dispersion  measure  known  as  of  late 
1971  . 

The  H157a  spectrum  from  the  direction  of  1858+03  is 
shown  in  Figure  IV-5.   The  center- vel oci ty ,  peak  tempera- 
ture, half-power  width,  and  the  integrated  line  power  of 
each  feature  are  brought  from  Table  I V- 3  and  shown  in  Table 
VI-3.   In  addition  to  these  parameters,  the  Doppler  temper- 
ature of  each  feature,  computed  from  Equation  11-37,  is 
shown  in  the  last  column  of  this  table.   The  Doppler  temper- 
ature is  the  kinetic  temperature  of  the  gas  if  the  line 
broadening  were  entirely  thermal. 

Table  VI-3 
Spectral  Parameters  for  1858+03 


vc 
km/s 

T(v    ) 
c  ; 

°K 

v 
km/s 

P 
°K    kHz 

TD 
°K 

38.08 
82.48 

0.013 
0.010 

30.16 
18.71 

2.34 
1  .12 

19    890 
7    654 
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Figure  V I -  3  shows  the  possible  positions  along  the 
line  of  sight  of  the  other  two  features  as  determined  by 
the  Schmidt  rotation  model.   The  near  positions  are  indi- 
cated by  dots  and  labeled  A  and  B.   The  far  positions  are 
indicated  by  x's  and  labeled  A1  and  B'  respectively.   A 
boundary  with  an  8  kpc  radius  is  shown. 

The  technique  to  determine  the  distance  to  the  pul- 
sar will  be  to  compute  the  expected  dispersion  measure  at 

points  along  the  line  of  sight.   The  distance  at  which  the 

_  3 
computed  dispersion  measure  equals  402  pc  cm   will  be  the 

distance  to  the  pulsar.   The  line  of  sight  passes  through 
four  different  types  of  region,  each  with  its  own  electron 
density  and  path  length.   They  are  the  region  outside  the 
ridge,  cold  H  I  regions,  the  diffuse  hot  gas,  and  the  H  II 
nebula.   The  arrangement  with  one  H  II  nebula  is  shown  sche- 
matically in  Figure  VI-4.   The  parameters  of  the  regions 
used  to  computed  the  distance,  R,  to  the  Pulsar  1858+03  are 
listed  in  Table  VI-4. 

Up  to  the  ridge  boundary,  the  accumulated  disper- 

_3 
sion  measure  is  73  pc  cm   .   If  a  nebula  is  encountered,  the 

dispersion  measure  will  rise  sharply.   If  the  path  length 

through  the  nebula  is  32.4  pc,  the  dispersion  measure  on 

_  3 
the  far  side  of  it  will  be  402  pc  cm   .   The  diagram  in  Fig- 
ure VI-4  purposely  shows  the  line  of  sight  not  passing  through 
the  center  of  the  nebula  to  illustrate  that  the  path  length 
through  it  may  be  less  than  its  diameter.   If  the  line  of 
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Table    VI -4 
Parameters    of   the    Regions    Shown    in    Figure    IV-4 


Parameter 


Value 


Reference 


Path  Length,   Sun 
to  Ridge  Boundary 


<N  >  Outside  Ridge 
e  3 


2.6  kpc 


0.03  cm 


Figure  VI-4 
Guelin  (1974) 


Diameter  of  Cold 
H  I  Region 


5  pc 


Baker  and  Burton 
(1975) 


Ne  in  H  I  Regions 


0.055  cm" 


Field,  Goldsmith, 
and  Habing  (1969) 


Average  Distance 
Between  H  I  Region 


333  pc 


Radhakrishnan  and 
Goss  (1972) 


<Ne>  in  Surrounding 
Hot  Region  (b=0) 


0.17  cm" 


Baker  and  Burton 
(1975) 


Path  Length  through 
H  II  Nebula 


adjustable 


Figure  VI-4 


<N„>  in  H  II 
Nebula 


10  cm" 


This  work 
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sight  misses  the  nebula  entirely,  the  distance  to  the  pul- 
sar must  be  15.9  kpc,  which  is  beyond  A'.   This  distance  is 
based  on  a  path  length  through  a  two-phase  medium  for  which 
the  dispersion  measure  increases  30.17  pc  cm   for  each  kpc 
of  path  length.   The  distance  of  15.9  kpc  is  also  consider- 
ably larger  than  that  derived  for  other  pulsars,  particular- 
ly those  for  which  their  lines  of  sight  encounter  H  II  re- 
gions (Terzian,  1972;  Prentice  and  ter  Haar,  1969). 

Although  it  is  possible  that  the  line  of  sight  en- 
counters neither  nebula,  it  seems  more  likely  that  it  pass- 
es through  part  of  one  since,  otherwise,  the  pulsar  would 
lie  beyond  the  solar  orbit.   Hence,  a  distance  of  2.6  kpc 
to  1858+03  appears  most  reasonable. 


The  3C396  and  COMP  Spect r a 


These  spectra  will  be  used  as  a  final  check  on  the 
model  which  has  been  developed.   The  spectra  are  shown  in 
Figure  IV-3.   The  results  of  Gaussian  profile  fitting  for 
the  3C396  spec trum  are  shown  in  Figure  I V- 9  and  the  para- 
meters of  the  fit  are  given  in  Table  IV-4.   The  reliabil- 
ity of  what  may  be  features  at  velocities  of  58  and  -119 
km/s  is  relatively  poor  since  the  peak  temperatures  are  on- 
ly about  four  times  the  RMS  noise.   Nevertheless,  by  compar- 
ing the  3C396  with  its  COMP  spectrum  in  Figure  IV-3,  it  is 
seen  that  the  peak  temperature  of  the  54  km/s  feature  de- 
creases by  more  than  a  factor  of  two  in  the  distance  between 
the  two  antenna  beam  centers  which,  is  30  arc  minutes. 
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The  absorption  spectrum  in  the  direction  of  3C396 
observed  by  Hughes,  Thompson,  and  Colvin  (1971)  consists 
of  a  wide  envelope  extending  from  -10  to  93  km/s  and  appears 
to  be  a  blend  of  many  features.   There  is  no  structure  at 
53  or  58  km/s.   The  21  cm  emission  spectrum  is  flat  and  ex- 
tends over  the  same  range  as  does  the  absorption  spectrum. 

A  means  of  checking  the  model  with  the  3C396  spec- 
trum is  to  compute  the  optical  depth  at  80  MHz  and  compare 
it  with  that  measured  by  Dulk  and  Slee  (1972).   The  line 
emission  measure  is  computed  from  Equati ons  V-l 0  and  11-41 
graphically  by  using  Figure  V-17.   Using  the  graph  sets 

T   =  7500  °K.   If  a  non-LTE  gas  with  N   =  10cm"3  is  con- 
e  J         e 

o  _  r  _  o 

sidered,  E.  =  2.30x10   pc  cm"  .   A  gas  with  N   >  100  cm" 

o  _  c 

is  in  LTE;  then,  E.  =  2.15x10   pc  cm    .   It  will  be  assumed 

tha t  the  recombi nation  line  emission  region  is  similar  to 

that  in  the  direction  of  3C391  (ie.  non-LTE  with  N 

e 

-  3 
10  cm   ).   Therefore  it  will  be  assumed  that  ET  =  1.1E.  = 

3      -3 
2.53x10   pc  cm    in  the  direction  of  3C396.   By  the  use  of 

Equation  V-2,  t(80)  -  0.237.   Dulk  and  Slee  obtain  t(80)  = 
0.49  (+0.3).   The  values  differ  by  a  factor  of  two,  but  are 
in  agreement  within  the  errors  quoted.   The  di ssagreement 
is  about  the  same  as  it  is  between  Dulk  and  Slee's  value 
for  the  line  of  sight  to  3C391  and  that  derived  from  the  H  II 
Nebula  model  as  shown  in  Table  VI-1.   In  this  case,  how- 
ever, the  agreement  with  Chaisson  (1974a,  b)  is  good.   It  is 
possible  that  Dulk  and  Slee's  values  for  t(80)  may  be  con- 
sistently large  by  a  factor  of  two. 
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Distribution  of  the  Recombination  Line 
Sources  Within  the  Ridge 


The  major  objective  of  this  study  has  been  to  de- 
termine what  kind  of  region  is  responsible  for  radio  re- 
combination line  emission  observed  from  the  galactic  plane. 
The  best  model  is  a  distribution  of  small,  low-density  H  II 
nebulae.   The  last  property  of  the  model  to  be  considered 
in  this  study  is  the  distribution  of  these  nebulae  within 
the  ridge.   The  recombination  line  surveys  which  have  been 
mentioned  can  be  used  to  obtain  a  rough  picture  of  this  dis- 
tribution.  The  accuracy  of  this  picture  is  severely  limited 
by  the  resolution  of  the  spectra  —  which  prevents  the  separa- 
tion of  blended  line  features  —  and  non-circular  orbital  mo- 
tion (e.g.  streaming)  of  the  interstellar  gas.   Even  if 
circular  motion  is  assumed,  placement  of  the  nebulae  along 
their  lines  of  sight  is  confused  by  the  two-fold  ambiguity 
of  the  radial  velocity -distance  relation. 

In  spite  of  these  limitations,  a  "first  order"  ap- 
proximation to  the  distribution  can  be  drawn  from  the  radio 
recombination  line  observations  in  the  ridge  from  longi- 
tudes between  5°  and  130°.   The  surveys  not  included  in 
this  summary  are  those  from  regions  closer  than  5°  from  the 
galactic  center  made  by  Lockman  and  Gordon  (1973)  and  Matthews, 
Davies,  and  Pedlar  (1973).   They  are  excluded  because  there 
is  no  radial  velocity  from  which  a  position  can  be  obtained. 
All  Doppler  shifts  are  due  to  streaming  of  the  gas.   Both 
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surveys  show  considerable  variation  of  the  line  profile  over 
a  distance  of  a  few  tenths  of  a  degree. 

Table  VI-5  lists  the  observations  made  in  the  range 
5°  < 1 < 1 30  ° .   If  at  some  longitude  no  line  was  observed, 
"none"  is  entered  in  the  last  comumn.   For  several  longi- 
tudes, it  is  apparent  that  a  wide  feature  is  a  blend  of 
unresolved  features.   For  these  cases,  the  central  velocity 
of  the  end  features  has  been  estimated  by  using  the  turn- 
over points  of  the  blend.   At  the  bottom  of  the  table  is 
the  list  of  references— each  with  the  spectral  line  observed- 
indicated  in  the  third  column.   The  possible  locations, 
based  on  the  Schmidt  circular  rotation  model,  are  shown 
in  Figure  VI-5.   The  near  and  far  positions  are  indicated 
by  a  dot  and  an  x  respectively.   The  blended  features  are 
shown  as  bars  to  indicate  the  range  in  which  the  features 
are  located.   Since  Cygnus  X  has  been  observed  to  be  a  region 
of  extensive  radio  recombination  line  emission  (e.g.,  Gottes- 
man  and  Gordon,  1970;  Seacord  and  Gottesman,  unpublished  ob- 
servations) its  angular  extent  is  shown.   It  lies  along  the 
local  spiral  arm  or  spur. 

In  Figure  VI-5  the  ridge  boundaries,  defined  by  the 
extent  of  the  nebula  positions,  are  drawn.   It  is  seen  that 
the  radio  recombination  line  ridge  is  an  annulus  with  an  in- 
ner radius  of  4  kpc  and  an  outer  radius  of  8.7  kpc.   The 
distribution  shown  in  the  figure  is  quite  similar  to  that 
of  the  giant  H  II  region  as  shown  by  Mezger  (1970).   He 


210 


Table  V  I  -  5 
Radio  Recombination  Line  Sources  in  the  Ridge 


1 

b 

Reference 

Velocity    Features 
[km/s] 

5.0 

0 

MPD 

10 

7.5 

0 

MPD 

5  to   36 

9.4 

0.4 

GC 

12  to 

22 

13.5 

0 

MPD 

15  to 

40 

15.8 

0 

GC 

20  to 

60 

17.5 

0 

MPD 

18  and  41 

21.0 

0 

MPD 

None 

21.4 

0 

GC 

80 

23.92 

-0.32 

GG  and 

JK 

75 

25.07 

+0.01 

GG 

50  and   100 

25.07 

0.01 

GC  and 

JK 

70 

25.1 

0.0 

MPD 

50  anc 

i  no 

27.64 

0.14 

JK 

None 

29.18 

0 

JK 

58 

29.5 

0 

MPD 

95 

30.27 

0 

MK 

30,    100 

31.2 

0 

JK 

25,    100 

31.9 

0 

S,   CC, 

P 

41.9, 

80.7,    102.8 

32.1 

-0.2 

s 

42.2, 

83.4,    106.0 

32.3 

0 

GC  and 

JK 

None 

33 

0 

GBG 

100 

36.5 

0 

MPD 

None? 

37.2 

-0.6 

S 

38.0, 

82.5 

38.9 

0 

S 

58.2 

39.2 

-0.3 

S 

53.9 

40.0 

0 

MPD 

None 

40.5 

0 

GC 

None 
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Tabl e  VI- 5  -  conti  nued 


Reference 


Velocity  Features 
[km/s] 


44.2 
47.5 
52.2 
70.0 

90 

no 

130 


GC,  P 

MPD 

GC 

GC 

GC 

GC 

GC 


40  to  87 

None 

None 

-4 

None 

None 

None 


Key  to  References  and  Spectral  Line 

MPD  Matthews,  Pedlar,  and  Davies  (1973);  H166a 

GC  Gordon  and  Cato  (1972);  H157a 

JK  Jackson  and  Kerr  (1975);  HllOa 

GG  Gottesman  and  Gordon  (1970);  H157a 

S  This  study  ;  HI  57a 

GBG  Gordon,  Brown,  and  Gottesman  (1972);  H157a 

CC  Cesarsky  and  Cesarsky  (1973a);  H166a  and  H157a 
Cesarsky  and  Cesarsky  (1973b);  H92a 


Pankonin  (1975);  H248a 
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defines  the  "giant"  H  II  regions  as  those  which  satisfy 

the  criterion  S(5  GHz)-D2  >  400;  S(5  GHz)  is  the  6  GHz 

n  — 

continuum  flux  and  D   is  the  nearer  of  the  two  kinematic 
n 

distances. 

Relevant  to  the  distribution  of  the  H  II  nebulae 
is  their  relation  to  H  I  gas.   Comparing  the  21  cm  absorp- 
tion and  emission  spectra  and  the  H157a  spectra  of  3C391 
(Figure  V - 1 )  and  those  of  3C396  indicated  that  the  21  cm 
spectra  show  more  features  and  are  extended  over  a  greater 
velocity  range  than  is  so  for  the  recombination  line  spec- 
tra.  However,  it  should  be  noted  that  there  is  an  H  I  emis- 
sion peak  close  to  each  recombination  line  peak.   A  compari- 
son of  the  radial  velocities  of  H  I  and  H  II  in  the  galac- 
tic pi  ane  within  60°  of  either  side  of  the  galactic  center 
shows  that  it  is  generally  true  that  H  II  features  coincide 
with  H  I  emission  peaks  (Mezgere_taJ_.,  1970).   One  or  more 
H  I  regions  may  be  physically  related  to  an  H  II  nebula  in 
a  way  that  has  been  suggested  for  the  Orion  Nebula  (Balick, 
Gammon,  and  Hjellming,  1974). 

When  relating  radial  velocities  to  distance,  it 
must  be  noted  that  Burton  (1971)  has  shown  that  the  Schmidt 
model  does  not  adequately  explain  the  observed  H  I  velocity- 
longitude  maps.   Gas  streaming  causes  the  rate  of  change  of 
velocity  along  the  line  of  sight,  dv/ds,  to  change  across 
a  spiral  arm.   The  velocity -distance  relation  is,  therefore, 
very  complicated.   Consequently,  Figure  V I - 5  does  not  represent 
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the  actual  nebular  distribution  well  enough   to  relate  the 
nebulae  to  spiral  arm  structure.    The  only  conclusion 
that  can  be  drawn  from  the  figure  is  that,  except  for  Cygnus 
X,  these  nebulae  are    restricted  to  an  annulus  with  an  outer 
radius  of  about  8.7  kpc  and  an  inner  radius  of  4  kpc.   It 
is  possible  that  the  wide  spectrum  from  the  longitude  of 
44?2  may  be  a  blend  of  features  from  a  string  of  H  II  re- 
gions along  the  inner  edge  of  the  Sagittarius  spiral  arm. 
If  this  is  correct,  this  string  of  nebulae  and  the  Cygnus 
X  complex  do  indicate  a  possible  relationship  between  the 
radio  recombination  line  ridge  and  galactic  spiral  struc- 
ture . 


CHAPTER  VII 
SUGGESTIONS  FOR  FURTHER  RESEARCH 

The  results  of  this  study  do  not  conclusively  de- 
termine a  model  for  the  radio  recombination  emission  regions 
However  a  distribution  of  low  density  (Ne  =  10  cm   )  H  II 
nebulae  seems  to  fit  the  data  of  this  and  previous  studies 
most  consistently.   There  are  a  number  of  aspects  of  the 
problem  which  need  further  study.   These  will  be  considered 
in  this  final  chapter. 

Assuming  the  conclusion  that  recombination  line 
emission  arises  from  low  density  nebulae  is  correct,  the 
aspects  of  the  general  model  which  remain  in  question  are 
as  foil ows : 

1 .   Size  of  the  nebulae 

2  .   Radial  and  z-direction  distribution  of  the 
nebul ae 

3.  Thermodynamic  conditions  in  the  nebulae; 
LTE  vs.  non-LTE 

4.  Energy  source  within  the  nebulae 

5.  Relation  between  nebulae  and  low  frequency 
absorption 
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6.  Relation  of  nebulae  to  H  I  gas,  both  hot 
and  cold;  and  galactic  structure  in  general 

7.  General  relation  with  discrete  sources  seen 
by  Altenhoff  et  al_.  ,  1973 

The  size  of  the  nebulae  can  best  be  found  by  high 
resolution  radio  recombination  line  observations.   The 
Altenhoff  e_t  aj_.  (1970)  continuum  maps  show  considerable 
structure  in  the  plane.   To  avoid  confusion  with  recombina- 
tion line  emission  from  large  H  II  complexes  like  W51,  W49, 
or  M42,  observations  should  be  confined  to  the  regions  like 
those  observed  by  Matthews,  Pedlar,  and  Davies  (1973), 
Gottesman  and  Gordon  (1970),  Jackson  and  Kerr  (1975)  and 
Gordon  and  Cato  (1972)  which  are  between  the  continuum 
sources . 

An  approach  taken  in  making  the  observations  would 
be  to  use  a  large  single  aperture  telescope  to  observe  a 
high  frequency  line;  say  the  H  92a  at  8.3  GHz.   The  ob- 
servations could  be  made  in  a  rectangular  pattern  of  con- 
tiguous beams.   The  extent  of  the  grid  should  be  thirty 
minutes  of  arc  in  longitude  if  possible  and  about  thirty 
minutes  above  and  below  the  plane.   Several  of  these  grids 
should  be  observed  along  the  plane.   One  of  the  grids 
should  encompass  the  region  observed  by  Gordon  and  Gottes- 
man (1970)  and  Altenhoff  e_t  aj_.  (1973).   One  grid  should 
be  placed  at  the  outer  boundary  of  the  ridge  and  at  least 
one  placed  at  an  intermediate  longitude.   The  outer  boundary 
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of  the  ridge  is  not  expected  to  be  discrete;  rather,  as 
the  longitude  increases  the  number  of  nebulae  decreases 
until  the  line  of  sight  density  drops  to  zero.   This  longi- 
tude, which  is  expected  to  be  between  40  and  45°,  can  be 
estimated  by  observing  individual  regions  between  the  con- 
tinuum sources  in  this  longitude  range.   At  the  ridge 
boundary  there  will  be  a  minimum  number  of  nebulae  along 
the  line  of  sight;  thus,  there  is  a  better  chance  that  one 
nebula  can  be  identified  individually. 

This  technique  should  provide  a  better  estimate  of 
size  and  distribution  —  radially  and  normal  to  the  plane  — 
of  the  nebulae.   It  would  also  provide  a  comparison  between 
the  ridge  in  general  and  the  region  in  which  Altenhoff  et 
al .  (1973)  have  observed  recombination  lines  from  discrete 
sources  and  the  hot  H  I  gas.   Unfortunately  the  recombination 
line  power  decreases  as  the  line  frequency  becomes  large. 
For  that  reason  each  beam  position  would  require  at  least 
7  hours  of  integration  time  in  order  to  get  a  AT  of  0.002°K 
if  the  43-meter  NRAO  telescope  is  used,  making  the  entire 
program  a  large  undertaking.   It  would  have  to  be  broken 
up  into  units  in  order  to  fit  in  with  telescope  scheduling. 

More  information  on  the  thermodynamic  condition  of 
the  gas  in  the  nebulae  can  be  obtained  by  observing  sev- 
eral higher  order  ( A n >  1  )  lines  all  with  frequencies  within 
a  small  range.   A  set  of  a-,  3-,  y-    and  higher  order  lines 
should  be  observed  from  several  regions  between  continuum 
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sources  along  the  plane  to  ascertain  if  the  thermodynamic 
state  varies  within  the  ridge.   From  each  line  the  temper- 
ature of  the  gas  is  determined  using  the  assumption  that 
LTE  applies.   As  An  increases,  the  excited  states  of  the 
atom  are  more  apt  to  be  col  1  i si onal ly  excited  and  the  meas- 
ured excitation  temperature  approaches  the  electron  kinetic 
temperature  of  the  gas.   The  amount  by  which  the  measured 
temperature  is  less  than  the  kinetic  temperature  is  a  meas- 
ure of  the  departure  from  LTE.   This  technique  was  applied 
to  denser  nebulae  by  Hjellming  and  Davies  (1970)  from  which 
the  line  powers  are  considerably  larger  than  they  are  from 
the  nebulae  under  consideration  in  this  model.   Consequently, 
the  highest  order  of  line  observed  will  be  smaller  and  the 
kinetic  temperature  may  have  to  be  estimated  by  extrapola- 
tion.  The  results  of  this  experiment  should  be  compared 
to  the  work  of  Matthews,  Pedlar,  and  Davies  (1973)  and 
Jackson  and  Kerr  (1971,  1975)  who  find  the  radio  recombina- 
tion line  emission  regions  to  be  close  to  LTE.   The  results 
should  also  be  compared  to  those  of  Brown  and  Balick  (1973) 
who  determine  that  this  region  should  be  considerably  out 
of  LTE. 

Once  the  size,  kinetic  temperature  and  thermodynamic 
condition  within  the  nebulae  are  better  understood,  it  is 
hoped  that  a  source  of  energy  can  be  confirmed.   A  possi- 
bility which  has  been  suggested  in  Chapters  V  and  VI  is 
that  the  energy  source  is  one  or  more  OB  stars.   Perhaps 
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if  the  distribution  of  these  nebulae  within  the  ridge  disc 
can  be  determined,  it  will  provide  some  confirmation  of 
this.   However,  the  distribution  in  the  galaxy  of  OB  stars 
at  distances  greater  than  2  kpc  is  not  known. 

Obtaining  a  better  understanding  of  the  relation 
between  the  radio  recombination  nebulae  and  cold  H  I  regions 
may  come  from  observing  with  a  matrix  of  beams  at  a  high 
frequency  line  with  the  matrix  centered  on  a  source  for 
which  Radhakri  shnan  et_  aj_.  (1972b)  or  Hughes,  Thompson,  and 
Colvin  (1971)  have  published  a  21  cm  absorption  spectrum. 
The  velocity  resolution  of  the  observation  must  be  as  good 
as  the  absorption  spectrum. 

Confirming  or  rejecting  a  relation  between  low  fre- 
quency absorption  and  radio  recombination  line  emission  in 
the  ridge  will  require  collaborative  observations  at  low 
frequencies.   Presently  the  80  MHz  optical  depths  of  Dulk 
and  Slee  (1972)  are  open  to  question  (cf.  Chapter  VI). 

As  mentioned  in  Chapter  VI,  the  distributions  of 
these  nebulae  can  not  be  reliably  determined  using  circular 
galactocentric  orbits.   A  density  wave  model  like  that  de- 
veloped by  Burton  (1971)  or  Burton  and  Shane  (1970)   would 
produce  a  more  correct  distribution,  but  the  nebulae  must 
be  positionally  related  to  21  cm  emission  regions.   Such 
a  relationship  will  be  difficult  to  develop.   If  the  energy 
source  of  the  nebulae  can  be  determined  —  for  example  if  it 
is  an  OB  association  —  the  type  of  source  may  suggest  a  dis- 
tri  buti  on . 
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The  purpose  of  this  study,  as  stated  in  Chapter  I, 
has  been  to  determine  how  the  intersteller  medium  in  the 
ridge  differs  from  that  in  other  regions  of  the  galaxy 
such  that  it  is  a  well-defined  source  of  radio  recombina- 
tion lines.   Although  considerable  uncertainty  still  exists 
there  is  a  trend— which  this  study  s upports— toward  the  ac- 
ceptance of  a  distribution  of  distinct  low  density  H  II 
regions.   As  this  chapter  shows,  there  is  much  work  that 
can  be  done  to  reduce  the  uncertainty.   The  location  and 
proposed  model  of  the  ridge  suggest  that  it  is  an  extreme 
population  I  region.   The  relation  between  galactic  struc- 
ture and  star  formation  has  been  shown;  cf.  Shu  (1973)  and 
Biermann  (1973).   It  follows  that  the  ridge  may  well  prove 
to  be  an  important  component  of  spiral  galaxies. 


APPENDIX 

The  following  are  details  of  Equation  1 1  -  7  and  the 
derivation  of  Equation  11-16. 

The  Thermal  Continuum  Absorption  Coefficient 

The  general  form  of  the  thermal  continuum  absorption 
coefficient  (Equation  1 1  -  7 )  in  terms  of  fundamental 
constants  is 
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where  Z  is  the  effective  nuclear  charge  number;  e  is  the 
electron  charge;  M  is  the  mass  of  the  electron;  k  is 
Boltzmann's  constant  and  y    is  Euler's  constant. 

Population  of  Excited  Levels 

■k 

This  is  the  derivation  of  Equation  11-16.   If  N  .  is 
the  total  number  of  ions  in  the  i  *-  h  ionization  state  (e.g., 
number  of  hydrogen  ions,  N^+),  then,  taking  Equation  11-15 
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and  summing  over  all  levels, 

co  oo  ■* 
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N.  can  be  found  from  the  ionization,  or  Saha, 
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N   is  the  neutral  atom  number  density  and  u(T)  is  the 
o  o 

partition  function  for  the  neutral  atom.   The  ionization 
energy,  x-»  ^s  tne  energy  required  to  remove  the  last 
electron  to  form  the  i    state  of  ionization.   Setting 
i  =  0  in  Equation  A-3  and  dividing  the  Saha  equation  by 
it  gives 
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N   is  the  quantity  desired;  so,  the  LTE  number  density  of 
no        ^ 

n-level  excited  atoms  is 
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All  number  densities  are  in  units  of  cm 

e .   =  - ( y  .  -  y.  )  is  the  threshold  energy  which  is 
in     Ai    Ai n  aJ 

the  total  energy  of  the  electron  in  its  orbit  and  equal 

to  the  sum  of  the  kinetic  and  potential  energies,  e„   +  V. 

For  hydrogenic  ions,  electron  shielding  of  the  nucleus 

can  be  ignored.   If  a   is  the  Bohr  radius  of  the  n    level, 
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The  orbital  kinetic  energy  is  eK  =  1/2M  aw  ;  to  is  the 

orbital  angular  frequency.   Then 
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Finally 
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Also,  e.  /k  =  157890  Z  /n   from  substitution  of  the  values 

of  the  fundamental  constants. 
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